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Doctor of Philosophy in Occanography 
University of California, San Diego, 1974 


EProressor Carl H. Gibson, Chairman 


Measurements of velocity, temperature and humidity were made 
in the atmospheric surface layer over the open ocean from the 
mesearch Platform FLIP, to compare estimates of the fluxes of 
momentum, sensible and latent heat, using the direct covariance, 
inertial and direct dissipation techniques. Measurements were made 
during a series of occanographic cruises approximately 60 miles off 
۱۱۶ coast of Baja California, between August 1970 and April 1973. 
rect covariance estimates of the fluxes were corrected for 
instrument tilt using simultancous measurements of pitch angle, and 
Gissipation estimates were corrected for stability effects. 


Agreement was obtained between estimates of momentum 


XXIV 





flux by the dissipation techniques and the direct measurements within 
the uncertainties of the direct estimates (+ 25%), applying the 
assumption that total production of kinetic energy (mechanical + 
buoyant) equals dissipation. From the direct dissipation technique, the 
average value determined for the velocity subrange constant was 

E =0.55 0.05. Trends in the direct dissipation estimates indicate 
that vertical turbulent transport and pressure transport of kinetic 
energy cannot be neglected and may be approximately half of the 
estimates of dimensionless transport measured over land. Comparisons 
Ere made between directly measured momentum fluxes corrected for 
instrument tilt using measured pitch angles, and values obtained using 
estimates of pitch angle in lieu of measurement. Results of the 


comparison suggest that large errors in direct momentum flux 









estimates can be incurred if tilt corrections are based on estimates of 
pitch angle. 

For agreement between inertial dissipation estimates and 
directly measured latent heat fluxes, an average value of the humidity 
Subrange constant = = 0.21 20,05 was required, “The direct 
dissipation technique was not applied to E latent heat fluxes since 
he frequency response of the Lyman-alpha humidiometer was not 
sufficient for direct estimates of the dissipation of humidity variance. 
imensional humidity spectra were similar to horizontal velocity 


pectra, and water vapor flux cospectra were similar to momentum 





flux Eospectra. 
| To. obtain agreement between direct estimates of sensible heat 
flux and inertial dissipation estimates, an average value of the 
temperature subrange constant Bo - 2.م‎ 2 0.3 was required. 
Estimates of sensible heat flux by the direct dissipation technique e 
migher than direct measurements by a factor of 2 or more. From 
direct estimates of the dissipation of temperature variance (+ 20%) 
and sensible heat flux (+ 10%), dissipation was greater than production 
of temperature variance by as much as a factor of 4, not accounted 
for by uncertainties in the di mee estimates, or by temperature flux 
divergence., Direct dissipation estimates of Bo varied from 0,7 + 
mmo? to 2.4 20.2. Temperature specira from different cruises with 
different environmental conditions did not compare with one another, 
and were not similar to humidity or velocity spectra. Cospectra of 
EEnsible heat flux exhibited similar differences. The differences 
ween production and dissipation of temperature variance, and the 
behavior of the temperature spectra and sensible heat flux cospectra 
could be accounted for by consideration of sources of production usually 
Beglected in the 7 variance budget, and attributed to the 
combined effects of ocean spray evaporation and radiative heating / / 
۲۰٢٠ 

A wide range of bulk coefficients were obtained from comparison 


Of estimates of sensible heat flux by the direct covariance and bulk 
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aerodynamic techniques, indicating the bulk technique is not reliable 
for estimates of sensible heat fluxes. A value of (1.29 + 0.36) x Jos 


ES obtained for the latent heat flux bulk coefficient. 
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fee PRODUC TION 


Een and moisture transfer from the oceans to the marine 
atmosphere are major driving faetors of atmospheric and oceanic 
eireulation. Within the last deeade a vast àmount of seientific 
investigation and experimentation has been directed toward an 
understanding of the energy exchanges and physieal processes which 
oceur at the atmosphere - oeean interface. Aeeurate determination of 
the vertical fluxes of momentum, sensible and latent heat is of prime 
EUportanee. This researeh has as its primary objeetives tine direet 

N 
measurement of the turbulent fluxes and eomparison of various 
teehniques of estimating the fluxes indireetly from related statistieal 
12111101 0 5 

Perhaps one of the earliest attempts to estimate the transfer of 
momentum and heat from the oeean was made by G. I. Taylor in 1913 
from a whaling ship using balloon and kite observations of mean 

E 1 
temperatures and wind veloeities (Taylor, 1970). Due to severe 
instrument limitations, early methods of estimating the turbulent 
fluxes have eonsisted of semi-empirieal theories relating mean 
properties or gradients to the aetual fluxes. No direct flux calculation 
was possible. Many of these methods including the bulk aerodynamie 
method, integral method, and profile method are diseussed by Roll 


Z 
(1965). ۲ 





With the cvolution and development of sophisticated instrumen- 
tation and experimental platforms for use on the open ocean, accurate 
measurements of atmospheric and oceanographie paramcters of scales 
ranging from millimeters to kilometers have become possible. As an 
outgrowth of the recommendations of the Joint Panel on Sea-Air 
۱ C 
Interaction of the National Academy of Sciences (1962, 1966) the 
Bardados Oecanographie Moteorological Experiment (BOMEX) was 
conducted during the summer of 1969. The principal objective in air- 
sea interaction was the study of the fluxes of momentum, sensible and 
E 5 
latent heat deseribed by Davidson (1968) and Kuettner and Holland 
6 
(1969). Current results from BOMEX have been summarized by 
Holland (1972). i Results of direct flux measurements from the 
research vessel FLIP obtained using the covariance or eddy correlation 
method and the well recognized difficulties of the measurement are 
8 7 
ENSscussed by Pond et al (1971), Phelps (1971), Phelps and Pond 
10 1 SR : ; 
ESI) and Leavitt (1973). The dissipation technigue for 
1 
estimating momentum flux was used by Gibson and Williams (1969) 
from FLIP in 1968, and for sensible heat flux during BOMEX by 
LS 7 : 
Gibson, Stegen and Williams (1970) and Steven, Qibson andi Friche 
14 
2973). 
Experimental data for this thesis was obtained from measure- 


ments made from the Research Platform FLIP during a series of 


oceanographic cruises betwcen August 1970 and April 1973 off the 





coast of Southern California. Mean and fluctuating component ° 
horizontal and vertical velocity, temperature, humidity, instr: ent 
Exotion, and mean sea surface temperature were measured 
simultaneously. Fluxes of momentum, sensible and latent heat were 
obtained by the direct eddy correlation technique corrected for 
instrument /platform motions and compared with estimates of fluxes by 
the direct and inertial dissipation techniques. The influence of 
Eunmospbheric stability on the determination of the fluxes by me latter 
two methods was determined, Latent and sensible heat fluxes were 
also determined using the bulk aerodynamic technique (including 

` 


stability effects) and compared to estimates of the turbulent heat fluxes 


by the direct covariance technique. 
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EN Direct Covariance Technique 

The optimum method of determining the turbulent fluxes of 
momentum and sensible and latent heat is to measure directly and 
appropriately average the covariances of vertical velocity with the 
suitable variable of interest. This method is known as the direct 
Ecvariance or eddy correlation technique. The technique is based 


on the definition of the vertical fluxes as given by 


Momentum Flux T = -pP uw 
Sensible Heat Flux Hg = PC, w 6 (1) 
: TL Heat Flux Hr = L wq 


where u, w are the horizontal and vertical "termes ene components 
of the streamwise velocity vector U , 6 is the fluctuating component 
Metemperature ("C), and q the absolute humidity fluctuation 

(u ۳٢ ہہ‎ The constants are p density, e specific heat 
capacity, and D, latent heat of vaporization. Overbars indicate 
timc averages assuming stationary flow. The covariances are 
obtained by integration of the appropriate cospectrum Re (where x 
is u, 8, ora) between low and high frequency limits in the same 
manner as the variances of the individual turbulent components are 
obtained from the appropriate spectrum T . The frequency range 


is determined to include all significant contributions to the integrals 





‘given by 


— — 


< 2 : 
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WX 1 ees ۳ ( ۰ ۰ x f gy ۳ی‎ (2) 


The low frequency ) or high pass cut-off is normally determined 
by removal of the record mean on data processing and is thereby set 
y the length of the data اسف‎ The low pass cut-off ڑگ‎ is 
initially limited by the bandwidth limitations of the instruments used 
to measure the parameters u, w, 0andq, and is subsequently set 
by low pass filtering during data analysis, usually adequate to 
determine the full value of thé appropriate flux. Filters are used to 
reduce high frequency noise and minimize aliasing in analog to digital 
conversion, In this work the frequency range of the integration was 
between about 1x w and JU Hz, 

Determination of the fluxes directly by this technique is fairly 
E cult in practice due to the effects of instrument platform motion 
on the measured turbulent velocity components, Pond et al (1971) 
obtained values of the covariance uw with the eddy correlation 
pEehnique using a sonic anemometer to measure horizontal and vertical 
Belocity from FLIP during Operation BOMEX. Employing an empirical 
merrection procedure, نس‎ corrections to uw were required to 
account for F LIP motion and effects on the flow field. Rotations of 
the principal axis of the Reynolds flux tensor were selected to make 


172 


the spectral correlation coefficient R () = @ /(@ 4 ) 
UW UW uu WW 





equal to a value of -0.5 for 0.01< iz/U< 0.1 as suggested by 


E 


Smith (1967 ans Weiler and Burling (1967). 7 Pond ect A found 


typical correction valucs for uw were 13% per degree tilt and 5% 
and 3% per Bocreo for wq and wô respectively, and 1t was الت‎ her 
the measurements were within 1-2° of the correet coordinate system. 
Montribuiions to the u, w spectra and uw cospectra by wave 
induced FLIP motion were not ineluded in the integration when 
determining variances and eovarianees, Despite the fact that no 
Eeoretical assumptions or approximations are required to employ the 

۱ 8 
rect covariance technique, Pond et al assumed that there was no 


distortion of the Reynolds stress due to the large tilt angles calculated 


(~ 10° ) which would affect the correlation between u and w. 





: 17 ۱ 
Kaimal and Haugen (1969) have found that the correlation 
b. uw | | 
coefficient r è a varies over a wide range from near zero to 
uw 
uw 


E0.5, particularly under unstable conditions. This suggests that the 
spectral correlation coefficient Rw may also vary with stability. 
Based on results from the Kansas experiment, Kaimal and Haugen 
suggest a need for + 0.1° accuracy in the internal alignment and 
mounting of sensors used to measure the direct covariance (uw) in 
the atmospheric boundary layer.’ 

For this research measurements of instrument platform tilt in 


the u, w plane were made using à vertical gyro located at a known 





angle with respect to the velocity sensors. Outputs of the vertical 
gyro were recorded simultaneously with the fluctuating velocity 
components and temperature and humidity fluctuations to determine 
corrected values of the covariances. Description of the instrument 
package and correction procedures are contained in the 


instrumentation and data reduction sections, 


2,2  Dissipation Techniques 

Both the inertial and direct dissipation techniques are somewhat 
advantageous over the direct covariance method in that they are much 
fess susceptible to errors "^ ed by instrument platforn? motion. This 
is due ta the fact that the flux estimates are based on measurements of 
the small scale structure of the velocity and scalar fields, Three 
variations of the dissipation method differ im the technique used to 
determine the rates of dissipation of velocity or scalar variance are 
determined, The dissipation may be measured directly or inferred 
from رت‎ subrange of the appropriate spectra or of the 
Erructure function. 

The technique for relating the dissipation of velocity and scalar 
variance to the fluxes of momentum and heat was first suggested by 


ES J. Taylor Eon S 


This technique is based on the assumptions 
1 ۱ 
discussed in Lumley and Panofsky (1964) ? that mechanical production 


is equal to energy dissipation at the same height in the constant flux 





Ever for near neutral conditions, so that 


— dU = 
(production) -uw ne € (dissipation) [29 


A 





>? 
* . ۰ ۰ 一 
Mere € = mean rate of viscous dissipation = 2ve,, , and e = 
11 paire 
1] 


tu, fe ; V = kinematic viscosity; p density, U mean‏ کی 
l, jJ J»‏ 
streamwise velocity, z = vertical height above the ae sea surface,‏ 


and u and w are defined as previously. 


The same assumption is made for the production and dissipation 


۱۳۲ ۱121 variance, i.c., 
— d | 
(production) -yw "+ ٦ (dissipation) (4) 


where y and y represent the mean and fluctuating scalars (tempera- 

ture (T- 8) and humidity (q+ Gas X, - mean rate of dissipation of 
2 

scalar variance = use "a and 2 , the molecular diffusivity 


e the scalar. 


Assuming a logarithmic profile for velocity and the scalars, and 
implicitly invoking conditions of stationarity and horizontal 


homogeneity, we find 


E. E 
E ^ = luy (5) 
dz Kz 

n. 7 ۱ wy 

一 1 Do = = (6)‏ كك 





Biere u, is the friction velocity, K is von Kärmän's constant, 

normally assumed tobe 0.4, y. the scalar scalo, and 

meee (yw OU/6z)/(uw ty /d2z) = K /K the ratio of eddy transfer 
1 H m i 

Coefficients for the scalars and momentum, assumed unity for neutral 

conditions (Reynolds analogy). Also assumed is K, ^ K 


H E 


humidity, Expressions for the fluxes can be obtained by combining 


for 


equations (3) and (5), and (4) and (6) 


2 EN 
a” = K€ xj"! (7) 
2 & X 2, 
y ہے سا ہے‎ (8) 
= Zu ^ 


2.2a Inertial-dissipation technique 
8 一 一 s 
BU or determined values of € and Xo from the inertial 


nee of the second order structure functions for velocity and 


temperature. From the similarity hypothesis of Kolmogorov (1941) ae 


2 
Ber turbulent velocity fields, and Obukhov (1949) i and Corrsin an 


for turbulent scalar fields assuming local isotropy, the second order 


۳۰۰ ۲٦٢٢۶٢ functions in the inertial subrange are given by 


— s A 928 
Ec o-luwmEr)-u(x)] = a € r (9) 


uu ٠ ME 


Dagf?) !]اد‎ I) = ba x 3 1 (10) 


where U and T (or humidity) are measured at streamwise separation 





Miistance r, between L (energy scale) and n (dissipation scale). 
۱ ۱ ۱ 3 = 1/4 
n is also known as the Kolmogorov scale, defined as (v /e ) 


Dues Of € and y 7 can also be determined in ihe same 


manner from the inertial subrange of the one-dimensional spectra 


given by 
$ (k) = Q ne cde [orem 
uu u 
El Sec 3 
© dx) = € 12 
7 ا‎ )12( 


Where o and 2 are universal constants assumed known, ‘and the 
radian wavenumber k= 2m ۶۸0 by Taylor's hypothesis. This 
technique requires accurate values of the constants o and de [on 
a and ارر دا‎ , and the velocity and scalar spectra and structure 


functions to follow the LA E 


form of the inertia! subranges, 

In this work, examples of temperature spectra obtained from 
meen Ocean atmospheric data are presented which appear to exhibit a 
very limited, if any, -5/3 inertial subrange. Application of Eq.(12) 
to single frequency (wavenumber) measurements of the spectral 
tion could lead to large errors in the estimation of sensible heat 
flux if an inertial subrange does not exist, Similar care must be taken 
when applying Eq. (10) to single separation measurements of the 


Berueture function. A substantial portion of the temperature spectrum 


Should be measured (to ~ 100 Hz) to determine if a subrange does exist 








before this method is applied. Humidity spectra appear to exhibit a 
-5/3 inertial subrange behavior over a wider range of frequencies 


Eo c 20 Hz) than temperature 


END Direct dissipation technique 

This technique involves measurement of the mean square 

values of time derivatives of fluctuating streamwise velocity and 
temperature to calculate dis sipation ۶۹٢٢٢ 013 بت‎ ء۱٣۳۰۰٠٠۹‎ Durecteostiniotes 
of € and m in the atmospheric boundary layer over the ocean from 
FLIP in 1968 and during Project BOMEX are described by Gibson and 
E. II 3 E 14 
Williams (1969), Gibson etal (1970) ^ and Stegen et al (1973), 


Assuming local isotropy and Taylor's hypothesis (U d/dx = -d/dt), 


the viscous dissipation may be determined from the relations 











٩ 2 
€. Y ES (13) 
m | 
and 
چ ب‎ 03) 
= : 


ENDr2Sctce the appropriate time derivative spectra are integrated to 
Exam values of the mean square time derivatives, Difficulties in this 
technique involve the high spatial resolution required of the sensing 
probes to nearly the Kolmogorov scale (17), of the order of 1 mm 


in atmospheric flows. High frequency response (~ 2000 Hz) and very 





high signal-to-noise ratios are required in the sensor and associated 
circuitry. At the present writing a humidity sonsor with sufficient 
frequency response was not available to measure adequately mean 
square 5005 time derivatives to calculate the humidity dissipation. 
feo a result, the direct dissipation technique is applied only to 
estimates of the momentum and sensible heat fluxes. 

The direct dis E technique does offer several 
advantages over the other techniques. Estimates of the fluxes require 
relatively simple instrumentation and data analysis, and are not 
affected by sensor motion betause of the high frequency response and 
fine scale spatial resolution of the sensors used, Measurement of 
۱۳ ۰ Over the entire range of frequencies (wavenumbers) are 
readily obtainable to determine spectral shapes. The spectra 
combined with dissipation measurements provide a means of accurately 
determining the universal constants. 

2.3 Stability Considerations 

Application of the dissipation techniques using Eqs, (7) and (8) . 
does not account for effects of stability on the flux estimates. These 
equations are modified by consideration of the turbulent budget 
equations for kinetic energy and scalar variance, and the flux-profile 
relationships in the surface layer, within the framework of the Monin- 
Obukhov surface-layer similarity theory (Obukhov, 1946 and Monin- 


Obukhov, 1954). = 





In the first tens of meters of the atmospheric boundary layer, 
phe fluxes of momentum and heat are assumed constant. Recent 
experimeníal support for this assumption has been given by Haugen 
et al (1971)? and Dyer and Hicks (1972). a In this "constant flux 
myer, according to similarity theory, the turbulence structure is 
determined by the surface shear stress >» the surface hest Wu Hes 
the buoyancy parameter g/T, and z the vertical height. The latent 
meat flux E, is included to account for humidity efíects 3 buoyancy 
in the Puna layer over the ocean. From these parameters 
Sharacteristic scales for velocity (u,.) , and scalars (d = T,,4,) 
are defined (as in Eqs. (5) and (6)). The length scales are z and the 
Monin-Obukhov length,defined as 

e T 


AT 


Dec یہ‎ (15) 
Kgw 


Mere the mean and fluctuating virtual temperature, T, and 6 y> as 
1 
given by Lumley and Panofsky (1964), 0 include humidity contributions 


to buoyancy. They are written as 


T. FEIL OO OLEM) (16) 
and 


0, * 040.61 Tm (17) 


where M and m are moan and fluctuating specific humidity (gm/gm). 


The scale L is a key independent variable of similarity theory and 





determines the eme ss of the surface Jayer above which buoyancy 
factors are not significant. 
Applying the scaling parameters to the vertical gradients of 
velocity and the scalars Gi ue. humidity}, Eqs. (5) and (6) are 


rewritten as 


E u 

dU x 

o RA‏ ا 

dz Kz SO UL) (18) 
= Dis 

d ^ ax 

E AL 

dz Kza, gr’) (19) 


The functions © , the dimensionless wind shear, and MES ۱ 
n 
the dimensionless scalar gradient, ardfewaluated empirically. It is 
generally assumed that 9 u = e for temperature and humidity. 
27 : 
Monin and Yaglom (1971) conclude that the forms of v and Sy 


for unstable conditions (the prevalent condition over the open ocean) 


Ere best given by 


中 = (1-0 نت‎ for ISL O (20) 
m m 

۲ p on 
Py = (1 - 0, 2/D) (21) 


e 


A review of the constants O and Ou has recently been given by 
m 


28 
Busch et al (1973) | and fairly good agreement appears to exist 


一 


یه 


for the value of 2 . Businger et al (1971)2 3 suggests O = 15, 
n 





P | P0050 | 
while Paulson (1970), Bacdpleyetal (1972), en and Dyer and Hieks 
(1970) 6 do ت‎ 0. 6 Miyake et al (1970) = = 16 
7 m. H Dion a PS " 006 71 i 
8 10 u 
miso used by Pond ctal, Phelps and Pond, and Leavitt, to obtain 
results from BOMEX data. For neutral conditions (s = 中 u^ J ۲ 
m 
۱ e9 
۱3۹1۳۵61۲ et al (1971) find values of = = 1,35 ander UNES 


Eomparced to the often used values of Po = 1.0 and K = 0.4, and 


5 0 = : 
sugges H 9 


۲ 4 Bud get Equati ons 


If horizontal homogeneity and stationarity are assumed and 
N 


the effects of humidity on buoyaney are ineluded, the budget 


2 2 2 2 








equation for kinetie 6۰" Per unit mass, e =u tw Tvinna 
« 4 4-4- 1 ES 
be written as (Lumley and Panofsky, 1964) 
۳۰ کس‎ — — 0 2 IEEE ہو‎ 
۴ ۷ E E 9 £0.61 T wm] -5 = we DNE pw-€ =0 (22) 


= 
The non-dimensional form of the budget equation is obtained by 
multiplying Eq. (22) by K fa; and employing Eqs.(15) and (18) so 
that Eq. (22) ean be written as 


LEM i =. Dea O (23) 
m 6 p € 


۰ 


The terms represent the turbulent shear production, buoyant 
production, the flux divergenee of turbulent kinetic energy, the flux 


divergence of pressure transport and viscous dissipation 





جح۰ ۳ 

۱۱۱۱ approaches to sinplification of Egs. (22) and (23) have 
meen suggested, The first approach, given by Eq. (13), is the 
assumption that production and dissipation of mechanical energy are 
Euual for near neutral conditions. For Iz/L | < US 5115 لت‎ 1 

| ۰ ۰ ۰ ۰ 34 م / 

Panofsky (1968) suggest that dissipation is balanced by the sum of 
buoyant and mechanical production, and that the flux divergence terms 


Ens be neglected so that Eq. (23) simplifies to 


ou 72۔-۔‎ = Q : (24) 
m € 


S 
* 


McBean, et al (1971)? conclude from measurements over land that for 
near neutral conditions (in « 0.2) Ed. (24) is a good approxima- 
tion. However, for more unstable conditions, - z/L>03, 
dissipation appears to exceed production, and for z/L-^ -0.5 the sum 
E ھت‎ > o is about 30% of shear production and more than half 

of buoyant production. Miyake et al (1970) » and Pond ctal (1921)? 
apply Eq. (24) to calculate momentum flux, however, they do not include 
stability effects on the mean velocity profile as given by Eq. (18). Hicks 
and Dyer (1972)°° include buoyancy effects from Eqs, (18) and (24) to 


calculate the momentum flux by the inertial dissipation technique. 


A, recent approach to simplification of the turbulent energy budgct 


| S 
has been suggested by Wyngaard and Cote (1971) based upon 





measurements in the surface layer over land. From measurements of 





7 ۰ 
we at three different heights the turbulent flux divergence of e 8 was 
estimated and found to approximately balance buoyant production in 


the stability range - 1.0 5 z/L * 0,  Eq.(23) reduces to 


N 


A 7 25 
= (25) 


where the imbalance I, is attributed to pressure transport 

divergence. Under stable conditions the imbalance may not be 
substantial and under strongly unstable conditions (4/L = - 2) 
dissipation and total production did balance. Empirically, the 


expression for the imbalance term is written as 


/4 ٤ 


3 
I = (1-15 TE - )1 + 0, ظ‎ | 27/3, (26) 


a | 2 
Gl 


Mere the first term on the right-hand side is the Businger formula 

fOr dimensionless shear production (Eq. (20), oc - 15), and the 

second term is the empirical expression determined by Wyngaard and 
7 تہ‎ ۱ ۱ 

Bote for unstable dissipation data. Using this approach, 


momentum flux corrected for stability effects may bc calculated using 


the empirical expression for unstable dissipation written as 


2 





E 
2 
ET 


EOS IL (27) 


— 





EIS expression for the dissipation was employed by Stegen et al 
14 ۱ M 
EU) to calculate momentum flux by the direct dissipation 
technique. 
In recent work employing data from Operation BOMEX, Leavitt 
1 : : BER: 
(1973) found that total production approximately balanced dissipation 
in the stability range - 0.22 z/L* - 1.5, and that dimensionless 
turbulent transport as estimated from measurements at two heighis was 
only about half of the dimensionless buoyant production as measured by 


37 
Wyngaard and Cote (1971). 


Z 
In general the budget of temperature variance (0 ) ray be 


: à 
written (Monin and Yaglom, 1971) í 
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where repeated indices are summed and u, and x. may DE replaced 
by (u, v, w) and (x,y, z) respectively. Here © is the mean 
Potential temperature. The variance budget can be 
simplified by invoking the assumptions of stationarity and horizontal 


homogeneity; Eq. (28) can then be written as 


== b li ae 
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where the dissipation of temperature variance x o is defined as 





و ۲ 
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SR 2 
Ina similar manner, the budget of humidity varianca (q ) may be 


۱ / 2 
written as (Coantic and Leducq, 1969) : 


on 
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ot Ss SEV. us Dre Os Ss q dx, 
1 J J J 
۲ 2 
2D E + و‎ (31) 
q gud g 


which simplifies under the same assumptions of stationarity and 


Borizontal homogeneity to 
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۳ wq E SK ano (32) 
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and the dissipation of humidity variance is defined as 





= one A) ۱ (33) 


The first three terms of the simplified budget equations for tempera- 
Eure and humidity variance (29) and (32), represent production of 
scalar variance, transport of scalar flux divergence, and scalar 


variance dissipation due to the smoothing action of molecular 





diffusivity respectively. The temperature variance budget contains 


additional source or sink terms not included in the humidity budget, 





The term 20R represents the effec:s of heat transfer due to radiation 


E 
C 

P 
the temperature field due to the internal friction of the fluid, and S 





تد 


on the temperature variance, the effects of viseous dissipation on 
and Bu Eepresent generalized souree or sink terms to account for any 
phase changes of water vapor or interactions of the temperature and 


humidity fields due to radiation. 


2 3 
Monin and Yaglom (1971) 1 and Plate (1971) ? suggest «by 
hypothesis that the viscous internal heating term 20€/C is insignificant 
P 
40 
End may be neglected. Friehe (1973) suggests that this term may not 
Be negligible and that there may be a positive correlation of 6€ due to 


an increase in 6 produced by local, rapid adiabatic heating by viseous 


dissipation. Although a mean value for viscous dissipation € may 


2 


: ۱ ۱ d 
be estimated from the mean square of the velocity derivative sp 
2 


1 0 
(see Eq. (13)), à positive correlation of O e may not aeeurately 
reflect the actual correlation between fluctuating temperature and 
2 


the instantaneous viscous dissipation 8e; کے‎ is only one component 


of dissipation, used to ealculate € by the assumption of isotropy. 


In general, the effect of radiative heat transfer is to reduce the 


4] 
value of the mean square temperature fluetuation (Goody (1956), 


3 
Townsend (1958), L Plate (1971) 7 For internal radiative transfer 








ہس —— 


A ۱ 2 
۲ ۵۱۷۸5 010۵ indicates that maximum destruction of 6 occurs at 
very small scales of the temperature fluctuations since the most 
pronounced radiative transfcr occurs ovcr the shortcst path interval. 
For atmospheric conditions with high humidity content and strong 
جح‎ 0 
radiative transfer, Phelps and Pond (1971) suggcst that more 
absorption of long wavc radiation occurs at larger schlcs resulting in 
a more homogeneous large-scale tcmpcrature field and suppression of 
39 a 
low frequency temperaturc fluctuations. Plate ' also indicates that 
partial absorption of long-wave radiation by atmospheric water vapor 
could causc a change of statc, (conversion to sensible hcat) which 
would influence the tempcrature field. It is not clear whethcr this 
absorption and change of state of water vapor would act as a source or 
sink in thc temperaturc and humidity variance budgets, and could be 


۲۰۰٠۱١٦٤٤ i thc generalized source terms. 


Coantic and Seguin (197 iy have shown that radiative flux 
divergence duc to absorption of infarcd radiation (long-wave) by watcr 
vapor and carbon dioxide always acts in such a way to increase the 
absolutc value of the turbulcnt heat flux. They conclude that the constant 
heat flux hypothesis for the surface layer may not be valid, cspecially 
for conditions of moderate wind and high moisture content, even for 


heights above 1l meter. 


The gencralized sourcc term may also include the effccts of ocean 





spray on the temperaturc and humidity fields and possible convective 
conditions in the upper boundary of the surface layer as suggested by 


10 
helps and Pond (1971). 


To further simplify the budget equations of temperature and 
humidity variance it has been common practice (Miyake et a1 (1970), 2 
| 36 8 f 10 
Hicks and Dyer (1972), Pond et al (1971), Phelps and Pend (1971), 
/ PT 13 | 
Wyngaard and Cote (1971), Stegen et al (1972) | and Paulson et al 
44 
(1972) ) to assume production of scalar variance is equal to the 
dissipation of scalar variance (Eq. (4)). Neglecting source terms 
۹ 


in the temperature variance budget and assuming the transport terms 


may also be neglected نو را‎ 5 and (52) reduce to 


e dz 2 Ag en 
and 
En dq E 
o یمس‎ == Lo / 35 
BET 92 ۵ (55) 


18 
R. J. Taylor (1961) justifies the assumption of production of 


temperature variance equal to the dissipation based on negligible values 





2 
of the flux divergence term "Is pw © and the viscous heating term 
P 


E 45 
€ as compared to the production term.  Panofsky (1969) ^ also 
Suggests omission of the transport term in the temperature variance 


budget and that in general, this term is of the order of 10% of the 


Eduction or dissipation terms. Similar results were obtained by 
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Wyngaard and Cote (1971) based on measurements of the transport 

7 
term over land. Wyngaard suggests that the transport term is an 
order of magnitude smaller than the production for moderately 
7 E ۱ H 
unstable conditions -1s z/L< 0. Recent studies by Leavitt (1973) 
and results presented in this work indicate that in thẹ atmospheric 
boundary layer over the ocean, dissipation of temperature variance 
can greatly exceed production for moderately unstable conditions. 

e ا‎ 
Leavılt finds that the mean value of turbulent transport in the 
temperature variance budget is about 25% oi the production 
(representing a loss) for . 2.۰۲ z/L 2.75 +3۷۳ 
7 ۱ En 

Wyngaard's ,suggesting that for highly unstable conditions the 

46 
transport may become a gain, supporting Deardorff's (1966) 


suggestion that for a counter-gradient of potential temperature, 


positive heat flux can be supported due to turbulent transport. 





3. INSTRUMENTATION AND DATA ANALYSIS 


E Experimental Arrangement 


Data for this thesis was obtained from measurements made 
from the Research Platform FLIP during a series of oceanographic 
expeditions in August 1970 (SOMA), September 1971 (MITOS D, 
EEbruary 1972 (MITOS II), April 1۰۶ ۱۱۸1 و111 دلب‎ RE ۹ 
(OWAX). For cach cruise F LIP was towed to an operating Area) 
the vicinity of 31°40° North latitude, 118° 00’ West longitude; 
approximately 60 miles west of Ensenada, Baja California. After a 
stable vertical position had been attained, F LIP drifted freely with 
ecel le upwind. Wind direction during each cruise was 
predominately out of the north and northwest so that open ocean 
Eonditions were achieved. 

ihe vecontiguration of F LIP employed for the MITOs III and 
OWAX cruises is shown in figure 1. In this configuration an 
instrument package was located at the end of the 18. 3 meter port 
boom. The 10.7 meter starboard boom was used to mount a fixed 
and a trimmable sail to maintain F LIP in an upwind orientation and 
compensate for the drag on the port boom. The instrument package 
used for measurements of turbulent velocity, temperature, humidity, 
and instrument motion is shown in figure 2. It was possible to mount 


fae package rigidly at the boom height of 12.5 meters or to heights 
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Research Platform FLIP in the stable vertical position, 
with 18 meter port boom and compensating sails on 
starboard boom 
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FIGURE 2 
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Instrument package for flux measurements 


JR 7۳۲۶ ۶ ۲773۳۳۲۳۳7۳۲۹۹۲۲ 2۳۲ مہ‎ CMV TT EE TOY مجو عدم اجرج‎ RATA VIA وهار‎ BO OSITO سبج رجنب‎ OEY FP ڈو‎ OY رج جرب‎ OPTS 8 — sits nn M nd 3 


= 


me 
Ü oan tetra wae 0n Van rt وش جد ونث د ماداد‎ | 
١ 
t 
1 


۷ wipe tpa q.i y eratis nm an cma چو‎ 
| 
a 2 ) 


| 
ea fl 


: —À { 


& sam mes 
یر ہبی‎ 


حم 
si ۶‏ 
* 
p‏ 


Ab! 
: 1 ] ۱ j ۱ 
A dí tad athe د‎ > ba m مھ‎ Ay 


e! Ra Sas he tn en ei 
rae 
e lr cl rt m. عأ‎ 


3 p 
- “were 
وم‎ RAD) €^ 3 Irea miid 
re a ^ — ے‎ 
OA "Cv و وهی‎ e } 
مم‎ 
i = et 


^ AT qs mart a y 





۰ ٠ 
e E Sg a ۱ ڈیا‎ 
— PON Une 1 } 
۱ | j ¢ po re nomm nn. 
and ۳ T re 
LU e | li er 
. ， 5 ۱ نے‎ 
Van 
E 
* | 
i i y E- ۲ | 
| ne WS 7 ۱ 
۱ 7 
tae [ 
: ۹ | f 5 
A .. NW 
۰ Ner í ; 
٠ پا‎ | 
ھ0‎ | 
| i | 
i ^ 
SS | | 
1 $ 
| ) i 
y à 
0 0 | ! 
0 1 
۱ i i 
‫۰ . 1 
ET 二 : 





in Boe ee ell Ba a سو ی‎ ET ل‎ ACI ds a ae چھہا‎ == 
15. 5. 341 


"EE 





of 3 meters by means of a rigid extension pole conneeted to the end 

of the boom. Relative positions of the sensors on the instrument 
paekage for the MITOS III and OWAX eruises is shown in figure 2. 

ME sonie anemometer used for measurements of three veloeity eom- 
ponents and shown mounted to the instrument paekage in figure 2 was 
employed only for the MITOS III and OWAX eruises. Other veloeity 
sensors ineluded TSI or DISA hot wire and hot film anemometers for 
high frequeney measurements, Cup anemometers were used to 

E sore mean veloeity on the instrument package. High frequeney 
temperature measurements were made using thermistors and eold 
Wire probes, A Lyman-Alpha humidiometer was employed for small 
Beale S measurements, Preeision m 1161:10۵1 9 ۱ OMAN 
Ehermometers were used to measure mean air and sea temperature 
and wet bulb temperature for ealculation of mean absolute humidity. 
All instruments on the package were mounted rigidly and aligned in the 
laboratory so that any tilt sensed by the gyro would represent 
Bistrument tilt. 

During the February 1972 eruise (MITOS II) a failure of the 
large boom neeessitated removing the sails from the 10. 7 meter 
boom and using the shorter boom on the starboard side to mount the 
instrument paekage approximately 8 meters from FLIP. Due to the 
proximity of the Poni package to F LIP during the MITOS Il 


-一 


eruise and rearrangement of sensors, paekage motion was not 





measured and only a limited amount of high Írequency temperature, 
velocity and humidity data were available for analysis. 

The experimental configuration of FLIP for the first two cruises 
ol SOMA and MITOS I consisted of two 10.7 meter booms located on 
each side of FLIP. An instrument package similar to that shown in 
figure 2 except of lighter construction was attached to a rigid cable at 
the end of the port boom at heights between 3 and 8 meters. The 
package was stabilized by guy lines to minimize vibration and maintain 
mistrument orientation into the wind, The cable to which the package 
was secured was kept rigid by mcans of a weight attached at the end of 
fie cable and suspended 20 meters below the sea surface. Due to the 
Experimental configuration of FLIP, the location of the instrument 
package was closer to F LIP for the first three cruises than for the last 
two cruises. No pitch corrections were applied to direct covariance 
measurements from the SOMA and MITOS I cruises., It was still 
possible, however, to obtain high frequency velocity, temperature and 
humidity data for analysis and calculation of turbulent fluxes by the - 
dissipation methods. Data from the MITOS IIl and OWAX cruises 


Meremsed for tbe comparison of different methods of calculating the 


turbulent fluxes of momentum, and sensible and latent heat using pitch 


corrected direct covariances as a basis for the comparisons. 
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3.4 Data Acquisition 

The data acquisition and analysis system developed over the 
od ol the FLIP cruises is shown schematically in figure 3. This 
fipure 1s originally from Gibson (1973). zu For cach subsequent cruiso 
the sensors and instrumentation were slightly different with changes 
and improvements made based upon the knowledge ga ned from the 
experience and results of previous cruises. The basic system frame- 
work, however, remained the same. 

Calibration and test facilities included the wind tunnel of the 
Applied Mechanics and مو ٹس‎ Sciences Department and a TSI 
Portable, calibration air jet for velocity sensor calibration under 
laboratory conditions and on FLIP at sea. The cup anemometers and 
Sonic anemometer were calibrated in the laboratory wind tunnel before 
and after each cruise while the hot wire and hot film calibrations were 
performed on FLIP. The Lyman-Alpha humidiometer was calibrated 
im a Cambridge esters humaüdity test chamber. Thermistor and cold 
mare sensors were calibrated on FLIP using a calibration feature built 
i uhe temperature bridge circuitry. Temperature sensor calibration 
consisted of changing the current through the sensor in a balanced 
Bridse by an amount which caused the voltage across the sensor to 
change by the same amount as that resulting from a 1° C temperature 
change. Prior to and following the recording of data signals a series of 


calibration signals were recorded for voltage and frequency references. 








FIGURE 3 
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Analog instrumentation schematic 





Analog Instrumentation Schematic 
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3.¢a Temperature measurements 

Temperature fluctuation measurements were made using 
thermistors and platinum wire sensors (cold wire technique). For 
atmospheric eonditions of very low intensity temperature fluetuations 
(SOMA. cruise) thermistors were ideal since they have high sensitivity 
path a coefficient of resistance of - 3,2 (2/9/? C. The plantinum wire 
resistance probes have a coefficient of resistanee of 0.3 (/0/? C, 
The thermistors used for the temperature measurements are 1000 ohm 
Veco micro-bead thermistors manufactured by Victory Engineering 
Corporation of Springfield, New Jersey. The thermistors were of 
0.33 and 0,064 mm diameter with a flat frequency response out to 
ana 5 Hz respectively. The platinum resistance wire sensors used 
for fine structure temperature measurements were of 1mm length or 
less, 0.62 u diameter, with a frequency response to 1.5 KHz. The 
Cold wire sensors were made by Steve McConnell at UCSD employing 
mercially available TSI probe supports and etched fine platinum 
Eure. For the first FLIP cruise thermistors were used in a Tektronix 
BC66 AC bridge cireuit. For subsequent cruises the thermistors and 
cold wires were used in AC or DC (MITOS I only) Wheatstone bridge 
Circuits incorporating a low noise and high gain (80 db) amplifier for 
amplification of the bridge "error" signal. The bridge amplifier 
circuits were designed and built by Tom Deaton at UCSD. Equivalent 


resistance noise level for the resistance circuits was kept below 





0.001 C by employing impedence matching E coupling for 
the bridge and amplifier sections of the circuit for the low resistance 
platinum probes (between 30 and 150 ohms). The entire bridge 
amplifier circuit was powered by batteries and made small enough to be 
located physically close to the sensor, reducing the possibility of noise 
pickup in the cable between sensor and bridge. The short cable length 
between the sensor and bridge circuits also minimized imbalances in 
Ese bridge caused by cable capacitance. By incorporating a cable 
ESpacitance balancing circuit in the bridge circuit, any imbalance 
caused by the cable capacitance was removed from the amplifier, 

S 
leaving full dynamic range for the temperature signal. The sensor 
calibratión procedure described earlier consisted of changing the wire 
or thermistor current by 0.3% and 3,2% respectively, corresponding 
EE ue coefficient of resistance for each sensor. The calibration signal 
output from the bridge was a square wave of approximately 40 Hz with 


peak-to-peak voltage corresponding to 1°C. 


5. 2l Humidity measurements 


Humidity measurements were made with Lyman-Alpha 
humidiometers purchased from Electromagnetic Research Corporation, 
College Park, Maryland. The principle of operation of the humidi- 
Ometer is based on the strong absorption of Lyman-Alpha radiation 
N? 15. 6A) by water vapor. The radiation is generated oe a hydrogen 


Source tube, crosses the measuring path of approximately 1.5 cm 


fength and enters the detector, a nitrous oxide ionization chamber, 
The source tube and detector were equipped with special uv windows 
of magnesium floride for transmission of the sensing beam across the 
measuring path. These windows were much Jess susceptible to 
moisture degradation than the standard lithium flouride windows. Pond 
et al (1971) and Phelps and Pond 11971) ^9 describe successful 
operation of the Lyman- Alpha humidiometer Gur ine Operation ل‎ 7 
for humidity fluctuation measurements. During the BOMEX measure- 
ments a mechanical flow through path described by Phelps, Bond and 
Borner A670)” was employed. For this research no mechanical 
additions were made to the sensor and the measuring path was open to 
the mean airflow. 

Two different models of the Lyman- Alpha humidiometer 
Ee used during the FLIP cruises for this research. The BR model 
Es employed for the first two cruises and the BLR model for the 
Ebsequent cruiscs. The output voltage of the BR model sensor has the 
Form 
S qr 


V = V 
O 


(36) 
where À is a constant depending on the absorption coefficient for water 
vapor and the source currcnt, q the absolute humidity, and 4 the 
path length. For humidity measurements with this model the fluctuation 


calibration was estimated by linearizing the logarithmic relation of i55 


in the same manner as was done for the BOMEX measurements by 





Pond et al (1971)? so that 


a Em c کے‎ where me, + y (37) 


This method of calibration is probably good to within 10-15% since the 
fluctuations were very small and the mean humidity during the measure- 
ments was in a nearly linear portion of the instrument calibration curve 
M5" C dew point). 

The more recent BLR model humidiometer has a logarithmic 
amplifier with a 40 db range following the detector amplifier. The 
exponential variation of the detector output with water vapor is thus 
converted to a substantially linear variation of output voltage with water 
ESpor concentration. In this model an increased concentration of water 
apor gave an increase in output voltage while the opposite was the case 
for the BR model. Although some source tube aging was noticed, the 
fluctuation calibration for both models remained fairly constant, Any 
drift on the instrument was minimized by insuring that the transmission 
windows were cleaned before each run and the zero point (dry air point 
Bor Te BLR model) was set. 

The humidiometers were calibrated in a Cambridge Systems 
hygrometric test chamber located at the Naval Electronics Laboratory 
in San Diego. Due to the limited availability of the chamber only a few 
absolute calibrations of the sensors were performed. The BR model 


was calibrated in July and August of 1970 giving a fluctuation 





libration of 1l. 7 eem volt used for data from the SOMA eruise 
ےد‎ 1970( and the MITOS I eruise (September 197 ۱ 9 E 
model sensor was ealibrated in Mareh 1972 between the MITOS II and 
MITOS II c' se and again in February and April, 1973 before and 
Aiter the OVX cruise, The calibration results for these three tests 
were 4.70 movet 4.81 e and 4. 75 ne /volt 
respectively. These results also support the eonelusion of Pond ct A 
for the BOMEX results that the fluetuation ealibrations were fairly 
independent of source tube aging for souree tube eurrents which were 
used (> 0.25 ma). Absolute fluctuation ealibration of the E 
for the last three eruises USER. better than 10%. The frequency 
response of the humidiometer appeared'to be flat out to about 5-10 Hz. 
Measurements of mean air and sea surfaee temperature 
were made with preeision Hg thermometers and Hewlett-Packard Model 
2801A Quartz thermometers. The quartz thermometers were employed 
only in the MITOS Ul and OWAX eruises. Mean temperature data from 
he MITOS IH eruise indicated that the exposed quartz thermometers 
vere subjeet to radiation heating. For the OWAX cruise the thermom- 
ters were equipped with aspirated spirated radiation shielding. Absolute 
umidity was determined from measurements of dry and wet bulb air 
emperature using a sling psyehrometer. Dew Point temperature was 


lso measured with a Cambridge Systems dew point hygrometer. 
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3. 4c Velocity measurements 

Fine scale velocity measurements were made with 
linearized constant temperature anemometers and hot-wire, x-wire or 
hot-film sensors. The anemometers were purchased [rom Thermo 
Systems Incorporated (Model 1054) and DISA ۲س‎ ٢ (Model 55M01). 
The straight hot-wire, x-wire and hot-film sensors w ere also 
purchased from TSI or DISA, constructed of platinum or tunê and 
from 0.4mm to 1.25 mm in length with frequency response out to 
1,5 KHz or better. "The hot-film sensors were constructed of a 
E025 mm quartz cylinder ovér which a 0.5 mm long platinum film 
was plated with gold leads to the film. 

Düring the first threè FLIP cruises, precise measurements 
of horizontal and vertical velocity with the x-wire sensors often were 
not obtained due to 009 (~ 30% for MITOS I) drift of the probe 
calibration in the wind field, The primary source of the calibration 
drift appeared to be contamination of the hot-wires by salt spray. Use 
of the hot-film sensors during the last two cruises greatly improved 
the fine scale velocity measurements and sensor calibration drift was 
not a significant factor. 

A three-component sonic anemometer was also used for 
melocity measurements during the last two cruises. For the MITOS III 
Cruise (April 1972) an EG&G sonic anemometer was borrowed from 


Rome Air Development Center for velocity measurements over a hand- 
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width from DC to about 10 Hz. Measurement of velocity with this 
instrument is based on the principle of Doppler-shift of sound waves by 
the flow field in the measuring path. Since the instrument is an 
absolute sensor no calibration was required on FLIP, Unfortunately, 
due to a failure of a portion of the EG&G associated read-out circuitry 
only a limited amount of horizontal and vertical ات‎ data was 
obtained during the cruise. ror the OWAX cruise, a similar EG&G 
three-component sonic anemometer was purchased, and the sonic 
anemometer readout system was designed and built by Tom Deaton 

and Mark Barnum at UCSD, This instrument and readout circuitry 
were thoroughly tested and calibrated in the 30 inch square, low 
velocity wind tunnel at UCSD up to a maximum wind speed of 20 meters/ 
sec, and was successfully used during the OWAX FLIP cruise. Some 
minor distortion of the vertical sensing path of the sonic anemometer 
Was discovered during the testing and calibration procedures, This 
distortion was due to a slight misalignment of about 1.1? between 

the vertical sensors but was corrected for in the calculation of 
momentum and heat fluxes during subsequent data analysis. 

Mean velocity during the cruises was also measured with 
Teledyne-Geotech precision cup anemometers and associated readout 
Circuitry designed and built at UCSD. The cup anemometers were also 

alibrated in the UCSD wind tunnel and used for comparison to the 


ean velocity obtained from the sonic, x-wire and hot-film sensors, 
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p 2d Analog instrumentation 


The analog instrumentation employed during the F LIP 
eruises and for data analysis consisted of analop monitoring, 
mea@eessing, record and playback equipment as shown in figure 3. Data 
signals were conditioned for FM analog reeording by means of ampli- 
fieation and filtering to aehieve desired signa] levels, Several 
techniques were employed to improve signal-to-noise ratios and 
optimize the dynamic range of the tape reeorders. These ineluded the 
use of pre-emphasis (de-emphasis for playbaek) eireuits designed and 

^ 
built by Tom Deaton, and Tektronix 3A8 operational amplifiers as 
differentiators for analogue pre-whitening of fine scale veloeity and 
perature signals. Buck and gain amplifers designed by 
Bruee Williams at UCSD were used for removing known DC levels and 
Subsequent amplification of the bueked signals. The 3A8's were also 
used for signal amplifieation. Burr-Brown sum and differenee 
Bimplifiers were used to obtain horizontal and vertical veloeity data 
from x-wire probes, and also for removing tape reeorder flutter noise 
during playbaek when neeessary. Krohn-Hite and Hewlett-Packard 
filters with eut-off rates of 12, 24, and 48 db/oetave were used to 
minimize any high frequency noise eontamination. Signals were sent 
through a multiplexer/demultiplexer eireuit or direetly through a 
ealibration and input seleetion pateh panel for reeording on FM analog 


mee. The pateh panel could be used to seleet either data or ealihration 
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^ 
signals for Momm. and also used as a signal input monitor terminal 
during recording and playback. 

Two types of tape recorders were used for FM analog data 
recording. A Sangamo 3500, 1/2 inch, seven channel FM tape recorder 
was employed for the first four cruises, and a Honeywell 7600, seven/ 
fourteen channel, 1/2 inch/l inch FM tape recorder for the OWAX 
Cruise. Both tape recorders were adjusted for a 50 db signal-to-noise 
ratio witha ] KHz, 5 volt peak-to-peak sine wave. Allowing for 
signal intermittency and discrete amplification factors available, data 
Signals were recorded as close to + 2.5 volts peak level as possible 
without overranging, so that typical signal-to-noise ratios for the data 
Pipes was approximately 40 db. Overload lights on the input selection 
panel and Krohn-Hite filters, as wellas Tektronix storage type 
EEUoscopes or chart recorders, permitted constant qualitative 
monitoring of input data signals. A hard-wired Ubiquitous Spectrum 
Analyzer with oscilloscope and/or x-y plotter output could be used for 
On-line checks of input signal frequency content and spectral shape. The 
Spectrum analyzer was also used as a valuable tool for determining 
noise contamination of signals which would have otherwise gone 


undetected, 


ENS Data ۸۶۸" ۶ ۶ 


- 


An IBM 1130 digital computer system was employed for data 


reduction and subsequent data analysis of the analog tapes recorded on 








A. detailed discussion of the laboratory computer system at‏ ضط کک 
UCSD and various digital techniques employed for turbulence analysis‏ 


is given by Gibson (1973) 1 and Clay (1973). 


Analog-to -digital 

(A-D) conversion was accomplished with a data acquisition interface to 
mie computer permitting 15 hit samples at up to 100 KHz for single 
channel operation. Multiplexed sequential sampling of up to 8 channels 
was also possible at correspondingly slower sampling rates. Digital 
data could be stored ona 500K data storage disk or on digital tape using 
a Datum Model 100 digital tape drive unit. The computer system also 
included an IBM 1442 card reader/punch and a Calcomp 565 plotter for 
displaying results. 

A high specd stand-alone digitizer /data logger ee was also 
developed for A-D conversion shown schematically in figure 4. This 
system consists of a Tustine Electronics AD 1500 A-D converter- 
multiplexer, 2 ما‎ 011617 11 695111361 and 1 چپ‎ 
a Kennedy 8109 digital tape drive unit. The advantage of this system was 
that.it could be used for double buffered ''gapless'' sampling at up to 8 
KHz for one channel operation or at correspondingly slower rates for 
up to 16 channel operation. By ''gapless!'' sampling is meant that there 
was no time delay in sampling the data caused by writing that data on 
tape. It was then EC to concatenate several tape blocks to make 


one computer record when analyzing tapes from the data logger system. 


Computer records consisted of up to 2048 words depending onu the number 
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: High speed standalone data logger 


ro OOO C VCOI GE 7 7 


HO I'IOHLNOO 


Sac 






da vL 


e a.b VINO 
CONV 

0٦‏ ۸ھ 

1 


۷ 1 ۵ 1 0ء 
CONV‏ 

ad TASANOD 
790 

ZeIL- DOIVNV 


e 


“y “9 ET 


SE VNS 
DOTMIV 





46 


of channels sampled (typically 4) and the type of analysis programs 
msec, 

An extensive, integrated group of time series analysis programs 
was developed by Gene Dial and Jerome Dayton for analysis of the 
turbulent data and plotting results. The majority a turbulence data 
analyzed was input from digital Bares, however, Sone of the earlier 
results were obtained from data stored on disk. Figure 5 shows a 
block diagram of the time series analysis system with programs 
Brouped according to their primary functions: data collection, analysis, 
and display.. This figure is also originally from Gibson (1973). id The 
data و یی‎ programs are used for digitizing operations to transfer 
data from the A-D interface to tape or disk, aud for data transfers 
between tape and disk. Analysis programs permitted examination of 
raw digital data, statistical analysis, and spectral and cospectral 
analysis using fast fourier transform techniques. The display programs 
could be employed for a variety of functions to organize and output 
analyzed data. 

A-D conversion consisted of amplification and filtering of the 
played back analog signals to make maximum use of the dynamic range 
of the A-D converter (+ 10 volts) or data logger (+ 2.5 volts), and 
prevent aliasing of the output spectral data. Input analog signals were 


filtered at the Nyquist frequency, i.e., at half of the sampling 


frequency. The analog instrumentation described earlier was used for 
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this off-line analog signal conditioning. The digital tapes were 
initialized with an alphanumeric tape label at the beginning, consist 
of many tape records in the middle, and have an end-of-file mark at 
the end. Specification of the number of channels of input data and the 
record length was made by control communication to the digitizing 
programs DATAP and DATIN through the computer keyboard/type- 
writer. Sevcral realizations over one section of analog eo were 
normally made at different sampling frequencies and corresponding 
filter frequencies so that د‎ estimates could be obtained over a 
wide range of frcquencies. 

The primary analysis programs employed for this work included 
E LOT for time series plotting and determination of کس‎ 
factors; the statistics programs STHST and RSTAT used for calculation 
of averages, variance, standard deviation, running mean, and max - 
imum and minimum values, and CHCOR and MOMNT for calculation of 
Eprrelation coefficients and multiple channel moments or covariances; 
mie spectral programs SETUP and COSPC were used for calculation of 
powcr spectra and cospectra. The SETUP program included options 
for boxcar or triangular data windowing and digital pre-whitening/post- 
darkening. For analysis of signals with steep power spectra and high 
power at lower frequencies (particularly true for horizontal velocity 
data) the triangular windowing was found to minimize the possiblity of 


"Meaka ge" from lower to higher frcquencics.  Turbulcnt velocity and 
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Scalar spectra fall off at a rate of slightly less than 6 db/octave. 
Predicted slopes for log-log plots of turbulent velocity and scalar 
power spectrum versus frequency would be -5/3 as even in Bas, (11) 
and (12). The boxcar window employs a direct transform of the time 
series while with the triangle window the data is multiplied by factors 
increasing frorn zero to one at both ends of the data record to the. 
middle. This means that much less energy is transferred from high to 
low frequencies in the convolution operation employing the triangular 
window since the side lobes B the triangle window fall off at 12 db/ 
octave rather than 6 db/octave. Another means of reducing leakage 
was to flatten the signal spectrum from -5/3 toa + 1/3 slope by 
analog differentiation or by digital pre-whitening (i.e. transforming 
the first difference of the digital data). 

Effects of leakage on atmospheric velocity data from the MITOS 
BE cruise and subsequent techniques to minimize the leakage are 
illustrated in figures 6a, 6b, 6c, and 6d. A velocity signal from the 
hot-film and its analog derivative were simultaneously sampled at 50 
End 500 Hz to obtain several hundred 1024 point records which were 
transformed and averaged. Figure 6a shows effects of leakage on the 
boxcar window transform of the "straight" signal. The spectral 
estimates were multiplied by the corresponding frequency (f) raised to 


the 5/3 power before plotting. Sampling at the two different frequencies 








FIGURE 6 


Illustrations of techniques for minimizing 
"leakage! in spectra 
Box car window transform, leakage 
b. Derivative pre-whitening 
c. Digital pre-whitening 


d. Triangular window transform 
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accentuates the lcakage problem for the high frequency run since the 


low sampling boxcar window is wider in time and narrower in frequency. 
So there is less low frequency leakage in the overlap region of the two 
spectra. Figure 6b shows the effects of analog differkntiation of the 
signal with a boxcar window transform. Spectral estimates were 
multiplied by f raised to the 5/3 power before plotting. Due to the 
low energy of the derivative spectrum at low frequencies there is little 
energy leakage to higher frequencies and consequently good agreement 
of the high and low spectra in\the overlap region. The fall off of the 
high end of the low frequency spectrum is due to filtering. Figure 6c 
shows Fhe effects of digital pre-whitening/post-darkening (correcting 
for the cffects of the transfer function of the differencing operation) of 
` the straight velocity signal with the boxcar transform. Results are 
similar to that of the analog differentiation pre-whitening. Figure 6d 
shows the straight velocity signal transformed with the triangular 
window and the spectral estimates multiplied by f raised to the 5/3 
power. Again there is very good agreement in the overlap region. 

Due to the time consuming nature of the digital pre-whitening 
technique, the triangular window transform or analog differentiation was 
employed for minimizing the effects of leakage during spectral analysis 
of high frequency data. Spectra computed with the boxcar window at 
| Sampling frequencies below 100 Hz did not appear to be greatly 


| affected by leakage. However, any leakage effects were minimized by 
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obtaining several realizations of spectral data at different frequencies 
and matching the high frequency spectra to the low frequency spectra, 
Reduced data from the analysis programs was initially stored in 
(x,y) lists on a special disk file where it could be modified for plotting 
or output to punched cards, The plotting routines consist of log-log 
(PLGLG), semilog (PLGLN), and linear-linear (PLNLN) coordinates. 
Data could be calibrated before analysis by entering the appropriate 
computer units per physical units into the analysis program or by 
modification of the output data with the FMUL program. The FMUL 
program allows we of the x,y lists. by specification of the 


An ables A,B,C,D according to the formulae: 


B 
4 = 2 d 
(new x) A. (old x) (38) 


D 
(new y) C * (old y) * (new x) 


Recalibrated values of x and y replace the old values on the disc 
file allowing for output data of the type shown in figuresba-d. The 
special disc file may also be reentered by inputing saved data cards. 
with the card read programs RANSR, RCDPS, and RCDDT, each 


program applicable for a different data card format. 


3.4 Velocity Pitch Corrections 
Data obtained during the MITOS III and OWAX cruises was used 
W calculate the turbulent fluxes by the direct covariance technique 


00 20641 fon instrument tilt. During these two cruises optimum 
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experimental conditions were achieved, thus making this caleulation 
pos sible. The instruments were mounted at the end of the long port 
boom 18 meters from the side of FLIP whieh substantially eliminated 
the strong tilt effect experienced during Operation BOMEX, Dirft of 
horizontal velocity sensors experienced during the first three cruises 
was substantially eliminated using the hot film veloeity sensors. The 
sonie anemometer was employed for limited horizontal and vertical 
veloeity measurements during the MITOS III cruise and for more 
extensive measurements during the OWAX eruise. The piteh compon- 
N 
ent of the instrument package due to FLIP motion was measured with a 
vertical gyro system developed by R. E. Davis of oeripps Institution 
of Oeeanography and used to convert measured vertical veloeity for 
this data during digital analysis. The vertieal gyro was calibrated in 
fae laboratory prior to the FLIP cruises and mounted on the instrument 
package so that its axis was exaetly parallel to the vertical axis of the 
sonie anemometer. Calibration of the pitch gyro was possible to 0, 1% 
using a speeially eonstrueted tilt platform, also used for the sonie 
anemometer calibration. 

To eorrect for instrument olatform tilt, the measured horizontal 
and vertieal velocity eomponents were rotated in the x-z plane by an 
amount equal to the pitch angle measured by the vertical gyro. This 
Borreetion is not tive to an inertial reference frame and as sueh 


does not correet for all eomponents of FLIP motion, but only for 





56 


instrument tilt due to FLIP motion. It was also necessary to correct 
the measured vertical velocity data from the sonic ancmometer 
employed for the OWAX cruise due to a misalignment of the vertical 
sensing path heads by 1, 1° 

Denoting the measured and corrected velocity components by 
Ev. and ENDE respectively, and making a positive rotation of 


the vertical axis in the downstream direction by an angle © 


(measured), the transformations are written as 


^ ۰ 
u =u cos + W Sin 
e m = m > 











; WW cosp- u (sin OF O ai 
+ Ta sin cos Q - d (sin هم‎ + ۵ )cos y 
m m S 
ی دہ‎ ums, 
Mere y = 6 or q, and e IS the Sonic IMISAalI ۷۷٣۶٣ 7 


order O (P - 1.19? for OWAX). Since the measured gyro angles 
S 


were typically less than 2° a small angle approximation can be made 


so that 


Ses 0 => 1 (40) 


sin p ™ © (radians). 





Un 
- 





Phe measured and corrected velocities and the measured gyro angle 


re functions of time and as such consist of mean and fluctuating 





omponents. To first order the transfer rnotions can be rewritten for 





( 
the fluctuating components and the small angle approximation as 
eee 2 2 
= wW + w - - 

lw 5 ie les ‚m 9 "m P "m s 

۱ 
)41( 
, = w - u - u 
E c^ m m 0 m Ps 





4. RESULTS 


Results from the data analysis of the SOMA, MITOS I, MITOS 111 
andOWAX cruises can be identified with two types of meteorological 
conditions. The SOMA and MITOS I cruises represent warm, humid 
conditions during late summer with measurements made at night and 
late afternoon. The MITOS II and OWAX cruises represent colder, 
less humid conditions with measurements made during late winter/ 
early spriug in the early morning or mid-day. Major environmental 
conditions and important MEC. results are given m Tables 1 
2, Stability data derived from the data analysis are given in Table 3. 
It will be shown from time series analysis and spectral and cospectral 
analysis that the MITOS I cruise represents conditions similar to 
mose experienced during operation BOMEX presented by Pond et al 
(1971), Phelps and B. (1971), i and Leavitt (1973). A The MITOS 
Ill and OWAX cruises represent conditions similar to those experienced 


by the above authors also (cold, less humid conditions). 


4. 1 Time Series Analysis 


Typical signals of horizontal velocity u, vertical velocity w, 
temperature T, and absolute humidity q, representing the warm, 
humid conditions of MITOS I are shown in figures 7, Ta, Tal and ۵۰ 
Similar signals for the colder, eae humid conditions of OWAX are 


Shown in figures 8, 9, 9b, 9bl, and 9b2. Also included in the OWAX 
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MITOS I-l time series, u, w, T; q signals 
"measured simultaneously, 12 minutes of 
continuous data 
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FIGURE 7 





a. First "magnification! of MITOS 1 time 
series (x5), 2.5 minutes of data 
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FIGURE 7 


Second ''magnification'' of MITOS l time‏ 1 م2 
series (x100), 14 seconds of data‏ 
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series (x50), 28 seconds of data 
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time series is the pitch angle measured by the vertical pitch gyro on 
the instrument package. For each cruise the time series begin with 
the large scale structure of each signal and show successively smaller 
scales of the original time series. The gyro pitch is only shown in 

the largest scale time series for OWAX data. The procedure employed 
for the time series analysis can be compared to increasing the 
magnification power of a microscope to look at ES Scales of any 
particular feature of the object. 

Measurements of absolute humidity during the MITOS I and OWAX 
cruises were made with the E CE ns humidiometer so the response 
En both cases will be similar. Temperature was measured with a 
"platinum cold wire for the MITOS I (figure 7) and OWAX 4 (figure 9) 
data and with a micro-bead thermistor for the OWAX 2 time series 
(figure 8). The nature signal for the MITOS I cruise was high- 
pass filtered at 0.001 Hz (6 db/octave) so that events occurring over 
a period greater than the time constant of the filtering process (159 
Seconds) will be attenuated. For this reason the period of time 
E ered by the first figure of the MITOS I series (figure 7) is of 
smaller scale than the first figure of the OWAX 4 series (figure 9), 
but can be compared to the OWAX 2 data presented in figure 8, 

Horizontal and vertical velocity for the MITOS I data was 
measured with an x-wire anemometer, while the horizontal and 


vertical velocity measurements for the OWAX data were made using a 
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FIGURE 8 


OWAX-2 time series, u, w, Q, T, q signals 
measured simultaneously, 8.5 minutes of 


: continuous data 
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sonic NE meter. Limitations of the frequency response of the sonic 
anemometer as compared with the x-wire can be seen for the smallest 
Scales represented in OWAX figures 9bl and 9b2. It should be noted 
that the centerline of each time series does not neccessarily Correspond 
to the zero or mean level of each individual signal. 

Berge scale trends in the MITOS I and OWAX data are shown 
in figures 7, 8, 9 representing 12, 8. 5 and 54 minutes of data 
respectively. In general, ramp-like structures of 100-200 m wave- 
lengths can be seen in the temperature and humidity signals with warm 
humid air associated with upward vertical velocity and negative 
horizontal posco Changes. Both the Toe and humidity rise 
gradually to à maximum value and then decrease sharply to a minimum 
value. The gyro pitch angle is in general out of phase with the 
horizontal velocity variations and in phase with the vertical velocity 
variations, with positive gyro angles associated with larger upward 
Excursions of the vertical velocity. 

Referrring to the first ''magnification'' of the time series in 
figure 7a and 9b representing 3and 3,5 minutes of real time, the 
ramp-like structure of the temperature and humidity signals can 
Clearly be seen. The temperature signal for the MITOS 1 data is 
more like the BOMEX temperature signals presented by Pond et ae and 
Leavitt. B The original Peine process should not have any significant 


affect at this scale for the MITOS Idata. The humidity and temperature 
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FIGURE 9 





OWAX-4 time series, u, w, 0, T, q signals 
measured simultaneously, 54 minutes of 
continuous data (¢ not included in 
"magnifications'') 
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“FIGURE 9 


Second "magnification" of OWAX-4 (x200), 
16 seconds of data 
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pus Third "magnification" of OWAX-4 (x100), 
32 seconds of data 
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MAP interfaces are not only Me with maximum values but also 
with minimum values of ternperature and humidity and with what has 
been termed "'eold spikes" in the temperature signals. The cold spikes 
appear to be more pronounced in the MITOS I temperature signals than 
in the OW AX temperature signals. The OWAX temperature and 
humidity signals appear more intermittent than the MITOS I signals, 
even though the height of the sensors for the OWAX data was 350 em 
as eompared to 500 cm for the MITOS I data. Froni Counce ment 
research on BOMEX, MITOS I, MITOS III and OWAX data MeComell 
(1974), po has obtained Tui: values for dT/dx from MITOS Tand 
BOMEX which are lower by a faetor of 2 than similar kurtosis values 
ealeulated for OW AX and MITOS III data. This indieates more 
intermitteney in the temperature data for the later eruises. 

Figures 7al and 951 show the smallest seale magnifieation and 
eaeh represent approximately 12 seconds of real time. The ramp-like 
Structure of the temperature and humidity fields is still evident at this 
scale at wavelengths of less than 5m. The temperature and humidity 
appear well correlated exeept at the interfaces of the ramps where the 
largest eold spikes in the temperature oeeur. Figure 9bl shows that 
cold spikes ean oeeur during relatively quieseent periods between the 
large scale ramps. There is relatively more aetivity at higher 
frequeneies and more small seale ramps for the MITOS I data as 


eompared to OWAX data. Another feature peculiar to the temperature 
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signals for both OWAX and MITOS I data is the occurrencc of reverse 
ramps at scales less than 500 cm. These reverse ramps showing 
interfaces reversed to the mean flow direction are similar to the large 
Scale ramp structure of temperature in a heated jet as discussed by 
Gibson, Friche and McConnell (1974). is The reduction of freguency 
responsc for the sonic anemometer as compared to the x-wire can be 
Seen in figure 9bl and 9b2. 

Figurcs 7a2 and 9b2 show the existence of cold spikcs at ramp 
interfaces occasionally on the ‘upwind reversed interfaces, and also 
during more quiescent periods. The OWAX data show more inter- 
mittency than the MITOS I data and the occurrence of reverse ramps is 
in evidence, The MITOS I temperature signal also cxhibits a slower 
warming trend after an interface in which the temperature decreases 
than the OWAX data exhibits. 

The effect of intermittency of the OWAX temperature data and 
the effect of the different behavior of the OWAX and MITOS tcempera- 
ture signals at larger scales on spectral and cospectral shapes will be 
discussed in Sections 4.5 and 4.6. Differences betwcen temperature 
Spectra and sensible heat flux cospectra arc similar to those obtained 


10 ۳ 
by Pond et al, s Phelps and Pond, and Leavitt, 


E Standard Deviations of Turbulent Fluctuations and Stability 
Dependence 


The measured standard deviations for horizontal and vertical 
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vclocity, and temperature and humidity fluctuations arc given in Tablc 
BENE "ues ofthe ratios 0 /u,, € /u,, co /T, and © /q. arc 
u A W At 0 er q 2 
given in Table 3 for stability results and plotted versus -z/L in 
figure 10. Figure 10 also comparcs prescnt results for this work 
with other current results, For cach run the averaging time is given 
under "run duration! in Tablc 1. The variances and covarianccs werc 
obtained by integration under the spcctral and cospectral curvcs, 
Integration of thc u and w spectra to obtain varianccs for OWAX 
mata did not include the area under the pcak at 0. 1 Hz due to FLIP 
i 
induccd motion. An approximate corrcction was made by drawing a 
eared line through the affected region in a similar manner as was 
8 . 11 
done by Pond et al (1971) and Leavitt (1973). Any resultant error 
from this procedurc can be expected to be small sincc the peak at 
EN Hz in the w spectrum for MITOS III data was removed by an 
instantaneous pitch correction of the w signal using Eq. (39) of 
Scction 3,4. The u and w spectra from SOMA and MITOS I do not 
Show a large peak at 0.1 Hz since a large part of thc F LIP motion 
Mas removed by the cable suspension system of thcse two cruises. 
The u, values for OWAX and MITOS III were obtained from 
Mtegration of the uw cospectrum corrected for tilt effects using 


the transformation equations of Section 3.4. The uw cospectrum for 







MITOS III was corrected for pitch using Eq. (39) while the OWAX data 


Wes corrcctcd for pitch using the mcan angle approximations of 











۱ 
FIGURE 10 


Normalized standard deviations of turbulent 
fluctuations plotted versus - z/L 
a. o /u,,o /u, versus - z/L 

u 2 Ww 2 


b. oo T. ۲ Zur versus - z/L 
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Eq.(41). Any residual peak at 0.1 Hz not removed by the pitch 
| corrections was faired through with a straight line approximation. To 
Stain the entire covariance from the cospectra, integration was 
BErformed between f - 0.001 Hz and f = 10 Hz. Mhe lowest 
frequency point obtained from cospectral analysis was| between .003 
land 0.005 Hz, so that it was assumed that the cospectra went to zero 
at 0.001 Hz for the integration. The additional arca obtained using 
this integration technique amounted to an average of 7% for the uw 
covariance and 4% for the wy copatiances. The wey values torurc 
2 7 
SOMA and MITOS I runs were calculated from the relation 0381 Cp U 
۲٣۰٠٢٢ no correct pitch angle was available to obtain pitch corrected 
e 
(uw cospectra. An assumed value of C ات‎ eoo was used in 
the calculations of un . This value of CD is an average value 
| 52 
obtained from the results of Deacon and Webb (1962) who suggest 
| - = -3 
Cn IO U) x 10 (average 1,47 x 10 for SOMA and 


3 33 
MITOS I cruises), Hasse (1970)? and Miyake et al (1970) who give 


- -3 
values of Cy Sime cies O,24) x 10 2 and )۱,+ TO for 


54 
slightly unstable conditions similar to the SOMA cruises, Smith (1970) 
-3 8 
ho obtains a value of (l1l. 35 +4 0.34)x 10 , Pond et al (1971) who 
tains (1.5 + 0,26) x 107° for BOMEX results, and Mitsuta and 
7 ۱ 55 -3 PA 
Eitani (1974) who suggest 1.22 x 10 Írom Northwest Pacific 


EEcanand East China Sea data. Also included in the average is the 


present result from the MITOS III and OWAX cruises of (1. 18 € 0.20) 
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x 107? for pitch corrected data, 

The ratio Oe os. tor this work is 2.84 £0.53. This ratio 
historically has shown considerable scatter when plotted versus z/L 
and values of a fu, for neutral stability range from 1.8 obtained by 
Kaimal et al (1972)°° to 2.3 suggested by Monin » Yaglom 
۱ 1). = Their value is shown at - z/L = ome in figure 102٠ 
Leavitt روجو)‎ 7 obtained a value of 2.8 40.5 for 0 from 
BOMEX data in the stability range of 0.1«- z/L «1.0 which is 
mepresented by the dashed line in figure 10a. Mitsuta and he 
report a value of O 23.30 £ 0. 78 for -z/L.*-0.2 fron OT 
work in the Northeastern Pacific and East China Sea for measurements 
taken from a cruising ship corrected for ship motion. 

The predicted shape for the behavior of the ratios Sa fü, ; 


oo LT. and نات‎ versus - z/L is given by the relations 


d 0 rar inopi e 
0111 wW 
i «1/43 (42) 
- 0g 三 )بر‎ 2015) 
21/5 
E = CG fm L 
ی‎ e 2310 


where C , CT and C are constants and u, T., and g, are 
W q کے کے‎ > 


given by Eqs. (5) and (6). Values of the ratio AUN for the OWAX 


and MITOS III runs with pitch corrections show considerably less 
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١ 


scatter than those ratios calculated using an assumed C The pitch 


n° 
corrected values of UNAM agree well with the results of Pond et al, E 
Ey - 57 
| Leavitt and Monji (1972), are slightly higher than the results of 
5 
Wyngaard, Cote and Izumi (1971) 8 for Kansas data, and slightly 
Am "ois آ0‎ ۱ 
lower than the results of Mitsuta and Fujitani. Leavitt estimates 
E = 1,7 although his data show the trend represented by the dashed 
a dotted line in figure 10a for 0.1<-z/L< 1.0. Pond et a from 
BOMEX data obtain an average value for Oy = 1.32 + 0.09 
represented by the dark triangle in figure 10a. The dotted line of 
۱ ۱ 06 
figure 10a represents a result of Cy = 2.2 given by Monji from 
data obtained at the Bonneville Salt Flats in Utah. The Wyngaard et 0 
and Kaimal et T result of Su = 1,9 from Kansas data is represented 
۱ OD 
by the unbroken line of figure 10a. Mitsuta and Fujitani obtain a 
value of O ISSO Za for. -2/1.= 0.2. The ratios oN. and 
تا‎ Us show that the large scale velocity field is anisotropic. Spectral 


results will be given in Sections 4. 5 and 4. 7 illustrating the approach to 


isotropy for the small scale structure of the velocity field. 





The ratios - Te 2 anal - er obtained for this work are 

| Shown in figure 10b and follow E prediction of Wyngaard et ps for 
CT 20.95., Results of ی و تین‎ (dashed line with slope - 1/4) and 

Phelps and Pond (1971)? from BOMEX data shown in figure 10b are 


higher than their San Diego data and ratios from this work. The ratios 


E 2. for the MITOS I and SOMA data are lower than the results for 








20 


the OWAX and MITOS III data. In general for the SOMA and MITOS 1 
Eruises tbe values of 9 y T. are higher than the values of 7” as 
was m und by Phelps and اج‎ for BOMEX results, The OWAX and 
MITOS II ratios are more nearly equal as was found by Phelps and 
la for pre- BOMEX (San Diego) data which oM the Wyngaard 
pet E^ prediction more closely than the BOMEX data. The high values 
of oT. for the SOMA cruise are not surprising due to the uncer- 
tainty in the assumed CD value, the small sensible heat fluxes 
measured during the cruise, and the possible effects of non-station - 
Mty at low values of - z/L. Near z/L = 0,9, /T., versus z/L 


exhibits a cusp-like behavior which makes the ratio more uncertain and 


ficult to measure 


4.3 Pitch Bed Covariances andiki rot Analysis 


Movarianees determined by the eddy correlation technique are 
fairly difficult to measure as noted in Section 2.1. In addition to the 
ors caused by platform (sensor) motion, the measurements are also 
Sensitive to instrument alignment errors and coSpectral distortion at 
high frequencies due to spatial averaging and sensor separation. The 
covariances obtained for the MITOS III and OWAX cruises were 
corrected for platform motion errors due to pitch using simultaneous 
| pitch angle measurements of a vertical gyro mounted at a known angle 
with respect to the velocity sensors (sonic anemometer). Spatial 


averaging of the sonic anemometer occurs over a path length of 20 cm 








el 


so that the smallest wavelength resolved is 1.26 m . Effects of this 
spatial averaging ea the velocity spectra and cospectra occur at 

f z/U values of 10 for MITOS Ill and 2,8 for OWAX (discussed in 
bection 4, 54, Distortion of the scalar flux cospectra at high frequen- 
cies due to sensor separation is more pronounced since the sonic 
anemometer and temperature/humidity sensors were separated by 
~50 E euucing the high frequency cut-off to f z/U values of 

۲ 1.0 and 4 for the two cruises. No correction for spatial resolu- 
tion or sensor separation wa$ made since contributions to the 
covariances above 1 z/U sil are small. ا‎ 6 7 6 51113116 DE ۶٣ 
covariances is less than 5% for the OWAX cruise and ncgligible for 
Me MITOS III cruise. 

The sonic anemometer used to measure horizontal and vertical 
velocity during the OWAX cruise had a misalignment of 1.1? in the 
vertical sensing path which induced u- contamination in the measured 
ical velocity. The covariances which include vertical velocity 
Ese been corrected to account for this misalignment represented by 
pP. in the equations of Section 3, 4, 

The measured pitch angle spectra of the vertical gyro for MITOS 
MT and OWAX is shown in figure llas log Í E plotted versus 
EN vhere f is frequency in Hertz. The OWAX spectrum is an 
average from the 3 runs analyzed. Both spectra show a pronounced 


ekat {= 0.1 Hz due to F LIP induced motion. The MITOS III 








۹ 
FIGURE Ll 


MITOS III and average OWAX pitch gyro spectra 
plotted versus frequency (key to spectral plots 
Table 6) 
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spectrum also exhibits à peak at 0.01 Hz and in general has a higher 
spectral level than the OWAX spectrum for f« 0.04 Hz. Mean pitch 
angles and values of ےہ‎ for the MITOS III and OWAX cruises arc 
presented in Table 4 . Also given in Table 4 are the measured and 


corrected vertical velocities (w), the uncorrected ahd corrected uw 


covariances with confidence limits, and values of -R (f) = Y / 
uw uw 
E2 
Gu ® wy) the frequency band 0:01 <f <0. 1 calculated forech 


pitch corrected run. A comparison is also made with two other 
methods of estimating the pitch angle in lieu of actually measuring the 


Miele, Table 4 gives uncorrected and corrected scalar flux covar- 


lances with confidence limits. 


The mean pitch angle from the MITOS II cruise shows a much 
greater variation than those from all the OWAX runs, with 2d e ۱ 
For this cruise (MITOS III) an instantaneous DUI SL رر"‎ pitch 
correction was applied to the measured u and w signals by Drs. 
Carl Friehe and Frank Champagne using the first two equations of (59) 
re-written as (for تويك‎ 0): 

u) = u QV) cos « (t) + w sin Q (t) 


2 
w (t) = w cos o (t) - u) sin qo (t) 


where MIS U + u ۱ 
w(t) - W + o w' 
g(t) - ب + ب‎ ' 





TABLE 4 


Pitch Angles and Corrected Fluxes 





6 + 0 OID 0 472-023 1 012 0 J 11+ O27 
From R '* -05 oe 6° 4° 39 

uw ; 
From us W تا/‎ 4, 7? 0° ,. 7 - 1.4? 

In 1n 
W, Vertical Velocity (cm/sec) 
Measured W 54 qa 7.5 23.1 
a 

Corrected for 
Misalignment e 5 0 -1.05 -12:0 
Corrected for 
Pitch We 36 237 -8.9 Ze] 

4 
-uw (cm /sec) 
Uncorrected 352 96 160 251 
Corrected” 455 + 8 ME, 270 + 68 3243 + 7 
Measured Rw) 0, 42 0-27 0.29 0. 43 
From R ™-0.5 416 335 453 401 

uw 
o (Y/U) 447 149 209 239 

m ın 

wd (? C-cm/sec) 
Uncorrected 1. 36 ]. 26 1, 37 1. 08 
Corrected "0011 1 11176 1743220019 3917 20:2 0.08 
— 2 
wq (u gm/cm -sec) 
Uncorrected 2.79 2.39 3.48 4. 72 
Corrected 7 0.1122 2.246 2508298 9: 8652 054-5211 209 


"Measured pitch angle used 








and the bars and primcs denote mean and fluctuating components. The 
calculations were performed on the UCSD CDC 3600 computer. The 
ںہ کک‎ and covariances were then computed with the corrected u 
and w. For the OWAX data the corrected covariances were obtained 
by assuming a mean pitch angle with à small e approximation and 
using the equations of (41). This approximation for the OW AX data is 
probably valid since the measured angles and standard deviations were 
small. The uncorrected covariances given in Table 4 are those 
which would have been obtained if no pitch angle or misalignment 
Correction had been made, Corrected covariances (also given in 
Table ca RCE Collections. Phe confidence 7+ ۹:77۳ 
covariances include errors incurred from the following sources: 
Uncertainty in calculating the area under the cospectral curve (11%). 
2. Uncertainty of the pitch angle measurement, o. , FS , and the 
sonic misalignment correction (OWAX only) (8%). 
3. Possible errors incurred by fairing through the effects of FLIP 
induced motion on the cospectra (9%). 

The percentages given are average ical for the corrected uw 
variance, more sensitive to pitch effects and FLIP motion than the 
Scalar covariances. The confidence limits for the scalar covariances 
amount to an average value of 10% of the corrected covariances. 

The mean R iw value of 5 data is -0.32 + 0.08 for 0.01> f 


ES 0.1. mu! obtains a value of S = 509r his. 509p 


"TN ڪڪ‎ 
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data in the same frequency interval using a pitch angle correction 


Ealculated from the relation 


W 
A . “1 m 
Q = sin | = | (44) 





where u and U are the measured mean en and horizontal 
velocities. Values of © calculated from this relation and We 
(corrected for sonic misalignment) are shown in Table 4, as well as 
the uw values obtained using this correction, Pond et ae use an 
assumed value of سا‎ = -0,50 to estimate the pitch angle for the 
covariance corrections. Pitch angles calculated by this method are 
also Basen in Table 4 with the resultant uw au for this correction. 
The uw corrected value of MITOS III calculated using the instanta- 
ous correction procedure (for large dr ) is in fairly good agree- 
ment with the values obtained from the other two methods, even though 
mae pitch angles required for the correction are larger than the measured 
pitch angle. For the OWAX data, the values of uw obtained for the 
Borrection used by جا تک‎ are an average of 12% lower than those 
corrected for with measured mean pitch angle; the values of uw 
calculated using the R iw - -0.50 correction are higher than the 
measured pitch angle corrected uw values bysas: mueh a sesfacto eo 


for the OWAX 2 run. In his comparison of 8 m data using the سا‎ = 


| 11 
20.50 correction and the 30 m data, Leavitt ^ indicates the differences 
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in mean stresses at the two heights may be duc to errors in estimating 
the pitch angles, and that the differences would be removed by a 2,50 
rotation of estimated pitch angles at 30m , and an opposite rotation 
for the 8 m data. Although the uncertainties in the correctcd uw 
values for this work average ~ 25% the trend of the corrections which 
estimate (in lieu of measuring) the pitch angle to give higher or lower 
uw results for the OWAX data than obtained for the Suae pitch 
angle procedure, appears to support eee aa suggested method for 


removing the difference between the stresses at the two heights. 
` 


E Flux and Correlation Rcsults 
ME Cay measured fluxes of momentum, scnsible i latent heat 
are given in Table 2. The momentum fluxes for SOMA and MITOS I 
mere calculated from ut = CU g as discussed in the previous scction. 
Since the fluxes were calculated for only one height for each run it was 
not possible to test the validity of the constant flux layer assumption. A 
general trend in the heat flux results can be disccrned by noting that the 
calculated BOWEN ratio Rp = H/H, is lower for all the night runs 
than the day runs, implying radiation may affect thc amount of sensible 
heat measured by the eddy correlation method. The lowest value of the 
BOWEN ratio for the OW AX results (from OW AX 4) was obtained from 
an early morning - before sunrisc, Thc night run of OWAX 3 extended 


into the daylight hours and was of longer duration than the OWAX 4 run. 


The MITOS I-7 run, which had a BOWEN ratio of 0,25, was made 





39 


during the late REITER before sunset. Phelps and Pond IS S havc 
suggcsted that radiation may have an affcct on the larger scales of 

fluctuations in air with high concentrations of water vapor,‏ مو کک 
causing more absorption of long wave radiation and a morc homogcneous‏ 
large scale temperature ficld. This affect may be more pronounced at‏ 
night Bien ou, long-wave radiation emitted from the ocean is present.‏ 


Valucs of the overall corrclations bctwcen each pair of parameters 


x 
dcfined as r = 


are given in Table 2, The average overall 
xy 


Ze y 
correlation - for the pitch corrected results of MITOS III and 
OWAX is 0.26 € 0.08 (mean and standard deviation), in agreement with 
the AN pe Dux 203b ts T Pond et al. 2 The results 
Br Pond et A may be expected to have less scatter than the present 

200 18لا‎ Since the proccdure used to correct their data for pitch cffects 
assumed that the u and w fluctuations hada R de) = =0. 5 correlation 
mathe fr¢quency range 0.01 - 0.1 Hz. From fa ca © motor data 
the average value of Ec 0,26 0,06 0ء‎ 0 correctior لل نام‎ ۲۲067 
Leavitt dctermined the anglc required to align the sonic anemometers in 
the manner described in Section 4. 3 using measured mean vertical and 
horizontal velocites. Using this angular rotation the resultant vertical 
velocity dctermined by transformation equations similar to those of 

Eq. (39) had zero mean. 


A stability normalization can be applied to the Raw values 


obtained from the pitch corrected covariances, as suggested by the 
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results of Kaimal and Haugen (1969)1* who found that the correlation 
coefficient may vary with stability. To remove the stability dependence 
if E alized velocity spectra presented in Section 4. 5a, the spectra 
ire normalized with the expression for the dimensionless kinetic energy 
ıssipation De . If the normalized correlation SE is determined 


= — .2/3 , 2/2 
the relation r = uw ® / /(b 中 ) / (see Appendix 1), the 
uw € uu ww 


verage value of Tow = 0,25 +0.06 for the pitch corrected data using 


2/3 2٤ 9 2 
1 / SD - z/L) / (Eq. (24)). For © / 
[€ m € 
1 7 2 
verage value of -r = 0.30 £0.06, The value of -r =0.25 for 
uw uw 


: < )1 + )5 ۳ , the 


5 
lu. , agrees with the value obtained by Miyake et al (1970), (07 25) 


1 4 
md the standard deviation of the result is reduced, 


The overall values of Kae (0.48 +0.04) for the OWAX and MITOS 
I data corrected for pitch are higher than those obtained from SOMA and 
ITOS I runs (0.16 £20.11). The pitch correction to the wO and wq 
FOSpectra is small compared to that for the uw cospectra, so that the 
Bend in the r values is probably real. Phelps and Bord also 


w 9 


j:btained higher values of r for pre-BOMEX (San Diego) data than 


w 0 
hose obtained from the BOMEX data. 

Ihe average value of — from allthe cruises for this work is 

|. 43 20.11. There is a slight trend for the values of ae for OWAX 
nd MITOS III to be higher than the SOMA and MITOS I values, but is 

not ۳ as the trend in the N values, The overall correlation 
-oefficient 58 given by Phelps and De for BOMEX data (0.33) is 
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Comparable to the average value of SS for thc MITOS I and SOMA. 
Eruises (r EU Sb). 
. ۷۵ 


er values for each run are lower than the corresponding 


ug 


a values with the exception of OWAX 4 which shows the opposite 


behavior. It is difficult to notice any real trends in the 7 and E 
values and the u@ and uq covariances since they were not determined 
for the SOME cruisc. The anamolous high value of - e 20.78 fer 
run MITOS 1-1, which also has the highest wq and uq values may be 
due to drift or misalignment of the x-wire in the mean flow since the 
7 cs unaffectcd by pitch is comparablc to that obtained from 
MITOS 1-7, Thc value of uq for MITOS I-1 is comparable to values 
got uq obtained by Phelps and ر00‎ during BOMEX, 

Results obtained for n for this work agree with those obtaincd 


Ex Phclps and 0 from pre-BOMEX (To 0,8) and BOMEX 


(To , = 0,5) temperature and humidity data. Lower values of r 


9 q 

were obtained for the SOMA and MITOS I cruises (0.38 € 0.09) than 

E... obtained for the MITOS III and OWAX cruises (0.73 x 0.13). 
The correlation coefficients which include the scalars can also 


be normalized in a similar manner as was done for the r 7 correlation 
11 ] 


coefficient. The normalization for the D and zT. correlation 


6 
Coefficient is Ca ای‎ whcre y is equalto O or q, andthe 
-1/3 
normalization for the 2 corrclation is ey? / (sec Appendix 1). 
i q 


Thc normalized correlation coefficients are given in Table 5. 





SOMA 


MITOS I 


MITOS III 


OWAX 


بم 


HQ 


09 


wO 


0, 06 
05 
07109 


0. 24 
0.18 


125 
ے0 
105555 


Y 
wa 
Ora 
0. 29 
0.27 


0. 36 
2202 


0.29 
(000 
0.36 


0,14 
0. 14 
(- 59 


li 


ieee tt > = Normalized Correlation Coefficients 


0082 
02995 


0. 12 
0. 14 
0.29 


UW 


Se 
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-1/3 1/6 


The fi Irs D and $c were not included in the normalizations 


since the differences between the individual runs amount to less than 


2 s 1/6 


4% if P. were included. The trends noted for the 


measured correlation coefficients are still applicable, however, less 


pronounced. The average r 0 correlation coefficient for all the data 
wW 


59 


BS 0.21 0.1, comparable to the result of Miyake et al (1970b)' * of 


r 0.24 x 0.2 for near neutral conditions (-z/1 20.04). The 


w 8 


trend in the £ ga correction coeficients tor OWAX ٥۹ MITOS IJI 


E52 +4 0.08), and the SOMA and MITOS I cruises (0.25 £ 0. 04) 
-1/3 


would not be removed by including the factor ©. in the normalization, 


a 


4. 5 opectra 


In this section horizontal and vertical velocity (u,w), tempera- 
ture (0), and specific humidity (q) spectra computed from the four 
cruises will be presented and discussed. The method of spectral 
presentation will be to plot initially the individual run spectra in 
absolute units for each cruise as log fo (f) versus log(f), where f 
is frequency in Hertz. The individual runs are then averaged for each 


)/N 
XX. 


EEuise and plotted in normalized form as log NON 


= > 080 iate normalizi uantit 
( ود‎ ) , where D is the appropri nor izing quantity 


ENS. u., T, , etc.) versus log f, ( = fz, /U) , the normalized 


N 


frequency, following the trandition of atmospheric spectra presentation. 


‘Table 6 gives the cruises, runs, appropriate plot symbol and 








TABLE © 


scr caicema Cospectral Plot Information 


Dimensional Sampling 


Freguencies 


Absolute Spectra ۳۰۱ voy mbpo | 





SOMA 1 ۱ 十 16, ' 130, بر‎ 
2 3 16, 150, 208 

3 z 16, 130, 200 

MITOS I IL \ + 2 0000 1۳ 
07 > 2.2 On ٦٦ 

MITOS III 2* O 2.09 20 01 
OWAX 2 | + ئ ئھ۳0۶"‎ Ow 
3 > 2000. ۶ 7۲0۵ 

4 X 2.24 20,7 5 


Average Normalized Spectra 


SOMA F 
MITOS I X 
ETTOS III 4 


OWAX حا‎ 7 


*Figure 17 shows the MITOS III w spectrum corrected 


Bor pitch (<q) and uncorrected (+) sampled at 3,0 Hz. 


104 





105 


sampling frequencies used to digitize the analog data. The record 
length for each spectral computation was 512 words /record/channel. 
Hach spectral realization at different sampling frequencies was 
calibrated in absolute units, 20th decade averaged, and matched to the 
Hower frequency spectra. Not all the points from each spectral 
realization are shown in order " reduce overlapping points and because 
Bae end points of each realization are affected by m filter cut-off 


(Nyguist frequency). 


۷ 


tocity Spectra 

‚ Absolute horizontal velocity spectra for SOMA, MITOS I and 
MITOS III/OWAX are shown in figures 12, 13, and 14. The vertical 
velocity r. from these cruises are shown in figure 15 - SOMA, 
figure 16 - MITOS I, figure 17 - MITOS III and figure 18 - OWAX. Only 
mie w-spectrum of MITOS III is shown corrected and uncorrected since 
Bn instantaneous pitch correction was employed for this cruise. The 
averaged u and w spectra are shown in figures 19 and 20 normalized 
with u. and plotted versus fN . 

۱۱۱۳ large peak in the u er w spectra between 0.05< f 

€ 0.2 Hz due in large part to wave induced FLIP motion is more 
pronounced in the OWAX spectra (z = 350 cm) than in the spectra 
from SOMA, MITOS I and MITOS IM. For both the OWAX and MITOS 


Mi cruises the instrument package was rigidly attached to FLIP, while 


for the SOMA and MITOS I cruises the instrument package was 








FIGURE 12 


SOMA horizontal velocity spectra 
plotted versus frequency 
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FIGURE 15 





MITOS 1 horizontal velocity spectra, 
plotted versus frequency 
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FIGURE 14 





MITOS II and OWAX horizontal velocity spectra 
plotted versus frequency 
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FIGURE 15 


SOMA vertical velocity spectra 
plotted versus frequency 
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FIGURE 16 


MITOS I vertical velocity spectra 
plotted versus frequency 
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EPICURE 17 





MITOS III vertical velocity spectra plotted versus | 
frequency corrected for pitch (qd) and uncorrected 
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BIGURE Ys 


OWAX vertical velocity spectra 
plotted versus frequency 
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suspended on a cable. This suspension technique did appear somewhat 
successful in reducing the wave induced FLIP motion effects on the 
velocity spectra. Figure 17 shows pitch corrected and uncorrected w 
spectra from MITOS III to illustrate that the instantaneous pitch 
rection docs appear to remove the peak at 0.1 Hz for this data, 
although only one high 0 spectral value at that frequency is 
shown. 

Ihe peak at 0.1 Hz in the velocity spectra may also be 
due in part to the effects of waves on the turbulent flow field as 
suggested by the results of Lai and Shemdin (1971). ae These investi- 
gators E that simple mechanically generated water waves had a 
noticeable effect on the velocity spectra (u,w) and cospectra (uw) 

Bt the wave frequency. A maximum percentage increase in e and 
B. of 24% and 9% respectively was found at heights above the water 
between 10 and 50cm. The water waves did not appear to have any 
affect on the spectra in the inertial subrange. 

The sonic anemometer horizontal and vertical velocity for 
the OWAX cruise, and vertical velocity spectra for the MITOS III 
Cruise have not been corrected kan: spatial averaging and exhibit a fall- 
off in ۵ level E i er ne nee ۱ and ot een 
velocity spectra. With a path length of 20 cm the response of the 
Sonic anemometer should be accurate for f, 5 10 for MITOS III and 


N 


SN = 2.8 for OWAX velocity data according to the relation 
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FIGURE 19 
Average horizontal velocity spectra normali: 

2 

with u, plotted versus normalized frequency 
(f z /U) | 
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FIGURE 20 





Average vertical velocity spectra normalized 


2 
with u, plotted versus normalized 1۲ sS 





e 
e 


Û 


| 


0 


p 


: 1.6.4.540 





fz /U = z/2nt | (45) 


Encre z - 1250 cm (MITOS III) and 350 cm (OWAX), and 4 = 20 cm. 
The averaged u and w spectra of figures 19 and 20 

normalized with u are in fair agreement with the Kansas 

data of Kaimal et al (1972) ed by the dark line in both 

meures. In general the spectral values for the SOMA and MITOS I 

€ gher than the OWAX and MITOS III spectra for d 2, possibly 

due to uncertainty in the u, "values and also due to sonic response 

for the OWAX spectra. 


^ 


To remove the z/L dependence of the velocity spectra, an 


2 





additional normalizing factor is required, namely n , where 
o = = € from Eq.(23). The inertial subrange form of the u 
E 
EDectra given by Eg. (11) may be rewritten as 
er 
EN uu E) 5 u $ DS 2/3 
2 u 2/3 € N عد‎ 
ua (211K) : 


For this work an & value of 0.53 was used based upon direct 
u 


dissipation results from MITOS I, MITOS III and OWAX 4 in the 


50 
concurrent work of McConnell (19774). For an 0 ， - 0.53 and., nte t p is RSS 
E - 0.4, (46) becomes 
f 中 IE.) 
N uu N 2/5 
__ x 9 47 
2— 273 0.29 f. (47) 


u, 中 
“€ 
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FIGURE 21 


Average horizontal velocity spectra normalized 


/3 


2 2 
with u, and > (assuming total production 
1 
= dissipation of kinetic energy) plotted versus 


normalized frequency 
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FIGURE 22 





Average vertical velocity spectra normalized with 


2 2:25 
, and e / plotted versus normalized frequei 
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FIGURE 23 





Average horizontal velocity spectra normalized 


2 
with u“ and o 2 (Wyngaard and Cote a 


و 
expression for normalized dissipation plotted‏ 


versus normalized frequency 
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FIGURE 24 
Average vertical velocity spectra normalized w: 
2 2/3 ۱ ۱ 
u, and $ plotted versus normalized frequency 
2 
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Seesimilar relationship can be obtaincd for the w spectra by assuming 
isotropic conditions and making use of the relationship 中 = 4/3 09 
WN uu 


in the inertial subrangc: 


IN ww EN) - 0.391 -2/3 
رت‎ et. Ä (48) 
Ue 
















Figures 21 and 22 show the ts u and w spectra normalized with the 
20) 5 -1/4 21 

expression po ZEE) / TL] / 
1 \ 

dissipation of mechanical energy is balanced by buoyant and mechanical 


from Eq. (24), assuming 


production. This 2007017861172 8.۲16017 ıs compared to the results of ا‎ 
(dark line figures 21 and 22) who obtains from BOMEX data values of 
ENZS £0.05 and 0,37 £0.07 for the constants in Eqs. (47) and (48). 
Figures 23 and 24 show the fu and کت‎ spectra normalized with 


Ze 


رت 
the expression Pr (1 + 0. 5|z/L| / ) from Eq. (27) based on the‏ 


Besults of E aard and Cote 2n who used values of کے‎ 00 ٤ 
K = 0.35 to obtain their normalized velocity curves with values ot the 
Eonstants for Eqs. (47) and (48) of 0.296 and 0.394, 

Constants obtained from this work for Eq. (47) are 
0.29 40.06 for $^ 4 and 0.23 + 0.06 for ?e, and 0.40 + 0.09 
End 0.32 + 0 9 E Eq.(48). A more sensitive test of the two 
expressions for E. will be discussed in Section 4.6. A test of the 


M bl r 2 لع مدي‎ ae 3 ] ns ] 中 = 4 3 中 which 
Validity of the inertial subrange relations hip 0 / en 


Ssumes local isotropy is shown in figure 25. The average ratio 
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Spectra ratio @ /® plotted versus 
ww uu 


normalized frequency 
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& /中 for each eruise is plotted versus f . The OWAX and 
ww uu N 

MITOS IIl spectra have approximately the required ratio while the 
MITOS I data is higher and the 5OMA. data is lower than the required 
ratio, The high and low ratio values of the MITOS I and SOMA data 
may in part be due to x-wire drift experieneed during the two eruises 
and discussed in Seetion 2. 2c, due to salt contamination on the probes. 
Both MITOS Ix-wire runs experienced 30% drift in the velocity 


calibrations but were matehed with eup anemometer results whieh 


agreed well with the mean of the before and after run ealibrations. 


E 


. The SOMA x-wires were ealibrated only before eaeh run so that the 


| absolute drift is unknown and only run 1 was matched with a cup 


—— wu 


anemometer. For this reason the averaged normalized SOMA velocity 
Speetra are probably inordinately high and are not ineluded in the 


estimates of the eonstants in Eq. (47). 


4.5b Temperature and hurnidity speetra 


Absolute temperature spectra for SOMA, MITOS I, and 


MITOS III/OWAX are shown in figures 26, 27 and 28; absolute humidity 


| speetra for the cruises are shown in figures 29, 30, and an Miero- 


bead تیور توق‎ S were us e io measure er کے‎ e doe CIR TOT ne 


۰ 
d > 
ےھ‎ 


dic SOMA « cruise E 6 - OWAX runs 2 Sl 3. The epe EC 
of MITOS I, MITOS III and OWAX 4 were obtained from cold wire 
temperature data. Humidity was mcasured with the Lyman- Alpha 


humidiometer for all eruises. 
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FIGURE 26 
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SOMA temperature spectra plotted vorsus frequen 
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FIGURE 27 


MITOS Itemperature spectra plotted versus frequ 
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FIGURE 28 


MITOS IH and OWAX temperature spectra 
plotted versus frequency 
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FIGURE 29 


SOMA humidity spectra plotted versus frequency 
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۲ 
MITOS I humidity spectra plotted versus frequency 
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MITOS IlI and OWAX humidity spectra 
plotted versus frequency 
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Ihe temperature spectra of MITOS I appear to exhibit a 
high frequency peak at 0.8 Hz for run 1 and at f> 1,5 Dr forming 
i. Temperature spectra from BOMEX presented by ہے ور‎ 1 and 

10 

Phelps and Pond also exhibit a high frequency peak for 0.25 <f< 
0.9 Hz for 8 meter data. The temperature spectra from MITOS IJI and 
OWAX 4 do not exhibit similar peaks and also do not follow the 

-2/3 . 
predicted 1 mertial subrange behavior for 1<f< 100 JIz. itis 
difficult to ascertain the behavior of the thermistor spectra for f» 1 Hz 
Since above that frequency thé spectra roll-off due to thermistor 
response. 

The humidity spectra from all cruises are, similar and appear 

-2/3 Me ! 
to follow a 1 behavior for f 20.8 Hz and do not exhibit any high 
frequency peak. The MITOS III humidity spectra obtained at 1250 cm 
exhibit the most pronounced inertial subrange for f 0.2 Hz. 

Differences between temperature and humidity spectra and 
between individual runs for the temperature spectra are shown in 
figures 32 and 33 which present average temperature and humidity 

Z 2 
spectra normalized with T, and q, plottcd versus z/L. The 
7 -2/3 ۱ 
humidity spectra collapse quite well and exhibit a en behavior for 
IN 7 0.3. The temperature spectra do not collapse and only the MITOS 
Itemperature data exhibita -2/3 behavior for iN 2. The tempera- 


ture spectra of MITOS III and OWAX appear to follow a -]/4 slope 
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PIGURE 3727 


Average temperature spectra normalized with 


2 
p. plotted versus normalized frequency 
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FIGURE 33 


Average humidity spectra normalized with q, 


plotted versus normalized frequency 
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mer 2 1/2 decades for EN OSI rression for the temperature 


nd humidity spectra in the inertial subıange given by Eq. (12) can be 


ewritten as 


o (f 
16 (f) ۱ b o 71/3, -2/3 7 
Y. eee) ok ee N 


merc y= T or رو‎ Le is the dimensionless scalar production given 
IB Eq. (19), D tnessealar subrange constant, and K= 0 ne 
dimensionless scalar production of temperature and humidity variance 
has been assumed equal and in balance with the dimensionless scalar 
dissipation given by Eqs. (30) and (33) for temperature and humidity. 
fis assumption may not be valid as will be discussed in Scction 4. 7. 
Figures 34 and 35 show the avcraged temperature and 
humidity spectra normalized as in Eq. (49) and plotted versus IN 
Also shown in these figures are the temperature and humidity spectra 
10 TE 10 
of Phelps and Pond and Leavitt for BOMEX data, Phelps and Pond 
56 ۱ 56 
for pre- BOMEX (San Diego) data, and Wyngaard and Kaimal for 
Kansas temperature data, normalized in a similar manner. Data from 


| i | ۱ ue -1/3 _ 
this work has been normalized with 9 = (1-16 z/L) e - 


و _ 1/4- 


and o 
MES d 


sell 
[(1-16 z/L) 25] as were used by Leavitt for his normal- 


ized spectra, This normalization is not significantly different from 


l 5 -1/2 
that employcd by Wyngaard 2 who used p - 0.74 (1-9 z/L) , 


1/3 2/3 -1/2 
62 ) 


5( 
中 and K= 0.35. The icine and 
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250. 5| z/L]l 








FIGURE 34 


Average temperature spectra normalized with 


2 
T,, ® 


al 
اما‎ 


(dimensionless temperature gradient) 
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H 
and pi plotted versus normalized 
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EPICURE 35 


Average humidity spectra normalized with 


2 

q.s ® (dimensionless humidity gradient) 
-1/3 . "a 

e plotted versus normalized 


frequency 
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58 
Wyngaard Kansas data, represented by the solid line in figure 34, 
if normalized with the expressions in this work and K = 0.4 would be 
approximately 3% lower (based on - 4/L = 0.2) than shown. The 
temperature and humidity spectral data for BOMEX and pre- BOMEX 
cruises shown in figure 34 and 35 was obtained by faiting a mean line 
کت‎ | | | 
Enrough Phelps and Pond's ^ spectral results. Since their results 
2 2 Ber ۱ 
were normalized with So and ES , the normalization of'ficure 34 
2 2 T 
was calculated based on average 7 /T, -4.23 for pre- BOMEX and 
6.63 for BOMEX (this includes correction for Kin their definition 
一 一 2 二 7 
BOT. - w0/Ku,), E. /q,, of 4.84 for pre- BOMEX and 5.25 for 


ale 
Ld bd 


BOMEX, and average - z/L values of 0.165 (pre- BOMEX) and 0. 195 


à 2 | -1/3 ١ 
BOMEX) to calculate the values of ل‎ b and RT . The 
-1/3 ۱ 
expressions for $i and E. / used to normalize the Phelps and 


1ا 027 

ena’? data were the same as those used by Leavitt. The Phelps 
10 ۱ ee eae $: 

and Pond and Leavitt BOMEX temperature and humidity spectra 


۳۰۰۰ in figures 34 and 35 agree for f,_>0.05 and arê higher than 


N 
San Diego data of Phelps and ang" and spectra from this work. The 
normalized humidity spectra of Phelps and Senado differ by 10% due 
ito the different conversion factors. | This difference is greater by as 
Much as 25% for n > 0.2, probably caused by low pass filtering 
Of the San Diego data at a lower frequency than the BOMEX data as 


noted by Phelps and Pond. P The 10% difference in thc normalized 


10 | Mr 
average humidity spectra of Phelps and Pond may be due to scatter 





161 


Bethe Phelps and Ponce. spectra but appears to be real when compared 
Eo the humidity speetra of this work, 

The general trend in the normalized temperature spectra is 
for the warmer more humid BOMEX and MITOS I datalto be higher than 
the Kansas land (dry) spectra and the less humid San Dicgo spectra. 
Ehe temperature spectra from MITOS III and OWAX 4 do not follow the 

-2/3 


predicted iN behavior for fN >10 as shown by the Kansas data. 


The MITOS I temperature spectra have the same shape as the BOMEX 
. ` -2 [ss 5 
data and follow the predicted EN behavior for x >10. The San 
10 
Diego thermistor data of Phelps and Pond does not extend beyond 
iN = 8 and has the same shape as the thermistor spectra of the SOMA 
10 
and OWAX 2 and 3 runs. The Phelps and Pond data do not extend 
MS ا‎ 
beyond EN - 10 and Leavitt's | average temperature speetra shows 
scatter and response limitations above IN = 10, so that any comparison 
of spectra for fy >10 is tenuous, 

The normalized temperature speetra of figure 34 are higher 
in spectral level than the corresponding normalized humidity spectra of 
۱۱۶۱۲6 35 for hn > 1. No attempt was made to average all the tempera- 
ture and humidity spectra for comparison of the two on one plot due to l 
the large spread of the individual spectra from ۱۹ c rusos 


Comparison of the constant Bo and DE which can be obtained from 


a. (49) will be made in Section 4. 7a. 
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4.6 Cospectra 

Results of cospectral analysis are presented in this section. 
Cospectra from the four cruises are presented ina manner similar to 
the presentation of power spectra, with averages and plot symbols 
remaining the same unless otherwise noted. Two realizations Cor each 
cospectral computation were made at the two lowest sampling [requen- 
cies shown in Table 6 for each run. Since the correlations between the 


۰ 


variables being compared in the cospectral analysis are less than 
Emity, cospoectral cstimatcs En have larger uncertainties than spectral 
Estimates at the same frequency, particularly for low frequency 
estimates, Uncertainties in the cospectral levels at lower frequencies 
have been minimized by the long averaging periods. The number of 
estimates of thc lowest frcqucncy cospectral points at f - 0.004 Tlz is 
between 170 for the OWAX 3 run and 32 for the MITOS III data. 

The cospectra from the MITOS III and OWAX cruises have been 
corrected for internal misalignment of the sonic anemometer (OWAX 
only), instrument pitch, and wave induced FLIP motion using the 
procedures discussed in Section 4. 3. The uncorrected wy cospectra 
m ae MITOS 4 e SOMA اع‎ Ti also 00 ad: به‎ scalar 
lux cospectra are less sensitive to motion effects than the uw cospoctra. 
The momentum flux cospectra from MITOS land SOMA arc not 


' a Y qe o 1 
mecscntcd since no pitch angle was avajlable to obtain pitch correct: ri 
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uw cospectra, and the (uw) cospectral estimates are sensitive to 


drift of the x-wires experienced during the two cruises. 













4. 6a Cospectra of momentum (uw) 

Absolute (uw) cospectra from MITOS III and OWAX plotted 
versus frequency (f) are shown in figure 36. In the frequency range 
0, 01 > 1 > 0.1 Hz used to determine the value of the spectral 
correlation coefficient R w) , the spectra exhibita eee fairly 
constant peak. No cospectral estimates were made below 0.001 Hz. 
Cospectra of momentum flux presented by COEM i and extending to 
0.0001 Hz exhibited considerable scatter and negative values in the 
frequency ٣۰۰۰۹۱٦ 0۰0001 > ۲ > 0.001 Hz , so ۳۳۳ ۱ oniy 
small (< 5%) contributions to the uw covariances Eos 0,001 7 ۰ 

The averaged uw cospectrum of OW AX and the MITOS III 
cospectrum normalized with a and plotted versus log Eu are shown 
lin figure 37. Wyngaard and Cote 24 predict an inertial subrange 


behavior for the uw cospectra given by the expression 


TCE CE) 
Ben... رریں‎ )50( 
uw N 


u 
where the function G 12715) " 1.0 determined empirically from Kansas 
IT Kaimal et 2150 for 2 s =. > 0 i with a predicted m of 

he constant يان‎ 0.048 . The cospectral curve predicted by Wyngaard 


nd Cote 00 iS Shown by the dark line in figure 57 where the divergence 








FIGURE 36 


MITOS III 3nd OWAX momentum flux (uw) 
cospectra (corrected for pitch) plotted 
versus frequency 
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FIGURE 37 


Average uw cospectra normalized with 
2 
u, plotted versus normalized frequency 


ale 
in 


(comparison with other results) 
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of the curvcs for ar « 0.4 represents diffcrent low frequency shapes 
or unstable conditions )-2 5 z/L«0). The results of Pond et S 
Eom IA EX are represented by the envelopc dcfined by the dashed 


8 
lines o vc 37, Pond's data had to be converted by a factor of 


E 4 2 
E fu, - 1.74 since the data were presented 88 ل‎ with C ۰ 
i NV 


The cospectra of MITOS III and OW AX follow the predicted 


41/3 
SN E behavior for hn 20.2, but are lower in cospectral level 


than the Kaimal et ig 6076 E cospectra, The value ۴ 
E. 7 0.01 2 0.002 for this data also illustrates the lower cospectral 


6 ۷7۵15, Lower cospcctral levels for fn 20,2 forthe MITOS Il and 


BIVAX data cannot be completely explaincd by spatial averaging crrors. 


-4/3 


The two cospectra collapsc on onc another and exhibit a EN 


behavior ovcr nearly a decade in the range 0.2 SUE <2, below 









۲ = 2.8, where the effects of spatial averaging occur for the OWAX 
ENS. The pitch corrections applied to data of Pond et ur using an 
assumed R iw value arc much larger than equivalent corrections 
applied to cospectra of this work, but the corrections are only 
nificant at low frequencies (f «0.1 Hz). 
j formalized ww cospectra of Leavitt for 30 m data also ل‎ 
| | a | MM than T 8 T data, | 


4/3 


Owevcr, the fall-off in the 30 m data is greater than 9۹ ۱ and was 


exhibit lower cospectral levels abovc ES 


ttributed to noise in one channel of the 30 m sonic anemomcter, 


Ithough no cxccssive noisc was obscrved inthe u or w signals of 





۱: ۱۳ ۱ ۳ 38 


Average uw cospectra normalized with 
2 
u, plotted linearly versus normalized 


aly 
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frequency (comparison with other results) 
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the sonic anemometers used during the OWAX and MITOS III eruises, 
noise in a single channel may have eontributed to the faster roll-off in 
the Ea The behavior of the uw eospeetra at higher frequeneies 
En 7 0.2) may be due to effeets other than spatial averaging and high 
frequency noise. Both normalized eospectra from the MITOS III and 
OWAX compare with one another and have the same shape as the uw 
eospeetra obtained by Kaimal et Ia The spectra appear to be shifted 
to lower frequeneies, suggesting that produetion of kinetie energy may 
occur at larger seales over the oeean than over land under eertain 
eonditions. 

A eorreetion was applied to the OWAX cospeetra for the 


internal misalignment of the vertieal sensing path, however, this 


distortion may have eontributed to lower eospectral levels for EN O 


The normalized uw eospeetra are shown plotted linearly 
P. 11 
versus logí,. in figure 38 and eompared with Leavitt's averaged uw 
N 
sa 
eospectrum (cireles). Leavitt obtained a value of M UOT 
0.02 from the 8 m normalized eospeetra. The eospeetra from this 


| 11 
Pork are lower than Leavitt's data los aN pound show evidence ol 


possible = مار‎ noise contamination with nose ph u 7 


eospeetral values for the MITOS II 23 above £ = Ol The ee 


N 


frequeney eospeetral levels (f,,<0.01) are, in general, slightly higher 


N 
11 
than results of Leavitt, ان رم‎ ue to the diíterenees in piteh 


Borrection proeedure. 
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4,6b Cospecira of heat and water vapor flux (wO , wq) 
wO Cospectra 

Dimensional w@ cospectra are shown in figure 39 - SOMA, 
40 - MITOS 1, and 41 - MITOS Ul and OWAX, The SOMA wé6 cospecira 
do not exhibit similar shapes; possible reasons for this behavior have 
becn proposed earlier and include the low wö covariance values 
indicating near neutral conditions and suggesting effccis of non- 
stationarity, uncertainties in the w and O calibrations, and lack of 
EU ch corrections. The SOMA data also have the shortest averaging 
time and only extend to a low frequency cut-off of 0.015 Hz so that 
uncertainties in the calculation of cospeciral levels are equivalent to 
the other 000 data. 

Te MITOS I wO- cospectra have a peak at f — 0.5 Hz 
and exhibit a shape similar to the shape of ihe BOMEX wO cospectra 
shown in figures 42 and 43. The dimensional w8 cospectra of MITOS 
Band OWAX corrected for pitch are different than the MITOS I wO 
Bospectra and exhibit a broad peak in the frequency interval 0.01 «f « 
EL 5 Hz ۱ similar to the peak observed > the MITOS MI and OWAX uw 
cospectra. The cospectra are not greatly affected by FLIP we 
motion, with only one averaged cospectral estimate from ihe OWAX 2 
iun at £7 0,1 Hz significantly lower than other estimates (corrected 
Bor in the integration). The cospectral points circled at f ** 0.002 2 


۱ -3 ۱ ۱ 
werc used to extcnd the integration to 10 Hz as discussed in 
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FIGURE 59 


SOMA heat flux (w0) cospectra plotted 
versus frequency | 
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FIGURE 40 





MITOS I w8 cospectra plotted versus frequency 
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FIGURE 41 


MITOS II and OWAX wO cospectra plotted 
versus frequency corrected for pitch 
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Section 4. 2. For subsequent cospectral plots, those assumed points will 
miso be circled, 
Figures 42 and 43 show averaged wO cospectra from all 
thc cruises normalized with u,0, the negative hcat flws (= -w0) 
plotted logarithmically and linearly vs log AN : Differences between the 
MITOS I and MITOS HI and OWAX data include an apparent shift of the 
MITOS I data to higher fN values than the MITOS III and OWAX 
cospectra and a narrower peak in the MITOS I cospectra and lower 
cospectral levels for fy <0.05. The averaged SOMA data also exhibit 
a peak at Ey + 0.03, probably not real since the cospectra were not 
Eorrected for pitch effects which are most pronounced for low frequen- 
cies. The peak in the MITOS 1 data is probably real and is comparable 
to the BOMEX results of Pond et E shown in figure 42 (envelope 
shown by hatched lines) and the results of eee averaged wO 
normalized cospectra shown in figure 43 (circles). An estimated 
average of Pond's et ٦ pre-BOMEX (San Diego) results for wO are 
also shown in figure 42 and agree wen with the wO cospectra of MITOS 


EI and OWAX 


re. 
بر‎ 5 


۱ 61 . 
` “" The prediction of Wyngaard and Cote for inertial sub- 
range wO cospectra bebavior, using arguments similar to those used 
Er obtain Eq, (50)for the uw cospectra, is wrilten as 
. E 
f... d Ey ۱ ۵ 


UN wO 
ل‎ = (Y MEZUN 
uno. wO w N 


4 t "esae vf ای‎ F 








FIGURE 42 


Average wO cospectra normalized with 
u,T, plotted versus normalized frequency 
eS 


(comparison with other results) 
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FIGURE 43 


Average wO cospectra normalized with 
u,T, plotted linearly versus normalized 


frequency (comparison with other results) 
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where H_ (z/L) = 1.0 for -2s z/L=s0, and the constant «& ~= 0. 14 
b w 6 
E 56 ۱ 
iom results of Kaimal et al. Ihe wO cospectra curve suggested by 
۱ 00 CENE. 
Kaimal etal is shown (dark line) in figure 42, with the divergence in 
mae predicted curve for en <2 representing the variation of low 
frequency shapes for different unstable conditions; thellower line is for 
mi; = O0 and the upper line is for z/L - -2. 
T. M -4/3 
1۳ cccospect/a of MITOS Dexbibit à IN behavior for 
ME> 2, with w _=0,29#+ 0.04, larger than the g of Kaimal et 
N wO wO 
56 11 
۳ by a factor of 2, but comparable to Leavitt's ^ result of e = 
ESO + 0.08; The low frequency cospectral points of the MITOS I data 
11 
E 0.012 > en < 0.03 are lower than those of Leavitt since no pitch 
morrection was applied to the MITOS I data. This difference in the 
cospectra may also be due to uncertainties in the low frequency 
| 1 
"cospectral estimates of spectra from this work and Leavitt's ; cospe a 
- The MITOS III and OWAX pitch corrected w9 cospectra 
56 f 
meree well with predicted curve of Kaimal etal for Kansas data and 
8 
e San Diego results of Pond et al, except at higher frequencies, where 
there is evidence of spatial averaging effects and possible high frequency 
noise effects above E zl کت‎ for the MITOS ۱۳ data, and above 
IN = 0.3 for the OWAX data. Contributions to the wO covariance are 


not greatly affected by the falloff, illustratcd by plottiug the normalized 


cospectra linearly vs log EN in figure 43. 





The MITOS II value of A O = 0.06 +0.02 also indicates 
that the high frequency cospectral levels are low when compared to the 
results of Kaimal et al. = The w0 cospectral levels are higher than 
the uw cospectra for fy ~0.2, and lower than the ‚uw cospectral 
levels for Bn «0.2. This is particularly true for tha MITOS 1 w0 
cospectra when compared to the predictcd uw cospectra of Kaimal 
Et E and the MITOS III and OWAX cospectra. Assuming ہے‎ 0 
the vertical transport of heat is 6 times more efficient than momentum 
for MITOS I and only slightly more efficient for the MITOS a data for 
smaller eddies En > lje The difference in cospectral shapes between 
the MITOS 1 w0 جو‎ and the predicted and calculated uw 
spectra implies that the vertical transport of heat for MITOS I 
conditions is much less efficient than the transport of momentum for 
larger eddies En <1.0). The difference in the MITOS III and OWAX 
wO cospectra and uw cospectra at low frequencies is notas 
pronounced as the MITOS I wO cospectra. 

wq Cospectra 
Cospectra of water vapor flux (wq) in absolute units for 
the four cruises are shown in figures 44-SOMA, 45-MITOS I; and 46 - 
MITOS m id OWAX. The = Mons 5 کر‎ id 
Similar shapes and exhibit a broad peak over the frequency range 
Ol < £< 0.3 Hz, the same interval as the peak in the MITOS Iii 


mang OWAX we cospectra, and broadcr than the peak of the uw 








FIGURE 44 





SOMA water vapor flux (wq) cospectra 
plotted versus frequency 
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FIGURE 45 
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MITOS I wq cospectra plotted versus frequency | 
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FIGURE 46 


MITOS Il and OWAX wq cospectra plotted 
versus frequency (corrected for pitch) 
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cospectra. There is evidence of wave induced FLIP motion in the wq 
Eospectra, with low cospectral points occurring at f^ 0,1 Hz 
ected for prior to integration). 

The expression for the inertial subrange form of the 
normalized wq cospectrum derived using the same procedures as 


vere done for the uw and wO0 cospectra is written as 


-A © (|) 
ايم‎ - GC OT EDE ۱ )52( 
wq w N 


S ماه‎ 
id hd a a 
S 


where it is assumed the ھ+"‎ BN (NG) = H(z/L) ^ ] for unstable 
conditions, and A is the cospectral:subrange constant. This 
expression was also used by AE for presentation of normalized 
EN cospectral results. NE wq cospectra from all the cruises 
are shown in figures 47 and 48 plotted as log f Em versus 0 


N 


linearly versus log Er ; the negative sign in the normalized cospectra 
results itemm the definition of Ax» SO that A .E wq , the negative 
water vapor flux, The normalized wq cospectra collapse with the 
exception of the high frequency end of the OWAX wq cospectrum. 
Reasons for this behavior of the OWAX cospectrum have been proposed 
in the E us Pe errem. “tt should = EE that the high fecu edo ast l 8 3 
of the MITOS III and MITOS I cospectra are similar, indicating that 


Enly the OWAX cospectra are affected by high frequency noise or sensor 


separation errors. The subrange constant determined from the data 












FIGURE 47 


Average wq cospectra normalized with usd. 
plotted versus normalized frequency 7 
(comparison with other results) 
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FIGURE 48 


Average wq cospectra normalized with u,q, 


plotted linearly versus normalized frequency 


(comparison with other results) 
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(excluding OWAX)is 0 270706 0.02. This valu’ oro is 
wq wq 
higher than en - 0.01 obtained from the MITOS IIl and OWAX uw 
= 56 TNI 
cospectra, and Se = 0.048 suggested by Kaimal et al. Leavitt 
obtains a value of m Ull £ OU, higher 0 ۸۸۰۰۶۰۰1٠٣ 
different from the 7" for this work. eo wq cospectra 
also exhibit a sharp fall-off at EN =“ 1 for 8 m data so that the values 
of 6 may be slightly underestimated, but are much lower than the 
0 values for MITOS I (0.29) and BOMEX (0. 30). 
The wq cospectra oi Fond et Ae for BOMEX and San Diego 
ل‎ : 1 
data are shown in figure 47; Leavitt's | average wq cospectrum from 
BOMEX results is shown in figure 48 (circles). No attempt was made 
8 
me average the BOMEX results presented by Pond et al, and the 
cospectra are represented by the envelope of the two unbroken curves 
in figure 47. The pre- BOMEX results are approximated by the dashed 
curve of figure 47. All the normalized wq cospectra agree, with the 
Exception of the OWAX cospectra for AN ۸0ے ول‎ 01۳7016 020 0 ۰ 
be somewhat underestimated at high frequencies En > 2) because of 
noise effects and instrument response. 
oo. Following the results of Leavitt, the cospectral constants 
0 and 0 are related to the spectral constants مق‎ and f by the 
wO wq 0 4 


equation 


2 (53) 
[o wo 8 


Bo 
Pa 
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Meat assumed that Hiz/i) = M(z/L) and P 2 . Using the values 


from the MITOS I results, a _/a 266 cipio +65 
wo 9 q 
z 11 ۱ 0 
Leavitt's  resultfor e dicc - 2.7 implics that Bo = 7,3 = . The 


high frequency fall-offs of the w8 and wq cospcctra make this 
relationship betwcen Bo and de somewhat cl A morc 
Bensitive test of thc rclationship between thc scalar spectral constants 
will be madc in Section 4, 7. 

The w8 and wq normalized cospectra arc not signifi- 
EI diiferent with the cxception of the MITOS I cospectra. The 
MITOS I data suggest that under certain conditions heat is transported 
تک‎ tically more efficiently than watcr E 58777311627 scales, and 
less cfficiently than water vapor at larger scalcs, The crossover point 
at En = 0,1 corresponds to a scale size of 42 meters using an 
Bveragc height of 44 mcters for the MITOS I Srs 

00101 un values indicatc that the vertical transport of 
water vapor is relatively morc efficient than thc transport of momentum 
for smaller scalcs (y7 0.2), although the difference is not as 


pronounced as the vertical transport of heat and momentum. 


f the Diga RR behavior (0.2 < N < 4.0) ofthe 


9 ۰ ۰ 
که‎ ۰۰ 
M 


OWAX — E wq om is zen. it tote ties ine t: ہے‎ 
of temperature and humidity variance occurs at larger scalcs for colder, 


less humid conditions, than for the warm humid conditions of MITOS J 


and BOMEX, 
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4.6c Cospectra of horizontal heat and water vapor flux 

lug, ug) 

The dimensional u8 cospectra from the MITOS I and 
MITOS III and OWAX cruises are shown in figures 49 and 50. 
Corresponding dimensional uq cospectra are shown in figures 51 and 
827۰ The MITOS III and OWAX 18 cospectra and the uq cospectra 
from all runs exhibit a broad peak in the frequency interval 0. 01 < f < 
0.1 Hz and have similar shapes. The MITOS I uß cospectra have a 
different shape than the other cospectra and exhibit à peak in the 
frequency interval 0.1<f< 0.4 Hz for the MITOS l-I run, and 
between | ۱ 20.7 Hz for the MITOS I-7 run, The uo and uq 
cospectra have not been extended to 0.001 Hz since the low frequency 
behavior of the cospectra is not established over the ocean, and few 
measurements of horizontal heat and water vapor flux cospectra have 
been made. The u8 and uq cospectra presented by Pond and 
کک‎ and the Kansas results of Kaimal et aes for u8 are 
referred to in this section for comparison of results. 

Average uO and uq cospectra from the three cruises 


-have been normalized with Uy, 6 and u, پک‎ and سو‎ as ird " MS ef ) 


د 
zy 2 ,‏ 
n‏ ۰ ۰ 
7 € ۶و a = 5 2 ٩‏ .^ ۰ ۰ 
0 8 ۰ یں ٠ a‏ 
De * . o, 4 3 a m E^ i ; 4 e‏ ۴ 
9 8 ۰ ۰ 


versus log T in a figures 53 and 54, 7 وہب‎ linearly versus T: 5 


in figures 55 and 56. The expression for the inertial subrange 


behavior of the normalized u8 cospectrum following the prediction of 








FIGURE 49 


MITOS 1 horizontal heat flux (u8) cospectra 
plotted versus frequency 
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FIGURE 0 


MITOS lli and OWAX u0 cospectra 
plotted versus frequency 
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FIGURE 51 


MIT OS I horizontal water vapor flux (uq) 
cospectra plotted versus frequency 
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FIGURE 52 


MITOS III and OWAX uq cospectra 
plotted versus frequency 
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Wyngaard and Cote is given ie 578۵ 


f (f) 


N ug N 3/2 


M MU (54) 


IN 


ne‏ ماه 
aD‏ `„ 


where H (z/L) = 1 for -2 5 2/1 <s0, and oe is 5 7 
subrange constant. Kaimal et بر"‎ suggest a value of e = 0. 035 : 

The -3/2 cospectral slope is not as well defined for the uO 

cospectra as for the w8 cospectra, Wyngaard and Cote a predict a 
slope of -3 for the inertial subrange of ihe u9 cospectra from simi- 
larity arguments. A slope of -3/2 appears to fit the u8O cospectral 
behavior at high frequencies of the three Cruises and is 5 

calculate the a 6 values for comparison with Kaimal et 2278 results. 
The MITOS TI a cospectrum exhibits higher cospectral levels 
for دآ‎ 20.2 than the MITOS III and OWAX cospectra, and lower 
cospectral levels than the other two cospectra for AN > 7 
cospectra have different subrange levels, also shown by the subrange 
constants for each cruise: مر‎ = 0,08 (MITOS D, 6 = 0.04 (MITOS UI), 
and e = 0,012, with uncertainities of - 20% in the estimates mM 

|. differences may in part be due to spatial averaging effects and high E : 
En however, Te کر یں‎ to factors of 2 
extend over nearly half 2 decade between 0.3 ہت‎ AA O 


normalized u9 cospectra of Phelps and Donde is represented by the 


envclope of the hatched area in figure 53. The results of Kaimal 





FIGURE 53 


Average uO cospectra normalized with u,' T, 
plotted versus normalized frequency 
(comparison with other results) 
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FIGURE 54 


Average uq cospectra normalized with u,q 
plotted versus normalized frequency 
(comparison with other results) 
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et ° are shown by the solid curve. The Phelps and Eo 
Eospectra are comparable to the MITOS I cospectrum for A > 0.04 
and are eher than the MITOS I data for fy > 0.6. The indicated 
value of the cospectral constant for the Phelps and ond” data is 

E 0.1 fora slope of -3/2. The cospectra of Kaimal et Ns do not 
significantly differ from the MITOS 111 cospectra. The trend in the 

09 values suggests relatively more efficient horizontal heat flux 

E sport for warm humid conditions. 


Assuming a similar behavior for the normalized uq 


cospectra, the expression for the inertial subrange is given by 


En) ۱ 


r. 3 


7 


3/2 


oe (55) 


.1 
where Ea is the cospectral constant and it is assumed that the 

function M y ES M, = ] for unstable conditions (-2s$ z/L=0). The 

ug cospectra also appear to follow a -3/2 slope for fA 0.3 so the 

slope value is retained for Eq.(55). The Phelps and once data is | 

Eo: esented by the e. of the 7 ds of TA. 54, The uq 
cospectra also exhibit differences in the interval 0 « IN Rb 5 

Beiniilar to the diiferencés im thé uà cospectra. Tho subrangG constants ^ 7 
determined for each cruise arc 8 =0.022 (MIRO Tand 8 = 0.008 


(MITOS III and OWAX) with uncertainties of 2096 and 3095 


| | 10 
respectively. The constant determined from the Phelps and Pond 








FIEURE 5 


Average uO cospoctra normalized with u,T 


plotted linearly versus normalized frequency 
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FIGURE 56 


Average uq cospectra normalized with u,q, 
جه د ع”‎ 


plotted linearly versus normalized frequency 
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۳۳۱ 15 € - 0.0353. Thetrend in the a values suggest a 
uq uq 
relatively more efficient transport of horizontal water vapor flux at 
small scales for warm, humid conditions. 
The normalized uq cospectra have lower cospectral levels 
Baan the uO cospectra for EN 0.1, most pronounced for the MITOS I 
cruise which also has higher uq cospectyal levels than the us 
cospectra for IN > ( corresponding toa scale of 50 meters for 
z = 5 meters. This suggests that horizontal transport of heat is 
relatively more efficient than, the horizontal transport of water vapor 
Wi small scales. Conversely, for the MIT OS I cruise the opposite is 
true for large scales. 
The rclative efficiency of the horizontal and vertical 
transport of heat and water vapor can be estimated from the ratios 
9 la al /@  . The a ./a ratios for the three cruises 
ug wê uq wq ud wo 
Me 0.27 (MITOS I), 0.36 (MITOS IU) and 0,27 (OWAX); the 
B /a moos are 0.27 (MITOS ), 0.1 (MITOS UD) and 0.27 
uq wq 
(OWAX). These results suggest that the vertical transport of heat 
and water vapor is always more efficient than the horizontal transport 
56 ۱ 
for small END ۰ 2) scales. , Kaimal et al. suggest this behavior for. ... 
horizontal and vertical heat transport for conditions over land; the _ 


erse is also true for larger scales En .4. 9 





219 


4. 6d BT ira parepperature end humidity (Sq) 

lhe dimensional cospectra of temperature and humidity 
are سس‎ ٠۰۰16 ظ0231) 7۰۱5 ۱۹ ب07۰‎ ۸ 58 - MITOST, and 59 - MITOS 
manga OWAX. The 0q cospectra of SOMA, MITOS II and OWAX 
0000300 ا‎ 011131 shapes. The Oq cospectra of the MITOS I cruise are 
pot so typical and exhibit a 007 peak at the frequency interval 
۲ ظ‎ > ! > 1.0 Hz , similar to the peak in the MITOS 1 wO cospectra. 

The @q cospectra have been normalized with T,q, ; 
consistent with the normalization of the cospectra presented in Sections 
4,6a, b, and c. The expression for the normalized 8q cospectra may 


4 


be written as 


5/3 ارات 


11 OG N - 
— -go H D 56 
med. 6q كيو‎ En 0 


ae 


as suggested from the previous cospectral results. The z/L variation 

of He T is not known and is assumed equal to 1 for unstable conditions 

as was done for the previous cospectral stability functions, The choice 
Bethe -5/3 slope for the 8q cospectral subrange is based on 

Similarity arguments and dimensipnal analysis by Wyngaard (personal 
munication to c. Friehe and a. pico 1972), “The ae a 


6q cospectra are shown plotted as log fo 


q OT versus log EN 


9 


in figure 60, and linearly versus log EN in figure 61. The normalized 


q cospectra exhibit differences in the 6 -q relationship at high and 
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FIGURE ST 


SOMA temperature-humidity (0q) cospectra 
plotted versus frequency 
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FIGURE 58 


MITOS I Oq cospectra plotted versus frequency 
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MITOS III and OWAX 6q cospectra 


plotted 
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versus frequency 
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FIGURE 60 


Average 9و8‎ cospectra normalized with T q 


‘plotted versus normalized frequency 
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FIGURE 61 





Average 04 cospectra normalized with ا‎ 


plotted linearly versus normalized frequency 
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low frequencies between the MITOS I cruise and the other cruises, 
similar to thc differences observed in the normalized w6 and ué 
cospectra. Since the 0q cospectra arc presented in this manner, 
unccrtainties in the assumed and calculated u, values will also affect 
the cospectral levels, but will not change the shape of the cospectra. , 
The SOMA 6q cospectra may Ms be affeeted by tempcrature calibra- 
tion crrors and uncertainties in the wO estimates used to calcolate T 
due to x-wire (w) drift and possible non-stationarity effects discussed 
in Section 4.2. n 

The highest normalized cospectral levels are exhibited by 
the SOMA 6q eospectra reflecting the lower values of the we 
eovariance. The shape of the SOMA 6q cospectrum is similar to the 
shape of the MITOS III and OWAX eospectra, indicating that the differ- 
ences in the SOMA w08 cospectra shapes is probably due to not correct- 
Eu or piteh effects. The '"crossover' point exhibitcd by the wO and 
E cospectra of MITOS Iis also well defined for the ٤ cospectra of 
MITOS Lat f, - 0.15 («30-50 mcters scale size) n 
| The MITOS III and OWAX cospectra appear to cxhibit a 
A O a a o 
Me bectively, AAA a Scale Size of 230 m gai the OWAX | 
Seta and 540m for the MITOS II data. A slight pcak in the 8m 6 


11 
Borrelation presented by Leavitt is also apparent at fN ^ 0.01 


corresponding to a mean seale size of 530 m. 








MITOS I 
MITOS III 
OWAX 
E 
"Leavitt . 


Kaimal et al 
and 


Wyngaard 


and Cote 61 


TABLE 7 - Cospectral Subrange Constants 


Q 
uw 


9 


0.048 


w8 


0.29 


0.06 


(ES 


wq 


0.08 


ud uq 


0. 05 0.025 


0.04 02008 


0.008 0,008 


020935 - 


9q 
0. 19 
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The trend in the normalized wO0 , uO , wq, and ug 
normalized cospectra for the colder, less humid runs to show small 
contributions to the "tcc at higher frequencies is also exhibited 
by the normalized 0q cospectra. Differences in cospectral level occur 
for ET Iu eue crossover' point. A slope of approximately -5/3 
is exhibited by all the normalized 9q cospectra in the interval 
0.4< En < 1.5. The subrange constants for the last three cruises were 
۳۳ 0,19; MITOS IL, & “0.12; and OWAX, « = 

Oq 8۹ md 
0.09. The uncertainty in the subrange constants determined from 
figure 60 was 10-20%, 

Although no information concerning the phase relationships 
between temperature and humidity were obtained for this work, it should 
be noted that for their BOMEX 0q data, اص‎ and Phelps and 
ena?” show only small negative phases at lower frequencies than the 
BOMEX Oq phases. The negative phases indicate that the temperature 
leads the humidity. 

A summary of the cospectral subrange constants is given 
in Table و‎ | 


5/3 E | 
E? k Spectra and Comparison of Flux Techniques 


In this section comparisons are made be cen lhe different flux 
Estimation techniques discussed in Section 2. The assumptions invoked 
to simplify the budget equations of kinetic energy and scalar variances 


are also examined. The directly calculated fluxes of momentum, 








1 


sensible and latent heat used as a basis for the comparison have been 


presented in Tables 2 and 4, and discussed in Sections 4, 3 and 4. 4. 


4. tà Comparison of eddy correlation and inertial dissipation 
techniques 
Momentum F lux 

Two basic assumptions invoked for estimation of 
momentum flux by the inertial dissipation technique are that the budget 
equation for kinetic energy may be simplified to the forms given by 
Eqs. (24) or (25), and that fire eed enel horizontal velocity 
spectra is represented by a -5/3 power law behavior as shown in 
Eq. (11). Eqs. (24) and (25) suggest two forms for the dimensionless 


viscous dissipation giben by 


E = (1-16 ا‎ L (57) 


1 


E 0 ۶ 


2 


(58) 


The expression for D . is derived from the assumption 
1 


that total production (mechanical + buoyant) is equal to dissipation of 


۱ ۱ ۱ E " 27 35 2 : l - SC T um ب‎ 
. mechanical energy. Wyngaard and Cote. have obtained the expression ,-... 


For D. from Kansas data as discussed in Section 2.6. 
2 


The velocity power spectral relation of (11) can be rewritten 
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5/3 o 
ee 
و کہ کو رر رو رر رتا‎ 2 
NERONE: HE (59) 
u,, 27 


where k is the radian wavenumber determined from the relation 


k = 2nf£/U by Taylor's hypothesis, K = 0.4, 9. given by Eq. (57) or 


ل ن 


2 
(58), and u, = -uw. 


In figures 62 and 63, average velocity spectra from all 
Bene cruises are normalized as given by Eq. (58) divided by @ , and 
u 


plotted logarithmically versus log k, the radian wavenumber in units 


- 1 00 
cm . The u, values for the SOMA and MITOS I spectra have been 


calculated based on an assumed CD 2-173 ome , and those for the 


MITOS III and OWAX directly measured and corrected for pitch as 


discussed in Section 4. 3. In figure 62, e from Eq.(57)is used in the 
1 
normalization, and in figure 63, P from Eq.(58). The value of the 
E Z 
Bubrange constant used in the normalization is œ =0.53 from direct 
u 


dissipation measurements to be presented in the next section. Also 


shown in figures 62 and 63 is the hot-film velocity spectrum from OWAX 


4 for comparison with the average sonic anemometer spectra. A 


EN... m^ عي‎ 


cR 


Beonstant value of. 1 for the spectra indicate a slope of -5/3 and ^ «e uas 


me ement between the direct and inertial dissipation methods, It 

Should be noted that a value of 1 (one) is not an exact agreement between 
2 

the two methods since the directly calculated u, values were used to 


ERlculate the values of Pe : 








FIGURE 62 


2959 | 
k average horizontal velocity spectra 
! 2 -2/3 
normalized with & u, $c (Kz) plotted 
u * 
۰ 1 l 
versus radian wavenumber (ordinate represents 
agreement between direct covariance and ner se 
E ۱ 2 
dissipation estimates of u,) 
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FIGURE 63 
5/3 
k average horizontal velocity spectra 
2 -2/3 
normalized with o. DM © (Kz) / plotted 
n 


versus radian wavenumber (ordinate represents 
agreement between direct covariance and inertial 


"T ھا‎ ۱ 
dissipation estimates of u,) 
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Tuous figure 62 the constant spectral levels and the 
Eorrespondingemnterval for the spectra are SOMA (1.7) and MITOS III 
19. 93) Ds 0.02 «0.8 SUM O AOS tor 500002 >7 5 
AE mond MITOS I (I. TO: O. 1] <k < 1 NL . From figure 63 the 
constant spectral levels over the same intervals are SOMA (1,3), 
MITOS 1 (0.94), MITOS II (0.73), and OWAX (0.58). The MITOS 
I velocity spectrum cxhibits the narrowest inertial subrange; all the 
spectra have a gradual approach to an inertial subrange behayior for 
1 > ۵ RT , and exhibit effects of wave induced FLIP motion in the 
interval 0.001 <k<0.003 cm” * 

i The high values of n — levels for the SOMA 
and MITOS I data as compared to the OWAX and MITOS III spectra, 
suggest that the assumed value of CD - 1.3 × ۳ may be low. 
Another source of error in the normalized spectral levels of SOMA 
and MITOS I spectra is x-wire drift, As discussed in Section 4, 5a 
this drift in the MITOS I data was corrected for by matching x-wire 
calibrations to cup anemometer mcasurements and using the mean of 


before and after run calibrations. It was not possible to completely 


eorrect all the SOMA data so the spectral levels arc probably 


۶ ٩ g ۰ 
.. ۰ ۰ .. ۰ ^ ۰ .. ۰ ۵ 8 wie", ? ۰ 
k - wu. ] A مو ا‎ 3 . ^ = .. 4 E O 2 «>»? a : te x e ۰ ۰ ۰ 7 ۰ 
a. o I~ . t ee ےو ہہک‎ S : s cibos ore we » ¿Y y Te à É ۰ I ^ ۰ s ^. 


inordinately high. 
Matching the MITOS I x-wire to cup anemometer results 
may have also introduced an error contributing to high spectral levels. 


From analysis of the OWAX sonic and cup ancmometer mean velocity 
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measurements, the mean velocities agree within 4% except for the 
OWAX 3 run which indicated 5-7% cup over-speeding. Analysis of 
the MITOS II hot-film and cup data indicated up to 15% cup over- 
speeding. This over-speeding for MITOS III may be high since the 
hot-film was mounted horizontally with respect to the flow field instead 
of vertically as was done for the OWAX 4 hot-film. Wyngaard and 
22 ۱ 5 
Cote used a 10% correction for cup over-speeding for their 
Kansas velocity data based on analysis of the Kansas cup anemometer 
: 62 63 
measuremenuts by Izumi aud Barad (1970). Hyson (1972) and 
64 
mondo et al (1971), suggest that cup over-speeding is 1% or less, 
64. ۳ ۱ 
however, Kondo et al indicate that upto 4-7% cup over-speeding 
Ever land may occur in the daytime. If a value of CD A 
and a 10% cup over-speeding correction is applied to the MITOS I 
data, the subrange spectral levels would be 0.8 (figure 62) and 0,61 
Ze 
(figure 63). This correction does not include the effect on D due 
to changes in the z/L values but the differences are not significant 
E570) . 
Using the corrected values for the MITOS I Spectra and the 
Tos H I OBRAK result, the aver aus weis al Bevel for وت‎ 
is 0. 85. a 0. 07 add a use "63 is 0. PRO, 004 . The CR 
limits only reflect differences in the average corrected spectral levels 
and as such are probably low. Better confidence limits for the spectral 


^ 


levels would be ~ 30% referring to the confidence limits on u, in 
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Table 4, 
The trend in the values indieates that for this work the 
inertial dissipation method underestimates the value of Te = Ike 


underestimation is more pronounced using the ®_ expression 
€ 
Tn Me 
proposed by Wyngaard and Cote. Pond et al obtain good agreement 


à C 
between direetly measured u, values (corrected for piteh by the 


۱۲۳ (f) = -0. 5" method) and the inertial dissipation method using an 
uw P = 


a= 0255 and the expression 


2 


uy, = [K(cB) z] e i 


m 


derived from Eq. (22) assuming no flux divergeneec, where B= g 


—— — 


E 274 2 0 dT x 1077 wa) (equivalent to the second term ın Ea. (22). 
Pond et al, however, do not include the cffeets of stability on the mean 
veloeity profile (P 0 in their ealeulations. Ineluding the AM 
dependenee in their ealeulations would increase the u. values obtained 


by the inertial dissipation method by 25%. 


For near neutral conditions (-z/L <0. 1) Miyake et al 







(1970)? 


apply the inertial dissipation technique with the assumption that 
| production and dissipation are in equilibrium (Eq. (3)) and obtain higher a 
Fvalues of u, than those obtained by the eddy correlation technique. 
36 
From measurements over land Hieks and Dyer (1972) obtain a value 
of 0 =0,54 using the relationship given by Eq. (56) with the u, values 
u 7 


Obtained by the eddy correlation technique and confining the estimates 
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P 


of œ to conditions of u, 720 cm/sec and H,- 10 eeu 
a 5 


S 
The subrange spcctral levels of figurcs 62 and 63 are also 
an catión of the ratio of dissipation of mechanical cnergy to 
production (including stability effects for production). Using the 
inertial dissipation techniquc to obtain the dimensionless dissipation 
and Ir. -z/L) to obtain the dimensionless shear production, m m 
finds a 7% differcncc (production greatcr) Den the two for 
B 0.53 (a value of^0.96 on figure 61). 0 : also finds the 
dimensionlcss turbulcnt transport to bc approximatcly half of that found 
by Wyngaard and Cote. if this additional correction were applied to 
nues results on figure 62 it would correspond to à value of 
approximatcly 0.80 (average of 0.96 and 0.64) in agreement with 
the valuc obtained for the results from figure 62. It should be notcd 
Ehoraf a simalar correction werc applicd to data from this work the 
results would be lower than those obtained from figurc 62. 


Hcat and watcr vapor flux 
The assumptions required for estimation of hcat and water 
E. flux he india dis وت‎ E dc Similar to'thosc 
employcd for estimation of momentum flux when applied to the nz ھت‎ 
Manco ed ons ana کی ا 77ھ080‎ E. MEG eur 
assumptions that the cddy transfer cocfficicnts for hcat and momentum, 


and watcr vapor and momentum are equal and the dimcnsionless scalar 


gradients for temperature (2) and humidity (2) are equal. 
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Ihe expressions for the temperature and humidity spectra 


from Eq. (12) can be rewritten as 


RS P (k) وت‎ 
= - g (61) 
2 us 9 
2 E o 
and 
S 中 E 
= 62) 
2 -1/3 1 ( 
Ode Be e 


where the starred quantities are calculated by Eq. (6) using directly 
| ۱ 3ؤ"‎ 
measured covariances, e, = P z (1-16 z/L) , and e as by 


Eq.(54). The difference between 8 , and e will not be 
1 2 
Significant for comparison of the temperature and humidity results. 
Estimates of Ba varo cons laerably from sales tot 
8 pot 65 
0.40 (Pond etal ), Wyngaard and Cote » bPaegqunausasbond (19071 
11 
Bo values of 1.0 or larger suggested by the results of Leavitt of 
14 47 ۱ 
B. E OS Stegen etal (1973), and Gibson etal obtain Bo ct 
oe E . ۰ 3 7 66 
using the direct dissipation technique. Boston and Burling (1972) | and 
Boston (1970) suggest Bo = 0.8 for direct dissipation measurements 
yer 2 mud flat,” Values of 8. “from direct dissipation measurements 
or this work are presented in the next section and vary from 0.7 to 


ES. There is a dearth of values of TE in the literature from open 


cean measurements and include cstimates of Ba = 0,4 by Pacquin 
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and mend” and Pacquin (1972)98 using second and third order structure 
functions, and 7 = 0.25 suggested by .تہ‎ from flux profile 
measurement results of Paulson et al which imply 6 - © += © , 
E m a 
the Ba, and 2 values given above are based on the definitions of scalar 
wariance dissipation from Eqs. (30) and (33), 
Temperature and humidity spectra from the cruises 
normalized in the manner of Eqs. (61) and (62) are shown in figures 64 
and 65, plotted versus log OE . A constant value of the spectra 
indicates a spectral slope of -5/3. Due to the variation of Bo and 
E estimates as discussed above, the normalized cospectra have not 
been divided by the subrange constants as was done for the momentum 
flux comparison. The value of the spectral levels thus represent the 
implied Bo or 3 required for agreement ee 1116 71 
dissipation estimate of nd or a , and the value determined from 
direct covariance measurements. Since directly measured covariances 
are used to calculate Py and 9 ۲٦6 1371017694 ۲16011765 NOL exact For 
Biven values of Bo and e ‚ the spectral velg divided by those, 
values represent the ratio of scalar dissipation by the inertial 
dis sipation technique to the scalar production assuming p -do = 


۱۳۳ ۰. 2 e. D y E ۱ T را‎ vit RE Bing od a hee | 
| | eye E | a 
m-16 z/L) / : 

The normalized temperature spectra of figure 64 exhibit 


considerable differences. The thermistor temperature spectra from 


the SOMA cruise and OWAX 2 and 3 runs exhibit a roll-off in spectral 











wt 


FIGURE 64 


人 人 average temperature spectra normalized 
و‎ -1/3 -2/3 
with 2 Tr, کل‎ / (Kz) / plotted versus 


radian wavenumber (ordinate represents Bo 
required for agreement between direct covariance 


: : کرت‎ : e 
and inertial dissipation estimates of T.) 
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FIGURE 65 
در‎ uf 
k average humidity spectra normalized with 
2 -1/3 -2/3 
Zee DR (Kz) plotted versus radian 


wavenumber (ordinate represents e required 
for agreement between direct covariance and 


f وی‎ ate 5 2 
inertial dissipation estimates of ۱ 
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level for k> 0.03 omy due to instrument response. The cold wire 
spectra are still increasing, however, and do not exhibit constant 
spectral levels below k = 0.8 mr With these<ception- or tne MIT OS I 
data. The MITOS I spectrum exhibits a constant spectral level of 

2,1 0.2 (implied B.) over the range of wavenumbers 0.06 <k < 

1.0 me for a slope of -5/3, The MITOS III spectrum has a 
Borrower subrange in the interval 0.5 «k «1.0 qu with an implied 
Bo - 2.5 2 0.2 . An inertial subrange for the OW AX 4 spectrum is 

E -bonmonexistent. The OW AX 4 spectral level for 0,7 <k <1,0 m 
۱۳۰ ما‎ The slopes exhibited by the MITOS III and OW AX 4 
aeui for k <0, 9 e are -4/3 and -5/4 POS See Vel Sled 
of the inertial subrange slope of -5/3. If the assumed value of 

Cp 29710 3 X 1077 were increased to 1.5 x w for the MITOS I data 
as discussed for the momentum flux comparison, the implied Bo 
would be increased to a value of 2.5 , comparable to the MITOS III 
Besult. The implied Bo values of enna and Pond et a for BOMEX 
data. (4.2 +1), Phelps and Dances for San Diego data (-1,0), and 


e DAS. l | 
Wyngaard and Bee! and Kaimal et al (0.4 £0.05) are shown at 








10 
the left-hand side of figure 64, The Phelps and Pond results from 
me هم‎ ۰ . s v. 9 E Ta ide 5 e iei. : م‎ 2 moe ہت کر‎ : 4 97 I تفع‎ e. A is Ke E tre ne; ۰ 
San Diego are based on spectra obtained from measurements of 
temperature using a thermocouple with frequency response out to 


~3 Hz (3db point), and may not sufficiently reflect the temperature 


5 ) 
2 JO 
behavior at smaller scales. The Kansas result of Kairnal et al and 





250 


Wyngaard and Cote = was obtained from measurements of the 
production and turbulence transport of temperature variance and 
assuming the imbalance between production and transport equaled 
dissipation using 20 cm fine platinum wire thermometers corrected 
for spatial resolution effects, 

The spectra shown in figure 64 indicate that large errors 
(overestimation) may result when estimating the sensible heat flux by 
the inertial dissipation oe A value of B م‎ = 620° 15 a 
for estimates of wd comparable to the direct covariances using this 
technique over the ocean, in conjunction with temperature measure- 
ments out to scales as small as 6 cm in wavelength. The results also 
imply that dissipation of temperature variance may exceed production 
by a factor of 2 or more (for Bo ze Direct dissipation 
E -surements presented in the next section also appear to support 
this implication. 

The normalized humidity spectra from the four cruises are 
ES. in کت کہ‎ The 2۳01 CT O 0.870 aver 

mime interval 0.01 «k «0. 04 cie ۰ The MITOS Ul spectrum has a 
Ea value of implied ' B, = 0. 29 T 0.03 E a considerably wider ir etes 
range (0.002 «k «0.1 MAD than the other spectra, The SOMA and 
OWAX values are the samo with an implied a 2n 0005 Lhe 


lowest spectral levels are exhibited by the MITOS I spectra in the same 





interval as the OWAX and SOMA spectra with the implied 
-3 
2 = 0.13 € 0. 02 (0. 15 for e -1.5x 10 ). The lower spectral 
-] 
values for k 0.04 cm are due to instrument response. The high 
spectral levels of the MITOS III data for k < 0. 4 A may be due to 
high frequency noise in the humidity signal at lower frequencies than 
the filter cut-off (Nyquist frequency). The average implied EN for 
this work is 0.21 +0.05, which is lower than the value obtained by 
65 OS e 
Pacquin and Pond, and Pacquin of 0.41. The value of 8 = 0.25 
A q 
11 ۱ 
suggested by Leavitt is comparable with the value of 5 from 
figure 65, and suggests e. = 中 if production and dissipation of 
n 

humidity variance are equal. This implies water vapor is transported 
more nearly like momentum than heat. 

Referring to the results of Section 4,6 and the expression 
Er oc . 0 given by Eq. (53), the ratio of the cospectral constants 

wO wq 


for heat and water vapor flux implied B. Ze 77# . Using average 


values of Bo = 88 ٤٦ B. سے‎ ۱7۶٣٣٣۴ Bo lod Ba 


not significantly different than the result predicted by the cospectral 


Bubrange constants, 


x ۰ e. ‘ "T ۰ E ۰ E J^. c ^ Pa ۰ ۳ ۰ É ۰ نه‎ 
eo 5 * y ۱ 002۷0000 EN ۴ 
۾‎ 


4. 1b Comparison of eddv correlation and direct dissipation 


6ھ 


- 


Estimation of the fluxes of momentum and sensible heat by 


the direct dissipation technique have been made from measurements of 
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the mean square time derivatives of veloeity and temperature. The 
mean viseous dissipation (€) and mean dissipation of temperature 
varianee C ; ealeulated using Eqs. (13) and (14) from the MITOS I, 
MITOS III and OWAX 4 data are given in Table 8. No eorreetions for 
69 +۶ 
Hekstad's (1965) modifieation of Taylor's hypothesis have been 
: 737107 ہے‎ ۱ 
applied as was done by Wyngaard and Cote , sinee the correetions 
are typieally < 2%. 
The € and X o estimates for MITOS I were obtained by 
51 
Steve MeConnell by matching high frequeney speetra of velocity and 
temperature time derivatives to the eonverted veloeity and 


* 


È 2 
o gt) xf /(2n) ) shown in figures 14 and 28. 


The X6 values for MITOS II-2 and OWARX 4 from een. are 


4 


temperature speetra (® 


ealculated for shorter time periods than the direet eovarianees but 
are representative values for the runs. The MITOS III € value is an 


average from three estimates at the beginning, middle and end of the 


51 
run, and the OWAX 4 (MeConnell) value is from the beginning of the 


run. Allthe € and X values have eonfidence limits of approxi- 
mately 20% or better. 


The @ and By values shown in Table 8 were 
E Sti t u i ab. ld : ۰ 


EE." ٭‎  ہ‎ ۰ ۰ 4? 
7 er تو‎ 


p^ p 5 BE A m m ا‎ ae fetu m rtm a Pis i mos ie e stero, Ee dis z x N 5 ا‎ 3 
aleulated using the direetly measured € and Xo estimates and 


: 
i Y TUUS sU. OmU e Cee ف‎ iis و‎ o, "LE! 
مو"‎ ee pe ۱11-۰ 77 


employing the relations 


E LEES 
$.. (k)=0 € k (63) 


uu u 





TABLE 9 - Direct Dissipation Results 


€ [m a 


u 
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Bo 


ani Ni (SC see) 


MITOS I 1 re 8.6 x Mee 0.53 
7 108 8.6 x IONS 0. 50 
TA 
MITOS III 2 7 6.3 x 10 0.53 
3 


OWAX 4 433 2.6 x 10. 0. 54 


2.4 


i29 


0.68 


0.67 





€ k (64) 


where u and 8 are spatial derivatives by Taylor's hypothesis. 
Velocity and temperature derivative spectra from MITOS I-1 and 7, 


and OWAX 4 are shown in figure 66, 67 and 68 (from McConnell). = 


-2/ 00 


= 3 = 
The spectra are normalized with € (velocity) and Xo € 


(temperature) and multiplied by N so that the spectral levels 
plotted linearly represent or and Bo , anda constant vae indicates 
a + 1/3 slope. The 0 is given as log (nk) where y is the 
EN rov microscale orae mound 10 (7k) k Cor reference 
to figures 62 through 65). The derivative spectra exhibit a subrange 
behavior similar to that shown by the velocity and temperature spectra 
of figures 62 through 65, with the MITOS I temperature derivative 
spectra exhibiting an extensive inertial subrange and the OWAX 4 

E iste derivative spectra having only a narrow subrange. The 
Bo values from the MITOS Iderivative spectra are comparable to the 
` implied BE values. The MITOS III (Table 8) and OWAX 4 f j values 
are much lower (by a factor of 3) than the implied f 0 values, The 

1 ET values determined by the direct dissipation technique do not 
require the assumption that production is equal to dissipation. The 
direct Bo values are mainly dependent on the shape (area) of the 6 - 


۰۰٠٣٣۰٠ and the direct Xo and € values. Since the 6 -spectrum 








FIGURE 66 


MITOS I-l and 7 تفت‎ velocity derivative 

2 
spectra normalized with € 3 plotted lineàrly 
versus normalized wavenumber (nk) 


(ordinate represents o from direct 


E 
dissipation estimate) (from McConnell, 2e 


1974) 
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FIGURE 67 


MITOS J-l and 7 ice temperature 


derivative spectra normalized with 
-1/3 

X of / plotted linearly versus 

normalized wavenumber (77k) 


(ordinate represents Bo from direct 


dissipation estimate) (from McConnell, 


2972) 
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FIGURE 68 


-1/3 
OWAX-4 k Z velocity and temperature 


2/3 


derivative spectra normalized with € 


1/5 
/ (temperature) plotted 


(velocity) and Xo € 
linearly versus normalized wavenumber (nk) 
(ordinate represents o or Ba from direct 


5 
dissipation estimate) (from McConnell, ١ 1974) 
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for OWAX 4 is displaced to the right of the u - spectrum, a larger 
relative area is obtained to determine x which results in a smaller 
value of Bo as compared to the MITOS I O0 - spectra and Po values, 
13 
Gibson etal’ obtained a value of Bo = 1117 from BOMEX 
data using the same procedure as used above to calculate P. . At the 
time, the Bo value was considered high, and two possible explána - 
tions for the high value were suggested. The first explanation was that 
i ۱ : de 
an anisotropic temperature field (estimated by S T? » skewness of 
9) would enhance vertical temperature gradients if positive skewness 
values were observed, resulting in higher f, estimates, Skewness 
eo 
values for the temperature spectral derivative (from email E 
averaged 0.81 for the MITOS I data and 0.84 for the OWAX 4 data, 
not significantly different. Both skewness results indicate similar 
anisotropic conditions. The second explanation was that Ba would 
increase with increasing intermittency of €, estimated by the factor 
1/3 


G K) , where K is the kurtosis of du/dx. Kurtosis of du/dx 


(also from McConnell)" from MITOS Iand OWAX 4 were ~19 and e 
EB. E A nn E: w did not A ee to ۹ 
the n Bo values. The menci Bo values from this MS ro E 
= to the BN es of mc ka x Pad : 

to similar measurements over land (dry conditions), however, other 


effects besides anisotropy and € intermittency are required to 


explain differences in the Bo values from the different cruises, 
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Results of calculation of momentum flux and  cnsible heat 
Nb the direct dissipation method for MITOS I, MITOM III, and 


OWAX 4 are shown in Table 9 for comparison to the directly measured 





TE 2 
covariances. "The direct dissipation estimates of u, and wO were 


4 


calculated using Eqs. (7) and (8), rewritten to include corrections for 


stability effects as 





DCN S) (65) 








4 
D 
| 


(66) 


For Eq.(65) the two expressions for © given in ee "و و47‎ 06 
compared. The u, values used in Eq. (66) are the direct dissipation 
estimates for Z/L = 0. The expression for $ was used in Eq. (66). 
The direct dissipation estimates of d for MITOS I are 

all higher than the calculated direct covariances suggesting the assumed 
values of Cp DS 107? is low. If a value of CD z]1.5x 107? is 

sed anda 10% cup over-speeding correction is دا‎ in Section 

4, 7a) to the direct dissipation estimates, the results for MITOS l are 
A Pee o and OWAx A. I should bi nota O 
that the E gence limits for the direct covariances are +25% and 


for the direct € estimates + 20% so that comparisons between the two 


Methods will only refer to trends using thc values of Table 9. The 





26 3 








895 2 c8 'I UT T v XVMO 
16 *2 98 'I 097 "TI e 111 نله‎ 
Dj se 58*€ Ore L 
90 *2 ھ2‎ | 61 *1 1 I SOLIW 
lp S 8 5 5 0 3 RE 777 m کے‎ 
EET ہی‎ D 
.9 82 $ d 
Pe | = gm (0 = “| 2 | = gm om 7 
2/ ， LAN “azy] Oy , 


aa 








)5957/ 7۵ - Do) 13 9111 36 

















682 6S€ s che y XV AO 
565 PIS 287 ۰5 2 TSOIN 
PIS 99 [I9 )506 TE 
£2 L88 818 (249) I ISOLIW 
» S. 
دک‎ aja gp y = Siem: noe 
€/2 az ee ê €/z $2 -— z 
3 5 LX | 
c/o! @/ 924) = z“ ) WO) xn[j umjuouioJA = 
. 6 


5 3251017153 3, 11013120755701 15510 puer 11011215302 Ápp^ Jo uo stredwon.: 1 


6 oLIGVv.L شر‎ 





264 


same is true for the wd estimates with better confidence limits on the 
values obtained from the two methods. 
a, ۱ 
Best agreement between the values of u, is obtained 
using z/L=0 for the direct dissipation calculations} Including the 
2 
MITOS I results corrected as suggested above, the u, estimates 
using v are an average of 12% higher than the direct covariances, 
1 
and the estimates using E are an average of 13% lower. The 
2 
results indicate that the dimensionless turbulent transport is, 
20015 
approximately half of that suggested by Wyngaard and Cote , and in 
ee 
agreement with the results of Leavitt for turbulent transport from 


4 


BOMEX measurements. 

The direct dissipation estimates of w6 are significantly 
higher than the direct covariance measurements, If the 10% cup 
over-speeding correction is applied to the MITOS I wO estimates, the 
direct dissipation results for all the data are higher than the direct 
covariances by an average of 30% for z/L=0, and 77% when 
Btability corrections are included. Differences between the wO values 


for the two methods cannot be accounted for ۰۰۰٠٠٠٠۶ ۲٣۹۹۱١٢168 the 


pitch corrected covariances and the Xo and € estimates, 


۰ ۰ : 
۰ .. s 5 کا ان‎ ۱ ۳ 2 un RB ee, - 
e LA Qvid 50 > » ^ pod ا‎ A 4 wy? ت‎ 4 Se yo ۲ 1 : ۰ .و‎ 


cularly when stability corrections are included in the direct 


es 
۰ $ 


parti 
dissipation calculations. The results support the evidence from the 
previous section that dissipation of temperature variance may exceed 


mroduction by a factor of 2 or more. 
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14 

Stegen etal (1973) used the direct dissipation technique 
to estimate momentum flux and sensible heat flux from high frequency 
hot-wire and cold-wire velocity and temperature measurements. Using 
the stability expression e to OIT ect the direct dissipation estimates 

2 
of u, , Stegen etal obtain an average value of shear stress of 一 
e ۱ 
0.39 dynes/cm , somewhat lower than the eddy correlation results 
8 x 
meporeca by Pond et al of 0.49 (Univ. of British Columbia result) 
2 l 

and 0.44 (Oregon State Univ. result) dynes/cm . If the results of 

14 ۱ 
Stegen et al for 6.7 m. (-z/L* .26) were corrected for stability 
effects using the expression D (Eq.(57), the estimated stress would 

1 
2 
be 0.39 dynes/cm و‎ in better agreement with the direct covariance 
0550115. It should also be noted that the direct covariance results 
Were corrected for pitch using the R ae 0150۹۳۲601۱11100 to estimate 
u 
the rotation angle. The analysis of Section 4.3 suggested that this 
method may tend to overestimate the correction under certain 
conditions, however, not necessarily applicable to the BOMEX results. 
14... m 

Stegen etal did not apply stability corrections to the  . 0 el 
Birect dissipation estimate of the sensible heat flux since the 
variability in the o ‘estimates were considered greater than the errors 
due to neglecting stability effects. As a result, a value of He = 0, 74 

2 
mw/cm was obtained, lower than the results reported by Pond et al 
: | 2 ۱ : 

ERO (OSU) and 1.3 (UBC) mw/em from direct covariance 


measurements, also using the R= - OVSmcowecciion, If stability 
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effects are included in the ealeulation of Stegen et a (using 6.7 m 
end - 2/L = 0.26), the resultant He AS u , in agreement 
with the direct covariance measurements. This result is much lower 
than the OSU inertial dissipation estimate of He - 2.68 ری‎ using 
a Bo = 0.4; for Bo = 1, Better agreement is obtained with a 

corrected value of He = 1,1 for the OSU inertial dissipation estimate 


of the sensible heat flux. Although there appears to be agreement 


between the direct dissipation estimate of H, when stability effect 


o 
are included and the direct OE T: measurement, the absolute 
spectral levels for the OSU, UBC, and UW (Univ. of Washington - 
p temperature spectra did not overlap with the low frequency 
portion of the UCSD (Stegen et al S temperature spectra (personal) 
communication, Leavitt, Friehe). Differences between the two sets of 
spectra may have been due to spatial resolution effeets at the high end 
of the low frequency spectra (OSU, UBC and UW), and high pass 
filtering effeets (below — 2 Hz) at the low end of the high frequency 
(UCSD) temperature spectra, 

Zes'was mentioned earlier, the MITOS 1cruise had conditions 
similar Bose of BOMEX, and the high frequency temperature spectra |. ND 3 
و ا‎ used to obtain Ke estimates were matehed to the 


temperature spectra presented in Section 4. 5b so that any differences 


Aue to calibrations would be eliminated. 
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٠٣3167۲6 ۶۱18.1 6. 6101062608 recxamıned 
A basic assumption common to both the incrtial and direct 


dissipation techniques uscd to estimatc the fluxcs of heat and water 


| 


vapor, is that production of tempcraturc and humidity variancc is 


equal to dissipation (Eq.(4)). For opcn ocean mcasurcments, this 
assumption may not bc valid. Production and dissipation of 
temperature variancc were calculated for the cold wire runs and the 


۱9۳ ۱۱۳۲5 are given in Tablc 10. The 7 values obtained by McConnell 


(1974)? and given in Tablc 8 wcrc used to calculate the dissipation. 


Production was calculated from the relationship derivcd by combining 


A 


Eqs. (4), (6), (8) and (19) and written as 
à 
= To our 


dT ee 
-wO. dz 0 17 Pry ١ 5 








oimilar results for humidity were not obtained since X, could not be 
determined from humidity mcasurements, duc to the limitcd frcquency 
response of the Lyman-Alpha humidiometer. Calculations of production 
of temperature variance wcre bas ed on 4 = 1. = K |) ^ kiz 0. zu 
Er 7 (1-16 2/my 2 ۲ en T. aad ۳ E و یں سے‎ 
uui سک ایقوسلوعفمش٭خصہ ےہ امو‎ 
dis sipation to production CEPS given 0118 ٢ 

For all the cold wirc runs, dissipation of temperature 


variancc is grcater than production. The smallcst ratios for 


dissipation to production from MITOS I of 2.6 would bc incrcased to 
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۱۳۰ “+0 


Production and Dissipation of Temperature Variance 


Dissipation 


2 2 O TIC NEN ری‎ 
Production (°C /sec) Dissipation (°C /sec) Production 


4 4 








MITOSI-1 1.7 ۶ 0 DEC OS mm 
E ED 
E usc 0 ۲ 4.3 x 10 2.6 
` -4 8980 
IEAUPOS H2 0.83 x 10 Be 3.8 
-4 e 
OWAX 4 2.9 x 10 1.3 x 10 4. 6 
2 4 
= T u 
— d ماج‎ sk mn 
Production = -w6 7 = (IO AA / 
dz aKz 


.. Dissipation =. 
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about 3 ifa 10% cup over-speeding correction is included in the 
T culation of xe . The values of production for the MITOS I runs 
would not be significantly affected by an increase in ue if a value of 
Cp کک‎ 107? were used, since e, would also increase, and 
the increase in $ and u, would be balanced by a comparable 
decrease of "e 1 

Although dissipation of humidity variance was not directly 
measured, results of Section 4. 7a indicate that production of humidity 
variance may be larger than dissipation. Using a value of D -0.4, 
and S1 1 2 , 1 obtained es production and dissipation 
rates for humidity variance. The average implied Ê لمت‎ [rO 
comparison of the direct covariances and estimates of the inertial 
dissipation technique, suggests that for a given - 0.4 , production 
is greater than dissipation. If the expression for D Dr and 
EN => Pu as suggested by 0 flux profile analysis, the 


E 


implied B. would be more nearly equalto 0.4 if production 


equals dissipation. 


cono 117670177 61 2 inertial dissipation IC for 
۳ temperature and the inertial dissipation results for humidity, ao, | 
Reon sideration of the as E nes T eee the c "oem 
budgets seems appropriate. Eqs. (29) and (32) were obtained by 


assuming horizontal homogeneity and stationarity. To test the 


assumption of horizontal homogeneity for the temperature variance 
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budget, ! term - u8 = Mas lS De considered. Forthe OWAX 4 
run, a horizontal temperature gradient of approximately 5? C per 
100 m. would be required to balance production and dissipation. 
This gradient is far greater {han would be expected 7 the ocean or 
seen in the OWAX 4 time 510001773077 "58661166 41۰ 0 non- 
stationarity may have affectcd the temperature results for the SOMA 
cruise due to the small heat fluxes measured, non-stationarity effects 
for the other cruises can be expected to be small as noted by 
Wyngaard and Cote pil for ام‎ temperature measurements. 

: Two terms which appeal im the temperature variance 
budget and not included for humidity are the heat transfer due to 
radiation, 20R', and the heat transfer due to viscous heating, 
28/6, . Effects of the terms were considered in Section 2.4. From 
the previous discussion it appears that the viscous heating term may 
be neglected. Effects of radiative heating/cooling for BOMEX 
conditions have been discussed by m ced Phelps and Pond, aC and 
Bond ct al. d Phelps and E suggest that absorption of longwave 
radiation by the increased water vapor concentrations in BOMEX would 
1 tend to eee a more homogeneous EI Structure and cause +> 
suppression of large scale temperature fluctuations. This would not 
explain the larger differences between dissipation and production for 
the Tess hunid OWAX cruise. As mentioned in Section 2,4, maximum 


一 = 一 一 ~ 


2 سد‎ ! 
destruction of 8 occurs at small scales (dissipation scales) for 





42 bs 
: If 20R were 


internal radiative transfer (Townsend, 1958) 
included as a dissipative term for this case, the difference between 
production and dissipation would be increascd. 
E qe 
Leavitt suggests that the cffects of cooling by radiative 
flux divcrgence mcasured during BOMEX may account for differences 
betwecn the temperature and humidity signals, the spectra and the 
or M ۱ 
cospectra, Leavitt also suggcsts that the coolcr air at heights 
betwcen 50 and 400m (as evidenced in the BOMEX potential 
temperaturc gradient) due to radiative cooling, descends by negative 
buoyancy adjacent to the backs of warm plumes and causes the "cold 
Spikes", As shown in figures Y and 9, the same phenomenon (cold 
spikcs) are also observed for both the warm, humid conditions of 
MITOS I, and the cold, less humid conditions of OWAX 4. From the 
OWAX 4 production-dissipation results, radiative cooling alone could 
not account for the differences between production and dissipation for 
both types of conditions. 
‚Ihe scalar flux divergence terms for the temperature and 
humidity budget may amount to 25% and 10% of the production as 
B | | 8ھ‎ 230 
~ Suggested by the results of Leavitt, , and Wyngaard and Cote 
(tempcrature). The magnitudc of temperature flux divergence is not 
sufficient to account for the differences between production and 


dissipation. If a humidity flux divergence of 10% of the production is 


assumed, it would rcsult in a decrease of the implied ib values and 


lend more support to the suggestion of Poco that T مر‎ - © , 
m 


; 4 

Temperature profile results of Paulson ct al implied Pa = 0,4 
Bez ۱ ۱ ۱ 

۳16 z/L) . This expression for ze would result in a greater 
difference between dissipation and production than obtained using 

-1/2 رر‎ 
D = (1 - 16 z/L) , to balance the production and dissipation 

: 0 

results would require 2 ۶2:401 516۷9271 [( , highly unlikely. 

The remaining terms in the scalar variance budgets are the 
source/sink terms, Sn and 0 mel represent contributions or 
losses to the variances due tv interaction or possible effects of spray 
on the temperature and humidity fields. It is proposed that the 
evaporation of spray can account for differences between production 
and dissipation of temperature variance. Evaporation of spray in 

: ۱ ; 51 
conjunction with the ramp modol proposed by Gibson etal can 
explain the ramplike structure of the temperature and humidity 
time traces, and can also explain the cold spikes in the temperature 
signal and the associated behavior of the humidity signal. Upon careful 
examination of the temperature and humidity time traces, they may be 
represented as shown in figure 69. The associated vortex structure 


-shown below the temperature and humidity signal representations is 


TI S. £s 
NE v.v 
... 


E 5 ۹ H مات‎ 2. . an t t aS E کن ادا واه د‎ O A 
و‎ q us, FU سرت‎ - » E 3 ^ . 2 et aus ٠ 22 ۰ “A ¿en ۰ 5 ۰ E = 3 e ` ۰ x um .. a ۰ ۰ - ^. ‘ot + Jy rs ne "a ue e. 


e 5 : ۰ re " is ۱ à P l ۰ 51 Mr E FAS 
an end view of the horizontal roll vortices of the Gibson model. The 
wavelength of the ramps is shown as 100 moters representing the 
large scale ramps in figures 7, 8 and 9. The 5m. ramps shown by 


: 51 ; ۱ i . 
Gibson etal are not inconsistent with the fine scale structure shown 











FIGURE 69 


Ramp model for temperature and humidity 
over the ocean 
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in the magnifications of figures 7 and 9, 

Water droplcts of spray accumulatcd in fine sheets along 
Pvesintetfaces of the roll vortices evaporatc, causing the ''cold spikes" 
of the temperature signals. The additional source of evaporation 
which incrcascs thc tcmperaturc variance cau account for the difference 
bctwecn calculatcd valucs of production and dissipation of tempcrature 
variance, Since the air above the vortex is dryer than the air within 
thc vortex, cvaporation occurs on the outcr edge of the vortex cooling 
thc overlying air. The maximum accumulation of air cooler than the 
ambicnt temperature will then bc forced down at the leading interface 
of D vortex due to the negative buoyancy of the coolcr air and thc 
EE eburc ofthe roll vortex, The désccnding coolcr air at the leading 
edgc of thc vortcx causcs a cold spike, and the excess water vapor 
duc to evaporation causes a "rounding" of the corrcsponding humidity 
Minima at the intcrface. The small amount of "smoothing" of the 
humidity signal at the lower part of the ramp interfacc may also bc due 
to the Lyman-Alpha humidiomctcr mcasuring a combination of water 
Mor ana spray oen has na thc vortex lcading edge or 


Bifecis offrcqueney rcsponsc limitations. The amount of water vapor 


' requircd to account for the difference bctwcen production and 


dissipation is small due to thc high valuc of thc latent hcat of vaporiza- 
tion, Similarly, thc amount of water droplcts required would also be 


mall Differcnccs betwecn the tempcraturc timc traces for the two 
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cruises will be discussed in detail in Section 5. 
| : 70 

From a laboratory study conductec by Wu (1973) , the 

concentration of water droplets measured at velocities above 8 m/sec 
2 

indicate that less than 10 droplets per m (mean diameter ~ 150 p ) 
occur at heights between 5 and 20 cm in a wind wave tank (14 mx 
1.5 m) for wind velocities less than 8 m/sec. Although it is fairly 

70 7] 
well established (Wu, 1973, Krauss, 1967  ) that droplet concentra- 
tion decreases exponentially with hcight above the sea surface, absolute 


droplet concentration for a given height and wind velocity is not known. 


For the OW AX 4 run, a value of So = ~ |x ome ? C/sec 


4 


is required to account for the difference between production and 
: 3 
dissipation. If S (ugm/cm - sec) is used to represent the adiabatic 


rate of evaporation, then =: may be written as 


E. MN 
» Ss. Tue 68 
o“ E 6) (68) 
P 
where S = SS (the mean and fluctuating evaporation rates). I it 


and @ are perfectly anti-correlated 


ne‘ 


is assumed that 
S : . 
(r* = --1, an increase of S results in a decrease of 0), ۱ 
۰ ES 9 e O 1 0 x ۱ Me m ۳ 7 3 ات مت‎ ÉÉ a we. n 0 رت‎ N "m A i zi = E Ve 0 a gs 


then ME may be estimated from the relation 


D 
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-3 2 
BNO Y So => 1.0 x 10 °C /sec and Oo = 0.1 ات‎ value of 


-3 3 
a. = 4.0x 10 pgm/em' sec is obtained. To determine if this value 
of O., is reasonable, an examination of the time series is required, 
Ifanaverage AT of 0.1°C for the "eold spike' portion of the 


temperature ramp interfaee is assumed, the time required for 


evaporation to oecur may be ealculated from the relation. 





Hoe (70) | 


obtained by integrating the time rate of change of the speeies 
concentration equation (for Tp 0 - sensible heat) and assuming the 
evaporation source term produces the entire temperature change 
0 ۰ z -3 3 
OI C), ore gs 2 c. ^ 4x10  ugm/em - see, the At 
required is ~ 10 see. For the OWAX 4 run (U — 5 m/sec) this 
represents 50 meters, approximately half a ramp wavelength. For 
B. 2i <l or AT > 0, 1 the time and/or distance would be larger. 
S 

Assuming that evaporation oecurs along a considerable portion of the 
vortex interface, this result does not seem unreasonable. The resulting 

M "—- -2 3 | DS 
۵ و‎ due to evaporation is 4x 10 pugm/em , much smaller than the 
E nean water vapor concentration and yet possibly enough to causea  . 
smoothing of the humidity minima at the ramp interface, As mentioned 


previously, due to effects of frequeney response limitations of the Lyman- 


Alpha and possible sensing of water droplets, this smoothing may be 





E mental. 

For the MITOS I and BOMEX conditions, the evaporation 
of droplets may also be occurring in conjunction with radiative 
cooling as suggested by Leavitt. a The combination of the two 
phenomenon could account for the imbalance of the temperature 
wertance budget for the MITOS I data and the difference in the 
temperature signals of the MITOS Iand OWAX meos Radiative 
cooling in the air above the roll vortices combined with evaporation 
would tend to break up the ramp structure and lengthen the "cold 
spike'' as shown in figure 69 and in figures 7, 7a, tal, and (a2 99 c 
Eme between the two time series can also explain shape of the 
temperature spectra and the wO and uO cospectra from MITOS I, 
The spectral and cospectral results are discussed in Section 5, 


4,8 Application OI RE Bulk ACO mamic TCDD ge for Estimation 
of the Turbulent Heat Fluxes 


The bulk aerodynamic method relates mean wind speed and sea- 


21۳ temperature and humidity differences to the fluxes of latent and 


sensible heat by the relations 


(71) 


ar 
1 
2 
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where AT is the difference between sea surface temperature and air 
temperature at a specific height, and Aq is the corresponding 
absolute humidity difference (saturated at the sea surface temperature). 
2 
Roll has suggested that the constants C - C = C for near 
w 6 wq D 
neutral conditions. To include the effects of stability on the mean 


wind shear, and temperature and humidity profiles, the expressions of 


(71) may be rewritten as 


CC 


S p wo 
١ (72) 
۱ H, L Gr U Aq = 
where 
(c 7/2) * & 
PB +٠٣ 2153) : (73) 


The constants 0/4 and é result from linearization of the 
dimensionless scalar gradient TaS 9 ‚ and the dimensionless 
wind shear ® . The constant 2 ی‎ represents the "roughness" 
ا‎ a Summed equal to one tor the — RS 
calculations, Using values of 0 = 16 and € = 3.5 (Krause, 
MC oce eanetsky (1962) ^), "valuss U ATS S voci 
and U Aq 7 were calculated for each run and are given in Table 11. 
Also given are the corresponding w and wq Ao ues by the direct 


EDvariance technique, aod the C and C required for agreement 
Ww wq 


6 
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FIX to estimate the total heat flux. 


between the bulk aerodynamic method and the direct estimates. 

The average value of er is (1.29 + 0.36) x 107? > in 
agreement with the value suggested by Pond et 本 لا .۲ » و9‎ re) 

-3 
x 10 for BOMEX results. For the pitch corrected covariances 
(MITOS III and OWAX), C = I03, suggesting "Qo «NE 
wa wd D 
ee 8 

Similar to the result of Pond et al for pitch corrected covariances 
using estimated pitch angles, 


The large range of the oae values suggests that the bulk 


0 
aerodynamic method is not useful in estimating the sensible heat flux. 
Trends in the data indicate that for warm humid conditions the bulk 
ES ue method over-estimates the sensible heat flux, aud under - 
estimates the flux for colder, less humid conditions. Similar results 
were obtained by Pond et T with regard to lack of agreement between 
the direct and bulk aerodynamic techniques, however, the sense of the 
disagreement was reversed. 

Results of this work and those of Pond et ul indicate that 
estimates of total heat flux (sen sible B latent) based on one flux 
E سس‎ by the bulk بد کے سس سب‎ c an average BOWEN 


+ 
ratio should employ the latent heat flux instead of 


Sam D XN Ese Ji.‏ و 
S .. 34 ee e TON PF ۳‏ 
ا یت 
۰ 


the sensible heat | 


aod 





5, SUMMARY AND CONCLUSIONS 


The primary objectives of this work were to measure the 
turbulent fluxes of momentum, sensible and latent heat by the direct 
covariance method corrected for instrument motion, and compare the 
direct measurements to estimates of the fluxes by the inertial and 
direct dissipation techniques. Stability corrections were applied to the 


estimates of the fluxes by the dissipation techniques. 


` 
۳ Momentum Flux 


Directly measured uw covariances, corrected for internal 
sensor alignment and instrument tilt using measured pitch angles, had 
an average uncertainty of ±4 2595. Sources of errors were discussed 
in Section 4,3. The corrected covariances were compared to the uw 


values which would have been obtained using two different methods of 


estimating the mean pitch angle in lieu of actually measuring the angle. 


Results of the comparison suggest that for small pitch angles («2?) 


"With small variations کو‎ 0 5 0( ` Corrections based CG O~- Siùn 


W /U tend to underestimate the pitch correction, and those based 
m mnm ۰ ۱ 


e. > 
te. 


kon R y --0.5; for 0.01«f£«0.1 Hz , tend to overestimate.the WT AA 


pitch correction. The value of the spectral correlation coefficient for 
the measured pitch angle corrected covariances was Roy E) = =- 0.32 + 


0.08. For the MITOS IIl uw covariance with us > Q0 xc cies 


instantaneous pitch correction applied to the data gave a uw value 
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comparable to results of the other two methods. The estimated pitch 
angles, however, were higher than the measured pitch angle. 

The measured UNAM - 2.8 £0.5 , was comparable with the 

x nd: ۱ 
results of Leavitt from similar open ocean velocity measurements 
۱ se 
using estimated pitch angle corrections, and Mitsuta and Fujitani 
using inertial reference frame corrections for ship motion. Pitch 
D 
Corrected O My results followed the predicted (- z/L) behavior 
۱ rl 
and were comparable with the results of Leavitt ^ and Pond et a]. 
Values of ou, for this work were higher than the results of 
58 
Wyngaard etal obtained over land. 

Vertical velocity fluctuations over the open ocean may have more 
low frequency energy containing eddies than over land as suggested by 
the scalar ramp model and associated horizontal roll vortices shown in 
figure 69. Discussion of the time series in Section 4. 1 included the 
observation that warm, humid air was associated with upward vertical 
velocity and negative horizontal velocity changes, with the opposite also 
being true for colder, less humid air, The associated anisotropy of the 
velocity field due to the large eddy structure which determines the scalar 

1 ; : Dd 
field structure is supported by the results of Gibson et al, who found 

1 positive skewness values of about 0.5 for (0w/ bx) when conditionally TM 

sampled (i.e. when w~ 0), as compared to the isotropic value of zero. 


The average correlation coefficient for the pitch corrected results 


was r = -0.26+ 0.08, also in agreement with the results of 
uw 
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zen 8 
Leavitt ^ and Pond et al. lo remove possible stability dependence of 


Be correlation coefficient, the tr values were normalized with the 
UW 


expression for the non-dimensional dissipation of kinetic energy. 


23 
Using ® M = (2 -z/L), the resultant normalized |value was 
1 
E. 9 2 
r --0.25 +z0.06 , comparable to the result of Miyake et Tn for 


UW 


near neutral conditions. 
Velocity AN Spectra, memnalizew toduclude stability ۵ 
Ba TT 2 
agreed within the given uncertainties of the u, values with the results 
ED ; 56 -2/3 
of Leavitt and Kaimal etal, and followed the predicted fN 
inertial subrange behavior. Values of 0.29 +0.06 (horizontal velocity 


spectra) and 0.40 + 0.09 (vertical velocity spectra) for the کپ‎ - 


spectral constants were in better agreement with the predicted constants 


using the dimensionless dissipation expression for E (as above). 
1 
Lower values of the constants were obtained using the form of 
273 273 oT 
M 5 ۹ ۵ B z/L| / ) suggested by Wyngaard and Cote.” Scatter 


2 


in the spectral results for the first two cruises was attributed to x- 
wire drift, possible cup over-speeding errors, and uncertainties in the 
. ۱ ا‎ == ms P ا‎ M dcs کہ ار‎ io: a0: 
calculated u, values from an assumed رن‎ Sle). Spatial 


averaging effects and distortion of the vertical sensing path of the 


و ^ 
5 


sonic anemometer (OWAX) caused a steeper roll-off in the high 
frequency end of the u and w spectra for the last two cruises than 


the first two cruises using x-wire anemometers. 


The normalized pitch corrected momentum flux aN - cospectra 
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1 
werc similar to uw cospectra obtained by Leav.. and Pond et al. 
The uw cospectra subrange constant for this work, 0 - 0.01, was 
uw 
1 
lower than comparable estimates of 0 5 by Lcavitt : (0.07) and 
u y 
56 
Kaimal et al (0.05m The lower Q ٢ value may be duc to à combina - 
u 
tion of high frequency noise in the uw cospectra and possible distortion 


of the sonic anemometer (vertical path) as mentioned above. Howcver, 


-4/3 


the high frequency bchavior of the cospectra followed the predicted SN 


subrangc slopc over nearly a dccade for fN rl. This indicates tie 
lowcr Bm value may be rcal and suggests that production of kinetic 
energy may occur at larger scales ovcr thc open ocean than ovcr land 
under certain conditions, A shift of production of kinctic energy to 
larger scales is consistent with the ramp modcl, with well dcfined 
horizontal roll vortices for thc colder, less humid conditions as shown 
in the bottom portion of figure 69. 

Comparison betwcen estimates of the momentum flux by the 
direct covariance and inertial dissipation techniques was made using a 


3 
E spectral represcntation of the data, illustrating the approach and 


5/3 


- 


agreement with the predicted k- power law behavior required for'* 


application of the inertial dissipation technique. Agreement was 


+ 
A icing NE 
T 


و 


obtained between direct covariancc estimates of momcntum flux, and v" " 
inertial and direct dissipation estimates within the uncertainties of the 


direct estimates, applying the assumption that total production of 


kinctic energy (mechanical + buoyant) was equal to dissipation (2 hs 


1 
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E. 
mower estimates of momentum flux were obtained by the dissipation 


techniques using the expression suggested by Wyngaard and Cote 


(P D. The CIEE value Ol thE velocity jmcrtial subranpe'constant 
a 


Nom the direct dissipation tcehniquc was œ = 0.53 +0.03. Trends 
u 

in the direct dissipation estimatcs of momentum flux suggest that 

Buvrbulcnt transport cannot be ncglcctcd in thc marine boundary layer 


and may bc approximately half of the estimates of dimensionless 


3 
turbulent transport ovcr land (Wyngaard and Cotc D 


E. Batent Heat Flux 


The pitch correctcd wq covariances had an average uncertainty 


4 


of ~+ 10%. Larger uncertainties could be attributcd to the 
Encorrected wq direet cstimates which were probably determincd 


within +25%. Values of ul, for this work followcd the predicted 


-1/3 


۲ z/L) behavior with e = 0,95, Values of 0 were lowcr 


than corrcsponding D. valucs, and werc lower than the cele 
10 . 1 

results of Phelps and Pond and Lcavitt for BOMEX data, The 

average correlation eoeffieient for the water vapor flux results was 


r SORA 20: 11 with a trend in the r Values Trom the colder, 
wq wg 

e 
liess humid cruiscs to bc higher than results from the warm, humid 


b S 
مب و‎ TOS : : 


cruises. To removc trends due to stability effccts, normalizcd 


Pe ۵65 


Borrelation coeffieicnts (r J) wcrc ealculated, using C & 
wq - H 6 


thc normalization. The trend remaincd in the r valucs and may be 


duc to the better defined ramp and associatcd vortex siructure shown at 








Cor 


the bottom of figure 69 for the colder, less humid conditions, as 


compared to the warm, humid conditions. Values of r 8 and r 
wW uq 


for the warm, humid conditions were comparable to results of Phelps 
10 
and Pond from BOMEX measurcments. 
Normalized humidity ÍN - Spectra and watcr vapor flux cospectra 
collapscd well and were similar to velocity spectra and cospectra. 

TT -2/3 
Eumudity spectra followcd thc prcdicted SN behavior. Normalized 
wq and w8 cospectra from the colder, less humid conditions werc 
also similar, rcflecting thc similarities in the tcmpcraturc and 
humidity ramp structure for those cruiscs, as comparcd to the warm, 
humid conditions. Normalized water vapor flux cospcctra exhibited an 

-4/3 ۱ i 
IN subrangc behavior with an avcrage value of thc cospectral 
subrange constant of g a a1) (le se 10h Or 8 joven 8> 
7 | 

the OWAX valuc of y “0.03. Comparison of a and Q values 
wq wq uw 

suggested the vcrtical transport of watcr vapor is similar to transport 

ES momentum for small and large scales, cspecially true for the 

colder, less humid conditions. Normalizcd horizontal water vapor flux 


EU 


Bospcctra exhibitcd an e subrange behavior with a trend in the 







a values. suggesting relativcly more cfficicnt transport of horizontal 
water vapor flux at small scalcs for warm, humid conditions. This 

| result is also consistent with the ramp modcl shown in figure 69, where 
ihe aller scales are more pronounced for the warm, humid conditions 


‘(upper portion of figure 69). Comparison of the m and 0 results 
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indicatcd vertical transport of water vapor is always more efficicnt 
than the horizontal transport for smaller scales. 
۱ 5/3 
A normalized k spectra representation was uscd for 
comparison of latent heat flux estimates by the direct ee and 
inertial dissipation techniques, Due to uncertainty of the actual value 
of the humidity subrange constant i , the spectra were plotted to 
epresent the implicd requircd for agrcemcnt between the two 
estimatcs of the latcnt heat flux. An avcrage 8 = 0.21 ± 0.05 was 
required for agreement. The direct dissipation tcchnique was not 
applied to estimatcs of latent heat fluxcs since the frequency response 
of thc Lyman-Alpha humidometer, used to measurc humidity 
fluctuations, was not sufficient for direct estimatcs of dissipation of 
humidity variance. The implied ß was lower than B -0.41 obtained 
q q 
65 
by Pacquin and Pond from measurcments of 2nd and 3rd order 
E one T 
structure functions. For ha =0.41, the k humidity spectra results 
suggest o = = for production of humidity variance equal to 


dissipation, or production greater than dissipation for ®_ = 9 nahe 


E 
EN os n ک0"‎ ; E CS 1 men 11 
a for this work is comparable to 3 = 0.25 suggested by Leavitt 

4 
based on the flux-profile results of Paulson et al, و‎ and also comparable 
کت‎ E F رات‎ E DS E E >“ . ۲ TZ . ۲ تن مد‎ 
to the results of Smedman-Hogstróm (1972), who obtained a value of 28 # qmn n. 


B -0.292 0.10 from measurements over land in Northern Sweden, 


assuming production equaled dissipation. 








o Sensible Heat Flux 

Directly measured w6 mE cci covariances, corrected for pitch 
effects, were determined with uncertainties of ~ +10%. Uncorrected 
direct wô values for the SOMA cruise were subject tb x-wire drift 
(in w) , temperature calibration uncertainties, and BO ble effects of 
non-stationarity for near ESI stability conditions, The average 
SOMA and MITOS I direct w9 values are probably estimated within 
+ 25%, although this uncertainity may be higher for the SOMA 
cruisc. Values of G/T. for the MITOS 1, MITOS IlI and OWAX 


-1/3 
cruises follow the predicted (- z/L) u behavior with CT - 0.95 


as suggested by Wyngaard et al. i me ہہ‎ values for the warm, 
humid cruises are higher than 37 2۳ and more comparable to BOMEX 
results of Leavitt, x and Phelps and ome than the results of the 
colder, less humid cruises. High values of og Fs for the SOMA 
cruise were attributed to uncertainties in Ue values (from an assumed 
C and possible effects of non-stationarity. 

D. Significant differences were exhibited between average values of 
the MES ion ٠۴ 全 ۲ for the cr. less humid ON 


Borrected for pitch (r qu 0.48 4 0.04) , and the warm, humid 
le 595 e at ist 9 T e E ss 3 AN . ie EE 2 5 د یش‎ Y a Pai iiis as was WI rt. se. : 


ies nee 


ms ۰ 


Bruises (0.16 20.11), not accounted for by lack of pitch corrections 
Or the warm, humid cruises, The r 0 corrclation coefficients werc 
۱ wW 
normalized in the same manner as was done for the r ۲ values to 
NV 


remove any trends due to stability effects. Normalized correlation 
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eoeffieients, r exhibited similar trends as the non-normalized 


w9 


values, however, less pronounced., The higher values of r o and 
wW 
ما‎ for the eolder, less humid eruises ean be attributed to the well 


defined ramp strueture for those cruises, illustrated in the bottom 


portion of figure 69. The average value of the normalized correlation 


coefficients was 7 = 0.21 + 0.21 , eomparable to the result of Miyake 
5 ۰ 

et al Y Ivo =0.24 z0,2) for near neutral conditions. No significant 

trend in the Y 6 values werc Observed, 


Normalized temperature 3 - Speetra from different cruises 
with different meteorological conditions did not eompare with one 
another, and were not similar to humidity or velocity spectra. 
Normalized sensible heat flux eospeetra exhibited similar differenees, 
however, the wO eospeetra from the colder, less humid eruises were 
Somewhat similar to the corresponding wq eospeetra. Temperature 

- -2/3 

spectra from the warm, humid MITOS I eruise exhibited a f 
behavior and were similar to temperature speetra obtained by Leavitt, 

10 E i 
and.Phelps and Pond for BOMEX data. Temperature spectra from 
Ehe eolder, less humid eruises appeared to follow a -1/4 power law 
-. behavior over nearly 2 1/2 decades of normalized frequency above 
E 7 0,5 ." The normalized temperature spectra from the warm, humid “SHE 
Eulises in general exhibited higher spectra levels than the speetra from 


eolder, less humid eruises. Ifitis assumed that evaporation of spray 


is eontributing fo the suppression of low frequency temperature 
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Mimetuations (discussca in the next section), spectra obtained at heights 


fart r from the surface should have higher normalized low frequency 
5۱6 ۰:1 ۱6۷۵15 than spectra from lower heights. Evidence of this 

5/3 
cffect can be observed by comparison of the k temperature spectra 


Bere 63) from MITOS II (12.5 m,. ) and OWAX (3.5 for the 
momen, less humid conditions, and comparison of 8m. and 30 m. 
ا کے‎ 
tempcrature spectra obtained by Leavitt. In both cases, spectra 
obtained at greater heights above thc ocean surface exhibit more 
mE 
paonounced low frequency enérgy. Leavitt suggested that the higher 
NS frequency 30 m., spectral levels may be due to radiational heating 
of the thermocouples, however, he noted that not all periods of day - 
light showed the low frequency fluctuations, 
cu "n -4/3 ۱ ۱ 
Normalized WO cospectra cxhibited an AN behavior, with 
higher values of « o associated with warm, humid conditions, The 
W 
shape of the normalized w8 cospeetra andthe wa values indicated 
feat production of temperature variance is shifted to smaller scales. 


Sifor the warm, humid conditions, and that under these conditions the 


‘Vertical transport of hcat can be 6 times more efficient than 


ale i 


Minomentum for smaller scales. Thc tcmperature ramp structures | 
E Buen ie. en ER 9 A np DEE ee 
represented in figure 69 indicate more small scalc activity for the 


warm, humid conditions supporting the cospectral evidencc. 


^ 


Normalized horizontal heat flux cospectra exhibited differences 


similar to the w8 cospcctra. The ud cospectra followed the predicted 
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-3/2 
M 


N inertial subrange behavior with values of «a a which indicated 
u 


that horizontal transport of heat was also shifted to smaller scales 
under warm, humid conditions. The ''cross-over' point of the wY9 , uO 
and Oq cospectra was that value of ER above or below which the 
cospectra from the warm, humid conditions exhibited higher (> 2 
Ber tower (< لہ‎ cospeetral levels than cospectra Trom the colder, less 
humid conditions. The value of EN at the cross-over point represents 
approximately 50 meters, which corresponds to half a fare wave- 
length as shown in figure 69, "For larger scales than 50 meters 
temperature is less correlated with velocity or humidity for the warm 
humid conditions, and thc rcverse true for scales smaller than 50 
meters. The "cold spike" of the warm humid conditions is not as 
Ensrp'es the cold spike!’ for the colder, less humid conditions, 
reducing the similarities between temperature and humidity (and velocity 
presumably) for scales larger than half a ramp wavelength. Smaller 
Beales im the tempcrature field are more pronounced for the warm 

^" humid conditions as illustrated in the top portion of figure 69, associated 
with a "brcaking-up'' of the roll a as shown. The break-up of the 

Eroll Aces 5 P aii ramp Structure may be due کے‎ it 
effects of evaporation (less pronounced for warm, humid conditions) and 
descending overlying cool air due to radiative cooling as suggested by 
Eus for BOMEX conditions. The 0 end O 6 results also 


suggested that vertical transport of heat is always more efficient than 





Aio 


56 


the horizontal transport for smaller scales as found by Kaimal et al 
over land; the reverse was also true for larger scales. 
Comparison of the inertial dissipation and direct covariance 
techniques for estimating sensible heat flux was made in the same 
manner as was done for latent heat flux. For agreement between the 
two techniques a value of the inertial subrange constant for temperature 
of e was required. "The s temperature specira indicated 
that under colder, less humid condition over the ocean, temperature 
may not exhibit an inertial subrange (slope -5/3) for scales correspond- 
uto k < 0, 3 sp? . Temperature measurements should be made out 
to scales as smallas 6 cm in wavelength when applying the inertial 
Bsssipation technique. As mentioned above, lack of an imertial subrange 
mre veimperature, under cold, dry conditions may be due to evaporation 
of spray. The inertial subrange of temperature is most pronounced for 
the warin, humid-conditions, even with a shift of larger production 
scales toward the diffusive scales. The subrange only extends to the. 
point where the production scales of the associated we cospectra 
E. eas. humid BEL no longer are dominant (wo co OT levels | 
ES to zero). 7 DN E Oe A يي‎ al NO ۲ 
Estimates of wê by the direct dissipation technique were higher 
than direct covariance estimates by a factor of 2 or more. 
Corresponding estimates of the inertial subrange constant for 


temperature, f varied from 0.7 £0.07 to 2,440.2. From 


۵ 3 
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direct estimates of the dissipation of temperaturc variance (v. /2) 
p Xg 


and sensible heat flux (w0), with uncertaintics of +20% and + 10% 

respectivcly, dissipation of temperature variance was found to bc 
ul 

higher than production (- w8 me by as much as a factor of 4. The 

differenees could not be accounted for by uncertaintics in the direct 

estimates, or by temperature flux divergence which may be as much 

0 ۱ 1 ۱ ۱ 
as 25% of production (Leavitt ). Diffcrences between production 
and dissipation of temperature variance could be accounted for by 
consideration of sources of production usually negiccted in the 
tempcrature variance budget and attributed to the combined effects 
of ocean spray evaporation and radiative cooling. 

The high B. values from the direct dissipation technique (which 
does not requirc the assumption of production equal to dissipation of 
EEDerature variance) may be due to the ef£fcct of unternal radiative 

az 7 NN 

transfer (Townsend  ) which tends to destroy 8 at diffusive scales 

or spray droplcts ( more likely for the warm, humid conditions with 

‚higher Bo values) hitting the cold wire probes and reducing the high 

frequency response. There is cvidencc that spray can hit probes since 

i X- -wire Rp the MITOS 1 I apc ا‎ cruises aus: mes rd 
were و 6ٹ‎ by sek nalen اہو ۱۳ وت‎ drift in the 
calibrations. Both of the above effects would reduce the high 


frequency area under the 6 - spectra and increase the cstimate of Ra : 


From the inertial and direct dissipation results, usc of the 
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dissipation techniques for estimation of sensible heat flux over the 
open ocean is somewhat tenuous due to the inapplicability of the 
assumption that production equals dissipation of temperature variance, 
and the possible lack of an inertial subrange. 


5.4 Model of the Atmospheric suríace Layer Over 





In Section 4, 7c a model of the structure of the atmospheric 
surface layer over the ocean was proposed consisting of Bois on roll 
vortices with axis aligned perpendicular (or nearly so) to the mean 
flow, influenced by the combined effects of evaporation of ocean spray 
and radiative cooling. The siructure of the roll vortices is determined 
by the characte ristics of the wind shear which brings up warm, moist 
positively buoyant air, balanced by the negative bouyancy of cooler 
overlying air due to evaporation and/or radiative cooling (inversion 
layer). For cold, > conditions, evaporation of spray droplets 
entrained in the vortex interfaces is enhanced, resulting in a more 
organized vortex structure, exhibited by the well-defined ramp 
structure in both the temperature and humidity signals. Under these 
conditions, the signals are more intermittent, with relatively more 
quiescent periods. Larger production scales E temperature and 
_ humidity variance are dominant. arhe poll vortices are stretched out 
E. the NE iuc Bus تن‎ coolerlr sabe vob Xu ٥ 
evaporation acts as a cohesive agent. Cold spikes are well defined 
and 7 short duration since only a relatively small amount of mixing 


of warm and cold (humid and dry) air takes place at the leading (down 
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wind) vortex interface. The tcmperaturc and humidity signals are 
well correlated. Rclatively large T values arc obtained and 

onding lower Bo values due to the existence of well defined‏ تک 
dissipative scalcs. The Bo estimatcs arc higher than over land due‏ 
to antisotropy in the small scales. There is not a ae nounced inertial‏ 
subrange in the tempcrature spectrum sincc the intermediate (transfer)‏ 
Beales are damped bysthe effects of cvaporation. Under these‏ 
excitons, as thc height above the sca surface increases, the effects‏ 
of spray evaporation will be less pronounced since spray concentration‏ 
dccreases exponentially with height.‏ 

Under warm conditions, with higher absolute humidity levels, the 
MOE structurc and associated ramp structure of the scalags begin to 
break-up. Less evaporation of spray occurs and cffects of radiative 
interaction with the tempcrature ficld duc to increased watcr vapor 
conccntrations tends to produce a morc homogeneous temperature 
structure and cause suppression of large production scalcs. Radiative 

cooling in the overlying air (inversion layer) also contributes to the 
break-up of the vortex structure and associatcd temperature and 


humidity signals. Cold spikes arc of longer duration and more spread 


22 > 8 weis 
4 0004 MAN مد‎ 


1 


But due to better mixing within the vortex structurc, This results in 
less correlation betwcen temperature and humidity (and vclocity) at 
larger scales, and a shift of production to smaller scales, Relativcly 


smaller Xo values are obtained and corresponding higher Ba values 








due to the suppression of both large scales (long wave radiati 
transfer) and diffusive scales (internal radiative transfer). ‚ppression 
of the small diffusive scales may also be due to frequency response 
degradation of the cold wire probes due to increased water droplet 
concentration. The intermediate inertial scales are more pronounced 
resulting in a well defined inertial subrange for temperature. The 
imbalance between production and dissipation of temperature variance 
can thus be attributed to effects of evaporation and combined effects of 
long and short-wave radiative transfer. 

The exact effects of evaporation and M transfer on 
BE ction and dissipation of humidity variance is not known since no 
direct dissipation measurements could be made. There appears to be 
a Slight shift to smaller scales of humidity for warm, humid 


conditions with radiative transfer. Conversely, the larger scales 

are more pronounced with relatively higher rates of evaporation 

occurring under cold, dry conditions. From measurements over land 
73 

(flat alluvial plain), Martin (1972) | suggests that at heights from _ 

2-4 m. large scale fluctuations contribute more to humidity variance 


Ethan temperature variance, and that the characteristic scale of the 


at 


۰ ene . a , i 2 e*, s ٹر‎ 本 : : یھ‎ Ree ۰ 
NN Die مر یی وت‎ e gE ER aa ml RU 


humidity field may be larger than that of temperature by about a factor 
of two. 
The model proposed above for lower levels of the marine does 


not appear to be in conflict with recent models of the marine boundary 
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ayer. Within the last few years, studies of the structure of the marine 
boundary layer have become highly sophisticated with the development 
and on of ultra sensitive radar, reliable low level aircraft 
observations, and fine scale measurements of velocity, temperature 
and humidity. From aircraft measurements of momentum, sensible 
and Pteni heat fluxes during Operation BOMEX, Grossman and Bean 
(1973), and Bean et al (1972) 2 (water vapor flux only) have suggested 
Buspibe structure of the marine boundary layer at the 16-150 mi. 

level consists of a combination of linear, buoyant, convective plumes 
or horizontal roll vortices at large and small scales, and random 


76 


cellular convective elements. Work of Konrad (1970) and Lemone 

77 ۱ ۱ 
(1972) , based on radar soundings and aircraft data, appear to 

74 

support the model of Grossman and Bean , however, Lemone 
suggests that the horizontal structures observed at wavelengths from 
1.5 to 6.5 km.(near neutral conditions) may coexist with other scales 
of motion (large and small) with increasing instability. Dean etal 
suggest that from humidity measurements (at 18-125 m) smaller scalcs 


are observed on crosswind runs, indicating more than one scale may 


be dominant in the marine boundary layer. Their humidity time series 


۰ 8 - 
۰ 
e f = 


۰ 


r exhibit the ramp (sawtooth) structure for along wind runs (along wind 
is that measured from FLIP), and a "top hat”! structure is observed on 
51 
crosswind runs. Gibson et al note that buoyant plumes may produce 


the distinctive sawtooth behavior but that they are only a sufficient and 
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^N 
not a necessary eondition for ramps to exist, 
41 . 
Lemone has suggested that the horizontal roll vortices axes 
may be between -5? and 20? to the left of the geostrophie wind. 
Er o ۹ ۱ 
Leavitt suggests ~ 10° axes orientation and a wavelength of 
“1.4 km for the roll vortiees during BOMEX. The MITOS I and OWAX 
time series (figures 7-9) suggest that different seales can eoexist within 
an individual ramp strueture, from wavelengths of :100-200 m , down 
5 
to seales comparable to those of the Gibson et al | model or smaller. 
T 
It seems reasonable to expect that scales from 1-6 km (Lemone ) 
may also eontain smaller seales. 

Results of the MITOS I and OWAX eospeetra analysis suggest 
that the dominant seale size governing the transport of heat ean vary 
aeeording to the meteorologieal eonditions. The scale sizes ean also 
vary with height as suggested by comparison of the MITOS III and OW AX 
we and u§ results for eold, dry eonditions, and by comparison of 

11 ۱ 
Leavitt's 8m and 30 m temperature speetra and sensible heat flux 
cospeetra. In both cases (not as pronouneed for the eospeetra), the 
larger seales En < 0.01) have more energy (normalized) for Leavitt's 

11 

30 m data than for the 8m. Leavitt suggests the low frequency 
enhancement may be due to radiational warming of thermoeouples; n— 
however, he also states that not all periods of daylight exhibited the low 
frequeney fluetuations. 


mil 7 
لا‎ whus suresested that air-soea temperatur& differénccs 
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and sea Sara temperature fluctuations may also contribute to 
differences between temperature and humidtiy measured ina cold, dry 
environment as compared to warm humid conditions, The air-sea 
temperature differences for MITOS I and OWAX were comparable 
suggesting that this is not a significant factor. Possible coupling 
between internal wave events, 158 frequency sea surface temperature 
fluctuations, and air temperature fluctuations, S may be 
affecting the temperature structure in at low levels in the atmospheric 
surface layer. 3 

Measurements of temperature over land as compared to results 
Over the ocean indicate that the positive buoyancy generated by ز‎ 
differences (enhanced by sensible heat transfer) can sustain a vortex 
structure (ramp structure) much better when only evaporation occurs. 
Boston and Beine result of Bo = 0.8 over a semi-moist flat 
suggests that additional buoyancy of water vapor and corresponding 
evaporation may account for the higher Bo value due to the anisotropy 
of the temperature field with a ramp structure. The positive buoyancy 
^ and lack of evaporation 5 dry land te Kansas) may not en 


to sustain roll vortices when coupled with more severe surface 


5 C ۰ e 
۰ do: x à +, 
.'. 


5 ۳ emer e . = - P 25 
roughness condition, although Haugen etal have suggested that 
buoyant plumes can effect the heat transport process over land. 
In conclusion it appears that the temperature spectra (and Bo 


values) and sensible heat flux cospectra do not follow universal 
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similarity over the ocean. The humidity spectra and cospectra may 
require modified scaling than what is normally used now to follow 
similarity predictions and the influence of evaporation and radiative 
transfer on the spectra and cospectra is not clear. Differences 
between production and dissipation of — varıance, and the 
behavior of the temperature spectra and sensible heat flux cospectra 
could be accounted for by consideration of other sources of production 
usually neglected in the temperature variance.budget, and attributed 
to the combined effects of ocean spray evaporation and radiative 
heating/cooling in light of the model proposed for the lower levels of 
the UB oe boundary layer. 
5,5 Recommendations for Further Study 

Application of the dissipation technique for estimation of sensible 
and latent heat and the bulk aerodynamic technique for estimation of 
sensible heat requires further study. A quantitative evaluation of the 
influence of evaporation of ocean spray on both the temperature and 
humidity variance budgets and the spectra and cospectra is desirable, 
This bout include mu of the laboratory work onde éted by ur 


to field experiments over the ocean to determine droplet concentrations, 


e un g’ 


NU, 


variations with height, and rate of evaporation for various 
meteorological conditions. The possible problem of water droplet 
discrimination of the Lyman-alpha sensor should be investigated. This 


could be done in conjunction with simultaneous laser spray 








measurements. A dynamic frequency response calibration of the cold 
wire probes is needed; this work currently is underway in the 
research group at UCSD. 

The influence of sea surface temperature fluctuations and 
coupling between internal wave events, sea surface temperature 
fluctuations, and air temperature fluctuations also requires further 
study. The temperature cold spikes and ramp structure may also be 
associated with the above phenomena, and also correlated with surface 
waves. a 


The influence of radiation on both the ern and humidity 
melds neo require further study. A promising laboratory investigation 
is currently underway at the large water-wave facility at the Institut 
de Mechanique Statistique de la Turbulence, in Marseille, France, 
(Coantic and Favre (1970) °°, (1973) ^) to determine the effects of 
radiation, evaporation, waier waves, and stability conditions on the 
transfer phenomenon under controlled conditions. 

Application of various models of the marine boundary layer for 
the MT Rrahnding of T M ees sen occurring at the interface 
is indispensable. It appears that various proposed models may not 
y Ee n a A that — CES set Ge CUP TE 


simultaneously may account for the various model types. 
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Instantaneous corrections to the turbulent fluxes based on all 
components of motion relative to an inertial reference frame require 
further study. This work is currently underway at UCSD using outputs 
of a Litton inertial platform (see figure 3) and techniques similar to 

۱ 089 5 
those employed by Mitsuta and Fujitani. Recent work by Friche 

Dr t m 
(1974) | indicates that fluxes determined by the direct covariance 
technique using a sonic anemomceter to measure velocity, may be 
affected by contamination of the velocity measurements by temperature 
fluctuations. A major effect is on the shape of the cospectra, 


especially the horizontal heat flux cospectra. Effects of this 


contamination require further study. 
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NORMALIZED CORRELATION COEFFICIENTS 


Correlation coefficients presented and discussed in Section 4. 4 
were normalized to remove any trends due to possible stability 
dependence of the correlations as suggested by Kaimal and Haugen. 

The normalization may be derived from normalized variances 
and covariances by integration of the normalized spectra and cospectra 


of Sections 4.5 and 4, 6 and applying the relation (from Eq, (2) and 


Sections 4.5 and 4. 6) 


2 
(c) = × -f © (f) d£ =f DE (A.1) 
and 
= & 1( ص۳پ‎ 3 df A2 
xy = fb, (fy di = (Oy) af, (A2) 
Normalized correlation coefficients are then defined as 
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and similarly for CI . Using the expressions for the normalized 
a | 
forms of the spectra and cospectra of Sections 4. 5 and 4, 6 the 


normalized correlation coefficients are 
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